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INTRODUCTION 


The  problem  to  be  considered  is  the  analysis  of  rectangu¬ 
lar,  laminated,  fiber-reinforced  plates.  The  loading  conditions 
are  assumed  to  consist  of  normal  pressure  applied  to  one  lateral 
face ,  arbitrary  edge  forces  and  moments  applied  to  the  edges 
and  a  uniform  temperature  distribution.  Internal  boundary 
conditions  involve  the  continuity  of  the  three  components  of 
the  displacement  vector  and  the  components  of  the  stress  vector 
normal  and  tangent  to  the  individual  plies . 

Assuming  that  the  primary  interaction  is  between  plies, 
the  plate  is  modeled  as  a  laminated  plate  for  which  each  ply 
is  assumed  to  obey  a  stress-strain  law  with  twelve  independent 
elastic  constants  and  three  independent  coefficients  of  ther¬ 
mal  expansion.  These  local  properties  are  assumed  to  vary  dis- 
continuously  from  ply  to  ply.  Internal  continuity  conditions 
at  ply  interfaces  can  be  identically  satisfied  by  an  integral 
formulation  of  three-dimensional  stress  function  theory. 
(References  1  and  2.) 
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SECTION  I 


FORMULATION  OF  THE  PROBLEM 


Let  {x,y)  be  the  in-plane  coordinates  and  z  the  normal 
coordinate  to  the  mid-plane  of  the  laminated  plate.  It  will 
be  assumed  that  (1)  individual  laminations  are  homogeneous 
and  elastic,  (2)  lamination  interfaces  are  perfect,  i.e.,  the 
displacement  vector  and  the  normal  and  tangential  components 
of  the  stress  vector  are  continuous  across  all  interfaces  and 
(3)  all  dependent  variables  are  continuous  in  (x,y)  and  dif¬ 
ferentiable  to  any  order  required.  The  nature  of  the  con¬ 
tinuity  of  dependent  variables  with  respect  to  the  z  coordinate 
will  be  specified  by  stating  the  appropriate  function  class 
C^(z) ,  n  =  0,1.  Functions  which  are  discontinuous  in  z  will 
be  marked  with  an  asterisk.  Let  at.(z),  a'^iz)  be  the  elastic 
functions  characterizing  the  kind  and  distribution  of  material 
within  the  plate,  -h  <  z  <  h. 

A  sufficient  condition  for  the  continuity  of  the  dis¬ 
placement  vector  within  the  plate  (u^€c®(z))  is 

z 

Uj^  =  I  u|(x,y,t)dt  +  u^(x,y,-h)  (1) 

-h 

The  strain-displacement  relations  are 


+  u.  .) 
3  /I 


where  a  comma  denotes  partial  differentiation  with  respect  to 
=  X,  ^2  =  y  or  x^  =  z.  Substituting  equation  (1)  into 
equation  (2) ,  three  of  these  equations  can  then  be  solved  to 
give  the  displacement  functions  u|  in  terms  of  the  strains. 
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u?  =  Y  „ 

1  '  xz 


(x,y,t)dt  -  1^  U3(x,y,-h) 


uJ  =  Y 
2  ’yz 


(x,y,t)dt  -  ■—  U2(x,y,-h) 


ui  =  Y 

3  '  zz 


The  remaining  three  equations  can  then  be  manipulated  to  produce 


three  compatibility  equations. 


"''xx  ^  J  (x,y,t)dt  -  J  (z-t)  -—■I--  (x,y,t)dt 

T_  ""h 


-  (z+h)  U2(x,y,-h)  +  Uj^(x,y,-h) 


r  dy  f  3^Y 

^  J  (x,y,t)dt  -  J  (z-t)  (x,y,t)dt 


(z+h)  —J  U2(x,y,-h)  +  U2(x,y,-h) 


_  f  ^yz  ^xz 

^xy  J  3x  3y 

-h 


^  i 

dt  -  I 


(z-t)  ■■  (x,y,t)dt 


"  +  2  ^  U2(x,y,-h)  +  ^  U3L(x,y,-h) 


Three  other  compatibility  equations  can  be  derived  from  equations 
(4)  which  are  obviously  not  independent  of  equations  (4)  and 
which  are  satisfied  when  equations  (4)  are  satisfied. 

The  equilibrium  equations  are. 
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9t  9t  9t 

XX  .  xy  ,  xz  „  ^ 


!ls:  +  !iiL  +  !Iz£  =  0 

9x  9y  9z 


9t  9t  9t 

xz  yz  zz 

9x  9y  9z 


These  equations  are  identically  satisfied  by  a  representation 
in  terms  of  six  stress  functions  . 


*^22,33  ‘*’33,22  "  ^‘*^23,23 

‘*’33,11  ‘*’11,33  "  ^‘*’13,13 

‘*’11,22  ‘*’22,11  "  ^‘*’12,12 


X  s=  -<{) 

xy 


33,12 

^‘*’13,2 

+ 

‘*’23,1 

‘*’12,3^ 

,3 

22,13 

+ 

(‘*’23,1 

+ 

‘*’12,3  ■ 

‘*’13,2^ 

,2 

11,23 

+ 

(‘*’12,3 

+ 

‘*’13,2 

‘^23,1^ 

,1 

(6) 

The  stress  components  t  and  t  must  be  continuous  in  z. 

xz  yz 

A  sufficient  condition  for  the  satisfaction  of  these  continuity 
conditions  is 

z 

4)..  =  (j)?.(x,y)  +  (z+h)  (|)}  .  (x,y)  +  f  (z-t)  <()t .  (x,y  ,t)  dt  i,j=l,2 

1313  ^3  3  “^3 


'•’13  =  + 


I  't'|3(x,y, 


t)dt  i=l,2 


'33  '"33 
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The  stress  function  4123  may  be  discontinuous  since  it  does  not 

appear  in  t  ,  t  or  x  .  Using  equations  (7) ,  the  stresses 
xz  z  z 

become 


''xx  ■  *22  *33,22  ■  ^*23,2 


*33,11  *  *11  -  2*13,1 


z 

I  (z-t)  ‘*’22,11  "  ^‘*’12,12] 

^  J  x,y /t: 

-(z+h)[|j  T^^(x,y,-h)  +  ly  'yz'’''''''**'] 


33,12  ‘*’13,2  ‘*’23,1  “  ‘*’12 


z 

^xz  “  I  ['*’23,12  ^  '*’12,2  "  ‘*’13,22  "  ‘*’22,1] 

z 

^yz  “  j  [‘*’12,1  '*’13,12  "  ‘*’23,11  "  '*’11,2] 

-h  x,y,t 

(8) 

The  functions  (|)?j  and  ()»^j  have  already  been  determined  by  evalu¬ 
ating  the  stresses  at  z  =  -h.  Assuming  that  this  face  of  the 
plate  is  free  of  applied  stress,  we  can  delete  x  (x,y,-h) , 

XZ 

x„^(x,y,-h)  and  x_„(x,y,-h)  from  equations  (8), 
y  z  z  z 

Hooke's  law  for  an  individual  ply,  referred  to  the  princi¬ 
pal  material  axes  (y^^,  Y2>  is 
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’'^11  "  ^11^11  ®12'^22  ^13'^33  ‘^l'^ 


^22  “  ®2l’^ll  ^22''^22  "'■  ^23'^33  ■‘’  ‘^2^^ 


^33  ^3l’^ll  '''  ^32^22  ■*■  ^33'^33  ■‘'  ^3^^ 


^12  ”  ^44'^12 


^13  ^55'^13 


^23  ”  ^66^23 


In  order  to  transform  equations  (9)  to  the  (x,y,z)  coordinate 
system,  the  coordinate  system  must  be  rotated  through  an  angle 
6  in  the  x,y  plane.  The  transformation  laws  for  the  stress  and 
strain  tensors  are 


Y . .  =  k .  k . oY  „ 

'  i;]  la  ]  3  '  a3 


T7  .  =  k.  k . oT  o 

ij  la  33  a3 

where  k^ ^  “  ®i  *  ® j ’  primed  components  refer  to  the 

(x,y,z)  coordinates  and  the  unprimed  components  refer  to  the 
principal  material  axes  (y^/  Y2'  ^3^  *  direction  cosines 

of  the  transformation  are 


cos9  -sin0  0 


COS0  0 


By  carrying  out  the  transformation,  substituting  the  results 
into  equation  (9)  and  solving  for  the  strains,  the  stress-strain 
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law  referred  to  the  (x,y,z)  axes  is  obtained. 


Y  =  bn  T  +  b^-T  +  bi.T  +  6nT 

'xx  11  XX  12  yy  13  zz  14  xy  1 

Y  =  b_,T  +  b__T  +  b^-T  +  b_.T  +  3oT 

'yy  21  XX  22  yy  23  zz  24  xy  2 

Y  =b-,T  +b-.,T  +b-,T  +b,-T  +3_T 

'zz  31  XX  32  yy  33  zz  34  xy  3 

^xy  ~  ^41^xx  ^  ^42^yy  ^  ^43^zz  ^  ^44^xy  ^  ^4^ 

^xz  ~  ^55^xz  ^  ^56^yz 

Y  =  b^j-f  b^,T  (10) 

'yz  65  xz  66  yz 

The  transformed  material  constants  are  given  in  terms  of  the 
constants  referred  to  the  principal  material  axes  by  the  fol¬ 
lowing  relations : 

4  4  2  i2 

b^^  =  a^^^cos^e  +  a22sin^e  +  (a^2  ^21  ^  6cos  6 

4  4  2  2 

^12  ~  ®  ^21®^^  ®  "**  ^^11  ^  ^22  ”  ^^44^®^^  6cos  9 

2  2 

bi3  =  0  a23sin  6 

3  3 

b^^  =  ^^^12  ~  ^11  ^  a^^)sin0cos  0  +  2(a22  ”  ^21  ~  ^44^  0cos 

^21  “  a^2®i'^^®  ■*■  a2-|^cos^0  +  (a^^^  +  a22  "  2a^^)  sin^6cos^0 

4  4  2  2 

^22  “  a^j^sin  9  +  a22COS  0  +  (a^2  ^21  ^^44^®^^  Scos  0 

2  2 

^23  ~  ®13®^^  ®  a23COS  0 

3  3 

b24  =  2(a22  “  a2j^  -  a^^)sin9cos  0  +  2(aj^2  "  ^44)  ®i^  9cos 


2  2 
a^^^cos  9  +  a32ain  0 


more 


2  2 

^32  ~  ®  ® 


*^33  ^33 


2(3^2  ”  a2^)sin0cose 


(321  -  3ii)sin9cos0  +  (322  *•  ai2)  sin'^0cos0 


+  3^^sin0cos0 (1  -  2sin  0) 


3  3 

(322  “  3i2)sin0cos  9  +  (a2i  “  3ii)sin  0cos0 


-  3^^sin0cos0 (1  -  2sin  0) 


“  ^^23  ~  ai2)sin0cos9 

2  2  2  2 
^22  ”  ^12  ”  a2i)sin  0cos  0  +  a^^(l  -  2sin  0) 

2  2 

bcc  =  3ccCOS  0  +  3^^sin  0 

55  55  66 


(3,^  -  3t.c)  sin9cos6 
66  55 


2  2 
3-.cSin  0  +  3^, cos  ( 
55  66 


2  2 
UiCos  0  +  a2sin  9 


2  2 
ttiSin  0  +  a2Cos  0 


®3  '  “3 


^4  ~  ^^^2  ~  ai)sin0cos0 


Equstions  (10)  m3y  be  extended  to  spply  to  the  entire 
l3min3ted  plete  by  ellowing  the  meteriel  constents  bi ^ ,  $1  to 
be  discontinuous  functions  of  the  z  coordinete.  This  discontinuous 
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dependence  on  z  may  by  accomplished  by  allowing  9  to  be  a 
discontinuous  function  of  z  or  by  considering  that  the  con¬ 
stants  discontinuous  functions  of  z  or  allowing 

both  possibilities  simultaneously.  The  functional  dependence 
of  the  material  properties  will  be  indicated  by  marking  the 
constants  b . . ,  0.  with  an  asterisk;  b..  =  b?.(z),  0.  =  0*(z). 

By  substituting  equations  (8)  into  equations  (10) ,  the 
strains  may  be  obtained  in  terms  on  the  stress  functions  (()|j . 

It  is  clear,  however,  from  an  examination  of  equations  (4) 
that  the  strain  components  "*^yy'  continuous  func¬ 

tions  of  z.  The  remaining  strain  components  may  be  discontinu¬ 
ous.  Therefore,  we  shall  examine  the  strain  components  y  t 


^yy '  ^xy ' 


^11  [^22  ‘*’33,22  ■  ^‘*’23,2]  ^12  [‘*’33,11  ‘*’ 


*  -  2d)* 

11  "'^13,1 


z 

+  b*  J  22  *  *22,11  ■  2*12,12] 

+  44[-  *33,12  +  *13,2  +  *23,1  '  *12]  *  ^l'' 


‘’2l[*22  *33,22  "  ^*23,2]*  '^22^*33,11  *11 

z 

+  b*3  I  [‘J>11,22  ‘*’22,11  "  ^‘*’12,12] 

^  “24  [-  *33,12  +  *!3,2  ^  *23,1  -  *12]  +  er 


2*13,1 


more 
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'xy 


* 

11 


^4l[‘**22  '^33,22  “  ^‘*’23,2]  ^42  [*^33,11 

z 

I  (z-t)  [(|)*^^22  '^•‘^22,11  "  ^‘*^12,12] 

+  >=44  [-  *33,12  +  *13,2  +  *23,1  "  *12]  + 


^*13,1] 


+  b*3  I  (z-t) 
-h 


(12) 


If  these  strain  components  are  to  be  continuous,  as  they 
must  be  in  order  to  be  consistent  with  equations  (4) ,  then  we 
must  seek  some  decomposition  of  the  stress  functions  (})|j  in 
terms  of  b^ ^ ,  3|  and  a  new  set  of  stress  functions  which  are 
of  class  C^(z).  Let  us  first  eliminate  the  discontinuous 
temperature  terms  by  making  the  transformation 


4*1  =  ♦*!  4-  P*iT 


'*'22 


^12 


tJo  +  F*,,T 


'22 

I* 

12 


ip*^  +  F*^T 


22- 

?*  r 
12‘ 


(13a) 


and  choosing  FjjL»  ^22'  ^12  eliminate  the  temperature, 

The  result  is 


D*F* 

u 


D*P* 

^  ^22 


3* 

b* 

'^ll 

44 

b* 

°21 

^24 

^4 

>=41 

44 

"^12 

8£ 

b* 

14 

^22 

^2 

b* 

24 

CM 

Xi 

®4 

44 

(13b) 

more 
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b* 

12 

b* 

11 

H 

D*PJ2  =  - 

b* 

22 

b* 

21 

b* 

^"42 

b* 

^41 

b* 

12 

b* 

°11 

*^14 

D»  =  - 

b* 

^"22 

*21 

*24 

b* 

^"42 

Hi 

*l4 

The  expressions  for  the  strains  then  become 

■^xx  ‘*’33,22  "  ^‘*’23,2]  ^12  [‘*’33,11  ^  ’*’11 

z 

+  b*3  I  (z-t)  +  *^2^33  -  2*J2,12] 

— X  f  y  /  u 
+  bl4[-  *53, 12  *  *13,2  ■  ^h,l  ■  ''’12] 

^yy  ■  *’2l['*’22  ■*■  *33,22  '  ^*23,2]  *  *^22  [*33,11  ■*■  *11 
z 

*^23  /  [‘*’11,22  ■*■  ‘*’22,11  “  ^‘*’12,12] 

x,y,r 

+  *"24  ["  ‘*’33,12  +  ‘*’13,2  +  ‘*’23,1  ■  ‘*’12] 

^xy  *"4l[’*'22  ‘*’33,22  “  ^‘*’23,2]  *"42[‘*’33,11  ’*’11 

z 

+  b*3  I  (z-t)[»*3_33  -f  **3^33  -  2**3_33]  dt 

_h  x^y/t. 

*«[■  *33,12  *13,2  *23,1  ‘  *12] 


(13b) 

2*13,1] 


2*13,1] 


2*13,1] 
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Now ,  let 


>^^11 

=  =!Ai  + 

^12^22 

+ 

CJ3A12  + 

1 - 1 

! 

n2 

=  + 

^22^22 

+ 

=23*12 

^22 

^2 

=  - 

C*2A 

22 

-  =33*12 

^12 

(14) 

where  the  new  stress  functions  are  assumed  to  be  of  class 
C^(z)  and  the  discontinuous  functions  are  to  be  chosen  so 
that 


XX 

^22 

yy 

^  -^11 

=  A,  « 

(15) 

xy 

12 

This  leads  to  a  set  of  nine  equations  for  the  discontinuous 
functions  C? . . 


b*  C* 
^11^21 


^12^11 


4-  h  *  C  * 

°14^31 


+  ‘>i4=?2 


'=!2=13 


*^14^33 


b*  C* 
21  21 


+  ‘’22^=11 


+  b*  C* 
^24^31 


b*  C* 
^21^22 


^22^12 


^24^32 


b*  C* 
^21  23 


+  ‘>22=13 


+  b *  P* 
24  33 


b*  C* 
^41^21 


+  b*  C* 
42  11 


+  »44=31 


b*  C* 
^41  22 


4-  b*  C* 

^42  12 


*^44^32 


^ll^23 


+  ^^4243 


+  '^44''33 


0 

1 

0 

1 

0 

0 

0 

0 

1 
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These  equations  can  be  broken  down  into  three  independent  sets 

involving  •  ^^13 '^23 '^33^  respectively. 

All  three  sets  have  the  same  determinant 


A*  = 


42 

b* 

Dll 

44 

42 

b* 

^21 

44 

4*2 

b* 

°41 

44 

and  the  functions  C|j  are  found  to  be 


A*C*j_ 


b*  b* 
“l4°41 


b*  b* 
°11°44 


A*C* 

^  ^21 


b*  b* 
12^^44 


b*  b* 
“l4°42 


i*C‘^ 


b*  b* 
11^42 


b  *  b  ^ 
^12^41 


A*r* 

12 


b*  b* 
21^44 


b*  b* 
°24°41 


A*C* 

^  ^22 


^24^42 


^22^44 


A*C* 

^  ^32 


b*  b* 
22^41 


‘=2l‘’42 


A*C* 

^  ^13 


‘>!4‘=21 


A*C*3 


b*  b* 
°22“l4 


b*  b* 
12^24 


A*C*3 


^12^21 


^11^22 


(16a) 


(16b) 


The  stress  functions  are  determined  by  eliminating  all 
remaining  discontinuous  terms  from  the  expressions  for  y  , 
Yyy»  Before  writing  out  the  expressions  for  these  stress 

functions  however,  it  may  be  noted  that  th^re  are  only  three 
compatibility  equations  to  be  satisfied,  equations  (4) .  There 
are  already  three  continuous  stress  functions  available; 

^22'  ^12 ‘  analysis  of  the  isotropic,  laminated  plate 
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(Ref.  1)  it  was  found  that  the  three  stress  functions  4)^2 
disappeared  from  the  final  expressions  for  the  stresses  and 
could  therefore  be  taken  as  identically  zero.  It  will  tenta¬ 
tively  be  assumed  that  these  three  stress  functions  will  not 
be  needed  for  the  present  analysis.  By  taking  these  three 
stress  functions  as  zero,  the  expressions  for  are  greatly 

simplified. 


^11 


(b*3C*2  +  +  t>|3C|3)A33 


^22  “  “  ^^13^22  ^23^21  ■*"  ^43^23^^33 


Q*  = 

“12 


(*^13^32  ^23^31  ^43^33^^33 


(16c) 


For  convenience,  a  new  continuous  stress  function 
z 

^33  "  [  ^^“^^[‘*’11,22  ‘*'22,11  ■  ^‘*'12,12]  " 

(z-t)  [c*j^A^3^^22  ^12^22,22  ^13^12,22 


* 


^21^11,11  ^22^22,11  ^23^2,11  ^^3Ai,12 

+  2052^22,12  2C*3A^2,12  -  H3^*12  +  ^2*3^^  ^13^13) 

^33,22  ■  ^*^13^22  ^23^21  ^43^23^^33,11 

"  2  (*^*3^32  ■*■  ^23^31  ^43^33^^33,12  ^il'^'22 


^  /  Y  /  ^ 


(17) 


has  been  defined.  It  should  be  noted  that  equation  (17)  is  an 
integro-differential  equation  involving  all  four  continuous 
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stress  functions  ^22'  '^12'  ^33  temperature,  through 

the  terms  ^22  equation  (13).  The  results  of  this  de¬ 

composition  may  be  summarized  by  writing  out  the  final  expres¬ 
sions  for  the  stress  functions  (}>| j  /  in  terms  of  ^^22’  ^12' 

A22  and  the  temperature,  and  the  corresponding  expressions  for 
the  stresses  and  strains. 


ni 

=12*22  + 

^13^12 

*  *^23=11 

+  ^43^3 

^2 

=22*22 

^23^12 

- fb*  C* 
^^13^22 

+  ‘=23=21 

^43^23 

<*’12 

=  - 

=32*22 

^33^12 

.(b*3C*2 

‘=23=31 

^  ^43''33 

12- 


<0*3  5  0 


i  =  1,2,3 


(18a) 


XX 


<05  _ 

^22 


yy 


'^11 


zz 


xy 


xz 


yz 


=  A 


33 

-  <0* 


d. 

-h 

: 

-\ 

-h 


12 

,<*’12,2  " 

<^2*2, 1^ 

dt 

x,y,t 

‘^12,1  " 

^1,2! 

dt 

x,y  ,t 

(18b) 
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^^31^21  ^32^11  *^34*^31^^11 

^*^31*'22  *^32^12  ’*’  *^34*^32^^22 

^*^31^23  *^32^13  ■*■  *^34*^33^^12 

I 

-  [‘>31<'=13'=22  +'‘^23'=21  +  ‘’43<'23>  ’  >=33 
+  ^32  (*^13^12  *  >=23‘»U  '*^  '=43^13> 

+  b34(t'f3C32  ■•■  *^23>'31  '*'  '=43^33>]*33 
+  (b|^F*2  +  b*2F*^  -  b*^F*2  +  e|)T 


f  [- 


^31^11,2  ^32^22,2  ^33^12,2 


^*^13^32  *^23*^31  '''  *^43*^33^^33,2 

^2Ai,1  "  ^22^22,1  '■  *^23^12,1 


+  (b*3C*2  +  b*3C*3  +  b|3C*3)A33^3  +  F*2T^2 


^22^ 


+  *^§6 


^31^11,1  ^32^22,1  ^33^12,1 


+  (b*3C*2  +  b*3C*^  +  ^|3'*33^''33,1 
^11^11,2  "  *^12^22,2  ■  *^13^^12,2 


■*■  ^*^13*^12  *^23*^11  '''  *^43*'13^^33,2  ^12"*^, 1 


F*iT 


*^65  J  [“  ^3Ai,2  "  *^32^22,2  *^33^12,2 


+  (b*3C*2  +  b*3C*^  +  *^|3''33*''33,2 


more 


^21^11,1  ^22^22,1  ■  ^23^2,1 

^^13^22  ■*■  ^23^21  ^43^23^^33,1  ^12^,2  "  ^22^, iJ 


■*■  ^66  I  ["  ^3Ai,1  "  ^32^22,1  "  ^33^12, 
-h 


1 


+  r*  +  h*  r*  4-  r* 

^^13^32  ^  ^23^31  ^  ^43^33^^33,1 

_f*a  _p*a  _p*a 

^11^11,2  ^12^22,2  ^13^12,2 

+  /h*  P*  +  h*  P*  +  h*  P*  ^  A  +  P*  T  -  P*  T 

+  ^  ^23^11  ^  “43^13^  ^33,2  ^  ^12^,1  ^11^,2 

(19) 

At  this  point,  it  may  be  proven  without  much  difficulty, 
using  equations  (4),  (17)  and  (19),  that  A22  /  A3^2€'C^(z) 

and  A33€C^(z)  and  the  formulation  of,  the  static  problem  is  thus 
complete. 

To  obtain  the  governing  equations,  it  is  necessary  to 
substitute  the  expressions  for  the  strains ,  given  by  equation 
(19)  into  the  compatibility  equations,  equations  (4)  and  (17). 
The  expressions  for  the  strains  can  be  simplified  somewhat  if 
several  groups  of  terms  which  appear  frequently  are  replaced 
by  single  symbols . 


"1 


b*  C*  + 

“13^12 


b*  C* 
23  11 


+  b*  P* 
^43^13 


^2* 


K*  P*  4 

“13^22 


b*  C* 
°23^2l 


+  b*3C*3 


"3 


^13^32 


^23^^ 


"^43^33 


"4 


K*  p*  + 

^31^21 


b*  C* 
^32^11 


+  b*  C* 
^34^31 


"5 


^31^22 


b*  C* 
°32^12 


^34^32 


^6 


^31^23 


^32^13 


^34^33 


17 


dt 

x,y,t 


] 


■o  >1 


The  strains  can  then  be  written  in  the  form: 


XX 

"  ^22 

yy 

= 

xy 

^12 

zz 

11 

H 

> 

'xz 


^55  I  ["  ^31^11,2  "  ^32^22,2  “  ^33^12,2  “  ^2Al,l 
-h 

^22^22,1  "  ^23^12,1  ^3^3,2  ^2^33,1  ^  ^12^^, 2 


^22\l] 

x,y,t 

z 

^56  f  ["  ^31^11,1  "  ^32^22,1  “  ^33^12,1  "  '^11^11,: 


+ 


-h 


^12^22,2  "  ^13^12,2  ^3^33,1  ^1^33,2  ^12'^,! 


_  F* 
^11 


T  ]  dt 
'  ■'x,y,t 


Y. 


yz 


^65  \  [”  ^3Ai,2  ”  ^32^22,2  “  ^33^12,2  “ 

-h 

^22^22,1  "  ^23^12,1  ^3^33,2  ■*■  ^2^33,1  ^12'^,  2 

“  ^22^  l] 
z 

■''  ^66  I  ["  ^sAl,!  "  ^32^22,1  “  ^33^12,1  " 


-h 


(20) 

more 
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°12*22,2  ■  ^13^2, 2  ''3*33,1  ''1*33,2  ^12'',1 


-  ^ii 


^2]  ■'% 

x,y,t 


(20) 


Using  equations  (18)  and  (20)  ,  the  inverted  stress-strain  law 
may  be  written  in  the  form 

■  fe  •  ... .  fe. .  m 


*  (''23  +  -ir> 


'zz  r 22  I*  J 


xy  I ^  ' z  z 


'yy  =  ^2  ^  '^KX  +  (''il 


-w  'yy 


(=13  *  -f*-)'' 


-^1,  .fy. 

7  y  xy  ^zz  \^^11 


^xy  1^32  I  *  J  ^xx 


(=5i  ^ 


'yy 


/ 

(^33  T*~l 


-I  J  xy  "  ly  ^zz  "  (jl2  T^J 


zz 


1 

jT 

7 


4  '  y  y  6  ^  xy 


fzz  *  'r' 


-  bfcY 


xz  -  b|5b|6  -  bggb^j 


yz 


“^65^xz  ^55^yz 

-  ^*5eH5 


(21) 


The  quantities  C|j  are  the  moduli  of  the  individual  plies  only 
under  the  classical  plate  theory  assumption  that  the  trans¬ 
verse  normal  stress  is  negligible.  Otherwise,  the  correct 
expressions  are  those  shown  above. 
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The  governing  integro-dif ferential  equations  may  be 
obtained  by  substituting  the  expressions  for  the  strains 
into  equations  (4)  and  (17).  However,  this  will  not  be 
necessary  at  this  point;  it  is  only  necessary  to  note  that 
these  equations  involve  integrals  of  discontinuous  quantities 
with  respect  to  the  z  coordinate.  All  of  the  terns  on  both 
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SECTION  II 


THE  LIMITING  CASE  OF  AN  INFINITE  NUMBER  OF  LAMINATIONS 

The  integrals  involved  in  the  compatibility  equations  are 
of  two  types : 

z 

I  F*(t)f(t)dt 

-h 

z  z 

I  G*(t)  I  H*(s)f  (s)dsdt 
-h  -h 

where  F*,  G*,  H*  represent  functions  which  are  piece-wise 
constant  in  z  and  f  represents  functions  which  are  continuous 
in  z  (of  at  least  class  C^)  with  a  bounded  derivative  every¬ 
where  in  the  interval  of  integration.  The  function  f  is  in¬ 
tended  to  represent  any  of  the  functions  and  F*,  G*,  H* 

represent  combinations  of  the  material  property  functions  as 
they  appear  in  the  compatibility  equations.  If  the  laminated 
plate  consists  of  an  infinite  number  of  laminations  arranged 
periodically  by  repetition  of  a  lamination  subgroup  containing 
a  finite  number  of  laminations,  then  it  has  been  shown  (Ref¬ 
erence  2)  that  the  integrals  have  the  limiting  values: 

z  z 

Lim  I  F*(t)f(t)dt  =  F  J  f(t)dt 
-h  -h 

z  z  z 

ni~  [  \  H*(s)  f  (s)dsdt  =  G  H  I  (z-t)f(t)dt  (22) 

-h  -h  -h 

The  average  values  are  defined  as  arithmetic  means  taken  over 
the  thickness  of  the  plate.  By  applying  this  theorem  successively 
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to  all  terms  contained  in  the  governing  integro-dif ferential 


equations  for  the  stress  functions  '^22'  ^12'  ^33' 
be  shown  that  these  equations  are  reduced  to  a  form  in  which 
all  material  property  functions  are  transformed  into  constants 
which  are  defined  as  mean  values  of  those  property  functions. 
The  case  in  which  the  nximber  of  laminations  approaches  infinity 
must  therefore  correspond  to  a  homogeneous ,  anisotropic  plate 
since  the  integro-dif ferential  equations  for  the  homogeneous 
plate  have  constant  coefficients.  However,  it  is  not  known  at 
this  point  what  the  material  properties  of  the  equivalent, 
homogeneous  plate  are.  In  order  to  determine  this,  it  will 
be  necessary  to  develop  the  integro-dif ferential  equations  for 
a  sufficiently  general,  homogeneous,  anisotropic  plate  with  un¬ 
determined  material  constants  and  then  to  evaluate  these  con¬ 


stants  by  matching  coefficients  between  the  two  sets  of  integro- 
dif  ferential  equations.  The  form  of  the  stress-strain  law  for 
the  equivalent,  homogeneous  plate  must  be  determined  by  trial 
unless  the  most  general  form  possible  is  assumed.  The  latter 
possibility  is  too  tedious  to  evaluate  but  it  seems  reasonable 
to  expect  that  a  laminated  plate  in  which  the  individual  plies 
obey  the  stress-strain  law  given  by  equation  (10)  will  not  pos¬ 
sess  any  higher  degree  of  anisotropy  than  the  individual  plies. 
If  all  coefficients  in  the  integro-dif ferential  equations  for 
the  two  cases  can  be  matched  without  any  inconsistencies,  then 
this  assumption  will  have  been  proven  correct.  Accordingly, 
the  stress-strain  law  for  the  equivalent,  homogeneous  plate 


will  be  taken  in  the  form: 
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The  integro-dif ferential  equations  for  a  plate  composed  of  this 
material  must  now  be  derived  by  following  formally  the  proce¬ 
dures  carried  out  in  the  previous  section:  FORMULATION  OF  THE 
PROBLEM.  It  is  necessary  that  this  derivation  be  formal  since 
the  enforcement  of  internal  continuity  is  not  now  required, 
and  the  motivation  behind  the  development  of  these  equations  is 
thus  obscured.  It  will  be  clear,  however,  that  all  of  the 


equations  thus  obtained  will  be  identical  to  those  previously 

derived  with  b|j  replaced  by  d^^  and  3*  replaced  by  e^.  The 

quantities  which  are  defined  by  equations  (16b)  are  re- 

H 

placed  by  which  are  defined  by  equations  (16b)  with  b^^ 

replaced  by  •  Also,  the  quantities  I^  will  be  replaced 
H 

by  I^,  defined  in  a  similar  way.  We  are  now  in  a  position  to 
compare  the  integro-dif ferential  equations  for  the  homogeneous, 
anisotropic  case  and  the  limiting  case  of  the  laminated  plate 


as  the  number  of  laminations  approaches  infinity.  By  equating 
all  coefficients  of  the  two  sets  of  equations,  the  following 
set  of  relations  is  obtained.  From  equation  (17) : 


23 


From  equation  (4b) 


d 

^31  ^^65 


^31  ^65 


d 

^21  ^65 


=21  ‘’65 


4  <‘65 


^32  ^65 


c“  d 
^22  ^65 


^22  ^65 


pH  H 
^33  °65 


^33  ^65 


=23  ^65 


C  b 
^23  65 


^3  *^65 


^3  ^65 


^2  ^65 


^2  ^65 


^12  ^65 


F  b 
12  °65 


^22  ‘^es 


F  b 
22  °65 


d 

^31  °66 


^31  ^66 


c“  d 
^11  ^66 


^11  ^66 


d 

^32  °66 


^32  ^66 


pH  H 
^12  ^66 


^12  ^66 


H  =  P  Vi 

^33  °66  ^33  °66 


pH  d 
''13  °66 


^13  ^66 


^3  ‘^ee  ^3  ^66 


tH  ^ 

^1  '^ee 


h  ^66 


^12  ^66  ^12  ^66 


^11  '^ee 


^11  ^66 


The  remaining  relations  for  this  equation  are  the  same  as  the 
last  five  from  equation  (4a) .  Similar  relations  for  equation 


(4c)  only  duplicate  those  already  found.  The  distributive 
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property  of  the  arithmetic  mean  has  been  used  in  some  of  the 
above  relations.  Asterisks  have  been  omitted  from  the  quantities 


bi  j ,  and  1^^.  All  such  quantities  appear  in  terms  under 

a  bar  which  designates  the  arithmetic  mean  taken  over  the  z 
coordinate. 

The  task  now  at  hand  is  to  solve  the  foregoing  set  of 
relations  for  the  twenty  equivalent  elastic  constants  d. . 

1  J 

and  the  four  equivalent  coefficients  of  thermal  expansion  e^. 
Considerable  simplification  may  be  obtained  by  noting  that  the 
foregoing  set  of  equations  can  be  represented  by  the  more  com¬ 


pact  set: 


pH 

^ij 

if  j  =  If 

i  =  1,. 

^55 

"  ^55 

‘^56  ""  ^56 

^65 

=  ^65 

"^66  ‘  ^66 

F«. 

ID 

=  F.  . 

ID 

ifj  =  lf2 

(24) 


There  are  thus  twenty- four  equations  available  for  the  deter 

mination  of  the  twenty-four  unknown  coefficients.  Consider 

H  ““ 

first  the  solution  of  the  first  nine  equations  =  C^^  for 

H 

the  unknown  coefficients  . .  Remembering  that  the  j  are 

defined  by  equation  {16b)  with  b|j  replaced  by  d^^^  ,  the  C^^ 

then  satisfy  the  nine  equations  following  equation  (15)  with 

H 

C?j  replaced  by  C^^  and  b^^  replaced  by 
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*^11^22  ^12^12  ^1A^32  ^ 


^11^23  *^12^13  ■'■  ^^14^33  ° 


‘^21^21  ^22^11  ■'■  ‘^24^31  ^ 


^^21^22  ^22^12  ■'■  ^^24^32  ° 


^21*^23  ■*■  ^^22^13  ■'■  ^24^33  ° 


‘^41*^21  ‘^42^11  ^^44^31  ° 


*^41^22  ■'■  ^42^12  *^44^32  "  ° 


^41^23  ■*■  ^^42^13  ■'■  *^44^33  ^ 


These  are  readily  solved  in  sets  of  three  with  all  three  sets 
of  equations  having  the  same  determinant. 


^21  ^11  ^31 


^  "  ^22  ^12  ^32 


^23  ^13  ^33 


The  result  is : 


^  *^11  "  ^31^13  “  ^11^33 


^  ^12  “  ^21^33  “  ^31^3 


^  ^^14  "  ^11^23  "  ^21*^13 


^  *^21  ^12^33  “  ^32^13 


{25a 

more 


^  ^22  ^23^32  ~  ^22^33 


^  ^^24  ^22^13  ■  ^12^23 


^  *^41  ^11^32  "  ^12^31 


^  ^^42  ^22*^31  "  ^21^32 


^  ^^44  *^21^12  ■  ^11^22 


(25a) 


The  next  six  equations  l”  =  can  now  be  solved  for  the  six 
unknown  coefficients  ^2^/  ^43'  ^31'  *^32'  ‘^34*  Referring 

to  the  definitions  preceding  equation  (20)  ,  and  remembering 
H 

that  the  are  defined  by  replacing  the  b^j  by  d^^  and  the 


C^.  by  C . .  =  C. we  have; 
13  ^13  13 


‘^13*^12  ^  ^23^11  ^^43^13  ~  ^1 


dC  +dC  +dC  =I 
°13^22  “23^21  °43  23  ^2 


^^13^32  ^23^1  ■*■  ^^43^33  ^3 


^31^21  ^32^11  ■'■  *^34^31  ^4 


^31^22  ■'■  ^32^12  ■'■  *^34^32  ^5 


*^31^23  ■*■  ^32^13  ■'■  ^34^33  ^6 


These  can  also  be  solved  in  two  sets  of  three  equations ,  each 
having  the  same  determinant  L  as  the  previous  equations . 


^13  ^A2  ^2*^11  ^  ^3^14 


^^23  ^1*^22  ^2‘^21  ^3^24 


(25b) 

more 


28 


(25b) 


^43  ^1*^42  ^2*^41  ^3*^44 


^^31  ^4^21  ■*■  ^S'^ll  ^6^41 


^32  ^4*^22  ■''  ^5‘^12  ^6*^42 


d,.  =  I.d^.  +  led,.  +  led.. 
34  4  24  5  14  6  44 


These  results  are  in  the  most  compact  form  when  written  in 
terms  of  the  nine  previously  computed  d^^'s. 

The  transverse  shearing  constants  are  given  directly  as 


^^55  ^55 


*^66  "  ^66 


^^56  *^65  ^56  ^65 


Also,  from  =17'  ^8  ^  ^8' 


*^33  “  ^33  ^31^2  *^32^1  *^34^3  "  ^^31^2^  "  ^^32^1^  “  ^^34^3^ 


®3  ”  ^3  ■'■  ^^31^22^  ■'■  ^^32^11^  ^^34^12^  “  *^31^22  “  ‘^32^11  *^34^1 

(25c) 

The  equivalent  thermal  expansion  coefficients  are  then  obtained 
in  terms  of  the  previously  calculated  material  constants. 


vH 


.H 


[*'u'V9  - 

W 

■'■  ^22  **’2 

:°9 

-  '*3'*8> 

+ 

-  V6>] 

[*'ll<°4°9  - 

°6“7> 

■*■  ^22***1 

^^9 

-  ■>3'*7' 

+ 

fl2<*>3'*4 

-  °1*>6)] 

[*'ll<V8  - 

D5D7) 

^  *'22<°1 

.^8 

-  V7' 

+ 

*^12  <*>2**4 

-  D^D3)] 

‘*12 

■*11  '*14 

D«  =  - 

“*22 

'*21  **24 

■*42 

'*41  "*44 

(25d) 
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1  21^44  24  41 

°2  ^11^44  "  ‘^14^^41 

°3  ^  ^^11^24  "  ‘^14*^21 

°4  ^  ^22^44  "  *^24^42 

°5  ^  ^12^44  “  ^14^42 

^6  “  ^^12^24  "  ^^14^^22 

^^7  ^^22^^41  "  ^21^42 

^8  ^12‘^41  “  ‘^11^^42 

^9  ""  *^12^21  “  *^11^^22 


(25d) 


Since  all  of  the  equations  have  been  satisfied  without  any 
inconsistencies,  the  original  assumption  concerning  the  form 
of  the  equivalent  stress-strain  law  was  therefore  correct. 


SECTION  III 


THE  CORRESPONDENCE  PRINCIPLE 


Since  the  stress  functions 
least  class  C*^,  it  follows  that 


ID 


are  continuous  in  z  of  at 


A.  .  =  a". 

n-»-oo  ij  ij 


where  the  label  "H"  refers  to  the  homogeneous,  anisotropic 
plate.  The  temperature  and  displacement  fields  obey  similar 
limiting  laws. 


Lim  T  = 

n->oo 


Lim 
n->oo  '^i 


H 


u . 


Referring  to  equations  (18)  and  (17) ,  it  can  immediately  be 

seen  that  the  stress  components  t  ,  t  ,  t  are  identical 

xz  yz  zz 

in  the  limit  to  the  corresponding  stress  components  for  the 
homogeneous,  anisotropic  plate. 


/  ^ 

f  X  \ 

XZ 

xz 

Lim 

H 

n^oo  { 

^yz 

"yz  > 

H 

t 

T  _ 

\  ZZ  > 

^  ZZ  / 

The  remaining  stress  components  in  the  laminated  plate  are; 


Lim  ^ 
n^“  ^xx 


C*  +  C* 

^21^11  ^22^22 


+  C* 
^23^12 


I*A^  +  F*  T^ 
2  33  22 


Lim 

n-»-oo  "^yy 


C*  A^  +  C*  A^  +  C*  A^ 
^11^11  ^12^22  ^13^^12 


-r*»H 

^1^33 


+ 


F*iT 


H 


(26) 
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T  =  C*  +  C*  +  r*  —  T*A^  —  P* 

n^oo  ^xy  31  11  32  22  ^33  12  3^33  ^12-^ 

The  same  stress  components  in  the  homogeneous  plate  are: 

XX  ^21  11  22^22  ^23  12  -^2  33  22 

H  _  cH  H  H  H  ^  A®'  -  I^A^  + 

YY  ^11^11  ^12^22  ^  ^13  12  -^1^33  ^11 

tH  -  aH  +  a”  +  A^  -  iV  -  f“  (27' 

The  elastic  constants  in  equation  (27)  are,  of  course,  known  at 
this  point,  having  been  developed  in  the  preceding  section.  It 
is  now  possible  to  obtain  the  stress  components  in  the  lami¬ 
nated  plate,  given  by  equation  (26) ,  in  terms  of  the  stress 
components  in  the  homogeneous  plate,  given  by  equation  (27), 
by  eliminating  the  stress  functions  A?j  common  to  both  sets 
of  equations .  The  result  of  that  calculation  is 


n->«>  XX 


^^21*^21  ^22^11  ^23^^41^ '^xx 


^  ^^21^22  ^22^12  ^23*^42^  "^yy 


^^21^24  ‘^22^^14  ^23^44^  "^xy 


^^21*^23  ^22^13  ^23^^43  ^2^^zz 


[^22  "  ^^21^21  ^22^^11  ^23^^41^  ^22 


^^2^2  ^22*^12  ^23^42^^11 


^‘^2l‘^24  ^22^14  ^23^44^^12-'’^ 


Fi  Jt” 


more 
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Lim  ^ 

n->-00  yy 

=  <^11^21  ^ 

42^^11  + 

C*,d.T ) 

13  41  XX 

^^11^22 

t  C*2di2 

+  C*  d . „) 

13  42  yy 

+  C*  d.jT^ 

13  44  xy 

^  ^12^^13 

^!3‘'43  - 

^  [^Ii  - 

^‘^ll‘^21 

^12^11  ^13^41^^22 

^^11^22 

^12^^12 

^13^42^^11 

"^12^14 

^13^44^^12]'^” 

Lim 

n-^oo  ^xy 

''32^^11  + 

C*~d. , ) 

33  41  XX 

^^31^^22 

^32^^12 

+  C*  d.^)T® 

33  42  yy 

+  (C*j^d24 

^32^^14 

+  C*  d..)T^ 

33  44  xy 

^^31*^23 

^32^^13 

^33^43  “  ^V'^zz 

-  [^12  ^ 

(C*id2i  + 

*^32^^11  ■*■  ^33*^41^^22 

^^31^22 

^32^^12 

■*■  ^33*^42^^11 

(C*id24 

^32^^14 

^  ‘^33*^44^^12]^ 

Some  simplification  has  been  accomplished  by  making  use  of 
equation  (25b) . 

The  Correspondence  Principle  for  the  stresses,  as  sum¬ 
marized  by  equation  (28) ,  gives  explicit  relations  for  the 
stresses  in  the  laminated  plate  in  terms  of  the  stresses  in 


(28) 
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a  corresponding,  homogeneous,  anisotropic  plate  subjected  to 
the  same  boundary  conditions.  The  equivalent  elastic  constants 
and  thermal  expansion  coefficients  for  the  laminated  plate 
are  also  determined  and  are  summarized  by  equations  (25) . 

The  only  assiamption  made  in  deriving  these  relations  was 
that  the  number  of  laminations  was  indefinitely  large.  In 
that  sense,  this  theory  may  be  regarded  as  an  asymptotic 
theory.  The  error  involved  in  using  these  results  for  some 
finite  number  of  laminations  must  be  evaluated  by  the  calcu¬ 
lation  of  some  specific  cases.  One  caution  to  be  emphasized 
is  that  the  assumption  of  a  very  large  number  of  laminations 
restricts  these  results  to  balanced  laminates. 

A  computer  program  has  been  written  to  permit  the  calcu¬ 
lation  of  the  twenty  equivalent  elastic  constants  and  the 
three  equivalent  coefficients  of  thermal  expansion.  This 
program  was  written  in  Dartmouth  BASIC  computer  language  and 
is  listed  in  Appendix  I  together  with  the  results  of  several 
sample  runs  using  the  Dartmouth  Time  Sharing  Computer  System. 
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SECTION  IV 


NUMERICAL  RESULTS 

In  order  to  evaluate  the  use  of  the  Correspondence  Principle 
for  a  finite  number  of  laminations,  a  comparison  solution  is 
needed.  The  case  of  cylindrical  bending  of  an  isotropic,  lami¬ 
nated  plate  is  easily  solved  and  will  provide  the  required  com¬ 
parison  solution.  If  the  Uj^(x,y,-h)  are  functions  of  x  only 
and  the  are  functions  of  x  and  z  only,  then  equations 

(4)  and  (20)  reduce  to  the  two  compatability  equations 


As  the  number  of  laminations  increases  indefinitely,  the  strain 
distribution  A22  approaches 
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may  be 


The  integro-differential  equations  for  A^2  ^22 

solvea  approximately,  using  Vol terra  iteration  and  retaining  the 
same  number  of  terms  in  each  case. 


The  terms  containing  negligible. 

The  lower  order  terms  in  these  expansions  will  be  recognized  from 
classical  plate  theory.  The  higher  order  terms  are  corrections 
which  are  of  the  order  of  the  square  of  the  plate  thickness.  The 
boundary  displacements  are  determined  from  the  differential  equa¬ 
tions  which  result  from  applying  the  boundary  conditions  at  the 
upper  face  of  the  plate,  z  =  h.  In  a  plate  with  a  finite  number 
of  laminations,  the  longitudinal  strain  distribution  A22  IS 
seen  to  be  of  class  C°(z)  with  a  kink  at  every  lamination  inter¬ 
face.  As  the  number  of  laminations  increases,  the  kink  locations 
approach  one  another  and  the  strain  distribution  approaches  a 
function  having  a  continuous  derivative. 
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The  two  cases  may  be  compared  by  simply  comparing  the 
z-dependent  coefficients  in  A22  with  the  constant  coefficients 
in  ^22*  The  coefficients  of  the  lower  order  terms  are  identical 
so  that,  in  terms  of  the  assumptions  of  classical  plate  theory, 
the  Correspondence  Principle  gives  exact  results.  The  most  mean¬ 
ingful  comparison  is  thus  between  the  coefficients  of  the  bending 
strain  terms 

fZ  * 

C(z)  =  J  ^  (z-t)  (t+h)  dt 


1+v* 

E* 


-h  1-(V*)2 


(s+h)  ds  dt 


and 


which  are  "thick  plate"  corrections.  Results  are  shown  in 
Figures  1-3  for  the  case  in  which  the  plate  is  made  up  of  aii  odd 
number  of  laminations  of  identical  thickness.  For  the  odd  numbered 

g 

plies,  E=10  X  10  ,  V  =  .25  and  for  the  even  nximbered  plies, 

g 

E  =  1x10  ,  V  =  .45.  The  rapid  approach  of  the  solutions  for  a 
finite  niomber  of  plies  to  the  Correspondence  Principle  solution 
may  be  interpreted  as  indicating  that  a  condition  of  nearly  com¬ 
plete  mutual  constraint  exists  between  plies,  even  when  the  number 
of  plies  is  relatively  small.  Taking  into  account  the  equality 
of  the  classical  bending  terms  in  and  A22  /  the  maximum 

error  in  using  the  Correspondence  Principle  is  less  than  one  per¬ 
cent  for  three  plies  and  decreases  rapidly  for  any  greater  number 
of  plies. 
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SECTION  V 


SOME  SIMPLE f  EXACT  SOLUTIONS  FOR  EDGE 
LOADED  PLATES  AT  UNIFORI'^  TEMPERATURE 


If  the  lateral  pressure  is  zero,  then  solutions  of  equations 
(4)  and  the  boundary  conditions 


■^xz  ^  ==  Ty^(x,y,±h)  »  (x,y,±h)  =  0 


described  by 


A.  .  =  A.  .  +  zA‘l .  ,  k!  .  —  constant 

ij  ID  x-i  "  iD 


are  possible. 

The  edge  forces  and  moments  are  all  constant  and  the  trans¬ 
verse  shear  forces  are  zero.  If  the  edge  forces  and  moments  at 
the  edges  and  the  temperature  are  specified,  then  the  six  constant 

values  A^ .  are  determined  from  the  matrix  equation 
X  D 

C  A  =  N 

where 


C  = 


c?. 

21 

C^ 

^22 

cO 

^23 

rO 

^11 

^12 

'■13 

pO 

31 

C® 

^32 

„0 

^33 

C^ 

1 

C" 

21 

22 

23 

C^ 

^11 

cK. 

l2 

C^ 

^13 

1 

^31 

C^ 

^32 

C^ 

^33 

C 

^21 

C 

^22 

^23 

^11 

1 

^12 

C^ 

^13 

^1 

"31 

C^ 

^32 

<^33 

2 

C 

C2 

c2 

21 

22 

23 

C^ 

^11 

C^ 

*"12 

c2 

^13 

^2 

'"31 

c2 

^32 

c2 

^33 
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dz 


dz 


The  displacement  components  at  the  lower  lateral  face  may  then  be 
found  by  evaluating  equations  (4)  at  z=-h  and  integrating 

U3^(x,y,-h)  =  (^^2  ~  ^22)^  +  ay  +  g 

“  ^11  ■"  ^1:^^  ^^12  ”  ^1^^  -  ax  +  Y 

1  2  1  1  ^ 

U3(x,y,-h)  =  -  A22  ^  “  -^12^^  ”  ^11  ^  +  k  +  Xx  +  yy 

The  six  constants  of  integration  a,  g,  y,  k,  u  are  determined 
from  the  displacement  boundary  conditions.  If  the  displacements 
referred  to  the  mid-plane  are  desired,  they  can  be  calculated 
from  equations  (1) ,  (3)  and  (20) . 

This  set  of  solutions  describes  all  possible  cases  of  uniform 
edge  loading  and  thermal  warping.  The  transverse  shear  stresses 
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are,  of  course,  zero.  The  solutions  are  exact  only  in  regard  to 
the  matching  of  arbitrary  edge  forces  and  moments  at  the  bound¬ 
aries.  At  a  free  edge,  for  example,  the  edge  forces  and  moments 
vanish  but  the  stresses,  in  general,  do  not.  The  stress  distri¬ 
bution  near  an  edge  is,  in  general,  three-dimensional. 

In  the  case  of  unrestrained  thermal  warping,  if  the  edge  dis¬ 
placement  conditions  are  chosen  so  that  the  corner  (0,0)  of  the 
plate  is  fixed. 


u^{0,0,-h)  =  U2(0f0,-h)  =  U2(0,0,-h)  = 

9u,  Bu^ 

=  -gy  (0,0, -h)  -  (0,0, -h)  = 

.  .9U-  9u^ 

=  ■  (0,0, -h)  =  ■  "g—  (0,0, -h)  =  0 

then  the  displacement  field  is 


Ul(x,y,-h)  = 

(4 

^22^) 

/ 

1  X 

U2(x,y,-h)  = 

11 

1  y 

+ 

{42 

U3(x,y,-h)  = 

- 

^22 

2 

X 

T'  ■ 

^12 

xy 

- 

^11 

- 

- 


y 

X 


For  the  case  of  simple  tension  in  a  balanced  laminate,  the  elastic 
constants  computed  from  this  solution  agree  exactly  with  those 
found  from  equations  (25a) . 
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SECTION  VI 

APPROXIMATE  THEORY  FOR  THIN  PLATES 

The  governing  integro-dif ferential  equations  to  be  solved  are 
obtained  by  substituting  equations  (20)  into  the  compatibility 
equations  (4) .  The  unknown  quantities  to  be  determined  are  the 
four  stress  functions  f  ^22'  '^12'  '^33*  other  equation 

needed  is  equation  (17) .  By  evaluating  equations  (4)  at  z  =  -h, 
the  displacement  components  at  the  lower  lateral  face  are  obtained. 

A22(x,y,“h)  =  -^  (x,y,-h) 

Ai3^(x,y,-h)  =  -^^2 

2A^2(x/y^“h)  =  ^2  aY  ^1 

{x,y,-h)  =  -  —  U3  (x,y,-h) 

aA.,  ^2 

(x,y,-h)  - - -  U3  (x,y,-h) 

az  9y^  -3 

-yii  (x,y,-h)  .  -  ^  U3  (x,y,-h)  ,23, 

The  displacement  components  are  thus  determined  when  the  stress 
functions  have  been  specified. 

The  boundary  conditions  on  the  lateral  faces  of  the  plate  will 
be  taken  to  be: 


Txz  (x,y,±h)  =  Ty2(x,y,±h)  =  0 

Tz2  (x,y,-h)  =  0 

^22  (X/Yfh)  =  -p(x,y)  (30) 
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The  boundary  conditions  on  the  edges  of  the  plate  will  be  left 
arbitrary.  The  temperature  is  assumed  to  be  uniform.  The  order 
of  the  various  terms  in  the  compatibility  equations  can  be  ascer¬ 
tained  by  examining  equations  (20)  .  If  ^22  ^12 

assumed  to  be  of  order  h“  then  y^/  Y„„f  Yv„»  Y-r. 

" 

order  h°  and  j  and  y  are  o (h) .  The  transverse  shearing 

strains  are,  in  terms  of  the  displacements: 

Bu,  9u- 
o  1,3 

^  9  8 

Z  X 


SUr,  Bu_ 
—2  +  _! 
8  8 
z  y 


If  the  stress  functions  are  assumed  to  be  given  in  terms  of  power 
series  expansions  in  z,  then  the  displacements  are  given  by 
similar  expansions. 


u^  =  u  +  zU^ 


“2  = 


^3  ^ 

W 

If  the  transverse  shearing 

strains 

are 

^  ^xz 

= 

^1 

9W 

8x 

0 

2  Yyz 

= 

^1 

8W 

By 

0 

Therefore,  the  displacements  are 

"l  = 

u 

- 

BW 
^  Bx 

V 

3W 
^  By 

u.  = 

w 
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and  the  stress  function  expansions  in  terms  of  U,  V  and  W  are 


A„„ 

-  2 

22 

8x 

3x^ 

A,  , 

9V 

2 

^  3  W 

11 

3y 

^12 

^3U  9V\ 

,9y  axi 

3^ 

3x3y 


(32) 


According  to  equation  (31) ,  the  functions  U,  V  and  W  are  the  dis¬ 
placement  components  referred  to  the  middle  surface.  However, 
this  is  inconsistent  with  equation  (29)  which  identifies  W  as 
the  transverse  displacement  referred  to  the  lower  surface.  This 
inconsistency  should  be  regarded  as  resulting  from  the  truncation 
of  the  power  series  expansion  given  by  equation  (31) . 

If  equations  (32)  are  substituted  into  the  compatibility 

equations  (4) ,  it  is  found  that  these  equations  are  identically 

.  .  2 
satisfied  to  order  h,  that  is,  the  error  is  proportional  to  h  , 

which  is  consistent  with  equations  (31)  and  (32) .  Satisfaction 
of  the  boundary  conditions  (30)  then  results  in  three  differential 
equations  to  be  solved  for  the  three  displacement  components  U, 

V  and  W.  Equation  (17)  will  supply  the  expression  for  the  other 
stress  function  An  esjamination  of  the  integrand  of  this 

equation  discloses  that  the  terms  containing  are  small  com¬ 
pared  to  the  terms  containing  A22  and  Aj^2 

neglected.  A  suitable  expression  for  A^^  is  therefore 
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a  2  *  ;s2  *  ^2 

^ _  +  2C  _ ^ _  +  C  _ 

3x2  ^  ^"'32  3x3y  ^  ^12  .^2 


*  ^2  *  a2  ^  ^2 

^23  ^^33  axSy  ^  ^13  3y2 


The  governing  differential  equations  are 


Lj^U  +  L2V  -  L3W  =  PjL 


L.U  +  L^V  -  L.W  =  Po 
4  5  6  2 


^12  / 


L^U  t  Lg 


LgV  -  LgW  =  P3 


L,  =  C°~  — 
1  22 

T  _  1 


^2  ■  7  ^23 


/  o  1  o  \  3^ 

^  (^32  7  ^23) '^^5^ 

O 


1  r"  ^  „ 
x3y  ^  ^33  Sy2 


^  (^21  2  ^33)  ^ 


^2  ,  .2 
x^  +  ^31  ^ 


\  Ai-i.  ^ ^ j  ju  y 

^3  =  ^22  ^  i^32  "^2 -  (^21  ^33)“^ - 2 

3  22  9x^  '  dx^^Y  '  8x9y^ 

-1^ 

31  3y3 

L.  _  r*  9^  .  /p°  .1  N  +  i  r“  ^ 

^  ^32  3x2  ^  Vl2  ■?  ^33/  3x3y  ^  2  ^13  3y2 

9  2  2 

1  o  a  /  o  1  °  \  3  03 

^5^2  ^33  ■*■  (Sl  ■'■  2  ^13)  33c3y  ^11  ^ 

^22  (‘2t'2  ■'■  ^\-^)  “V“  +  (p\l  +  ^13)  - F 

^  \  23/  ax^^y  '  21  13/  9x9y^ 

.cl 

^11  3y3 

^7  ^  (^*^22  “  ^22)  ■'■  (^^^^32  "  ^^32  I  ^*^23  "  '2^23)  9x2ay 

■'■  (^^12  ”  ^12  ■*■  ^^33  “  ^33)  9x3y2  '*'  2  (^^13  "  ^13) 


(c]^  +  C^,)  -\-—  +  (Cg,  +  cj'g)  j 

\  12  33/  Sx^ay  '  31  13/  3xay2 


(33) 

continued 


1  (hcls  -  =«)  Ip 
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Equations  (34)  are  a  set  of  three  differential  equations  of  eighth 
order  in  the  three  displacement  components.  Before  considering 
the  general  case,  it  will  be  helpful  to  consider  two  special  cases 
which  are  much  simpler  and  which  will  provide  some  insight  into 
the  nature  of  the  solutions. 
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1.  Homogeneous,  Isotropic  Material 

This  is  a  single  ply  plate  of  isotropic  material.  The  elas¬ 
tic  constants  of  the  material  are 


^11 

^22  ^33  ^ 

1 

E 

^12 

^21  "  ^13  "" 

^31  ^  ^23  ^32  “  "  E 

^44 

^  ^55  ^66 
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E 
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""  ^24  ""  ^34  " 

*'41  “  *^4  2  “*’43  “  ^’Se  '  *^65  “  ° 

^11 

=  c  -  ^ 

^22  i_^2 

^12 

^21 

E 

^33 

1+v 

^13 

^31  ^  ^23 

C32  =  0 

fh 

c! . 
ID 

=  C . .  dz  = 

=  2hC,. 

C.^ 

ID 

=  0 

c?. 

ID 

2 

=  J  c,. 

-h 

dz  =  4  C. . 

1  J 

^1  = 

■  ^2  =  ^3  =  “ 
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au  av 
ay  ax 


Eh^ 

'a\ 

3^W 

3(l-v2) 

a^2 

ay2 

Eh^ 

■  2 
a  w 

2 

aw 

^  \  . - 

3a-v2) 

3y2 

ax2 

Eh^  a^w 

(i+v)  axay 


Eh^  a 
3(1-v2) 

■  A-' - -  V^W 

3  (l-v2)  9y 


and  the  stresses  are 


E  Z  ^  M 

l-v2  [3x2  J 

EZ  r  a^  .  1  TT 

r  iF  J 


EZ  aw 

i_^2  axay 

f  z^  zh^ 

V  6  T"  ”  3 


E(zW)  _1  ,2„ 

2(1-v2) 


,  EC^  -,.l^.  j_  2^ 

2(i-v2)  5y 


If  the  plate  is  simply  supported, 

TT  rT  •  mirx  .  niTv 

W  =  sin  — —  sin  , 

mn  a  d 

and  the  transverse  shearing  stresses  are 
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to  (jJ 
tr  ^ 
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2.  Balanced  Laminate 

A  balanced  laminate  is  defined  as  a  laminated  plate  for  which 
=  0.  The  material  property  distributions  are  even- 
symmetric  functions  with  respect  to  the  mid-plane  of  the  plate. 

In  this  case,  the  operators  and  Lg  are  zero  and  the  in¬ 

plane  displacements  can  thus  be  determined  independently  of  the 
transverse  displacement.  By  inspection,  the  particular  integrals 
of  U  and  V  are  of  order  h*^.  Substituting  these  into  the 
last  of  equations  (34) ,  the  equation 

LgW  =  -P  (l+0(h^)) 

is  obtained.  Since  the  terms  in  this  equation  arising  from  U 
and  V  are  of  order  h  compared  to  unity,  they  may  be  discarded 
with  the  result  that  the  transverse  displacement  is  approximately 
uncoupled  from  the  in-plane  displacements.  Since  U  and  V  are 
of  order  h^W,  they  may  be  ignored  compared  to  W.  However,  this 
does  not  necessarily  mean  that  the  effect  of  U  and  V  on  the 
edge  forces  and  moments  is  negligible;  this  remains  to  be  evalu¬ 
ated.  Contributions  to  various  terms  will  be  considered  negli- 
gible  only  if  they  are  of  order  h  compared  to  unity.  The  edge 
forces  and  moments  are 


N..  =  c:,  iX  +  c:.  ig  +  4 


X  '21  3y  22  3x  I  23  \^ay  Sxy 
h 


-h 


^2  "^zz  ^22^ 
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09V  o  9U 

N  =  r  -f  c  + 

y  ^11  ay  ^  ^12 


1  =  fm  w\  [•'  * 

I  (^3y  *  axj  J_j^^i 


^zz  ■"  '’ll'' 


N 


xy 


o  8V  o  ^1  o  y^u  W 

^31  3y  32  9x  2  ^33V3y  3x 


-h 


^3  "^zz  ~  ^12^ 


-  M 


X 


=  C 


2  .3^W  L  ^2  9^W  ^  ^2  9^ 

21  975-  ■*■  ^32  9x2  23  9x9y 


-  „  =  c?,  +  c?,  +  c? 


•11  3y2 


'12 


9x' 


'13  9x9y 


_  M  =  +  c2  ^ 

^xy  ^31  av2  32  9x2  33  9x9y 


“  Q„  =  ^ 


X 


9 

9y 


_2  9^W  ,  _2  9^W  2  9^W 

^31  9y2  ^32  9x2  ^33  9x9y 


“  Qy  =  ^ 


^2  9^W  ^  ^2  9^W  p2  9^W 

9x  r21  9y2  ^22  9x2  ^23  9x9y 


^2  9^W  .  ^2  9^W  .  „2  9^W 
^31  9y2  *-32  9x2  ^33  9x9y 


2  9^  2  9^  2  9^ 

9y  1  ^11  772  ^  ^12  9x2  +  ^13  9x9y 


Terms  have  been  discarded  from  the  moment  and  transverse  shear 

2  . 

force  expression  which  are  of  order  h  .  Since  the  in-plane  dis¬ 
placements  do  not  appear  in  the  moments  and  transverse  shears  and 
the  transverse  displacement  does  not  appear  in  the  edge  forces# 
the  problem  is  truly  uncoupled  at  this  level  of  approximation. 

The  transverse  pressure  terms  which  appear  in  the  expressions  for 
the  edge  forces  are  due  to  the  Poisson  constraint  effect.  The 
presence  of  these  terms  means  that  the  vanishing  of  the  in-plane 
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displacement  terms  does  not  ensure  that  the  edge  forces  vanish. 
This  must  be  taken  into  account  in  the  specification  of  the  bound¬ 
ary  conditions. 

In  the  determination  of  the  stresses,  the  dominance  of  the 
transverse  displacement  is  clearly  evident  since  the  in-plane 
displacements  contribute  only  to  O(h^)  compared  to  the  trans¬ 
verse  displacement  terms. 


,  _  fr*  sfw  *  9^ 

^xx  “  I  ^21  9y2  ^  ^22  3x2  ^  *-23  3xay 


■  -[■ 


■  *  a^w  ,  p*  a^w  .  p*  a\ 

^11  ^^2  ^12  3x2  ^13  axay 

*  a\  .  p*  3^w  .  p*  a^w 

*^31  3y2  ^32  3x2  33  axay 


Z  +  ^22  T 


Z  +  T 


Z  -  Fi2  T 


t(z-t)  C22  ^23) 


*  )  3 

ax^ay 


/  *  *  /* 

(^12  ^21  ^^33)  3v23v2  r^31  ^ 


*  \  a' 

13fe 


31  ^13/3xay 


,*  a^ 

'11  av‘^ 


W  dt 


ic 

3 

aw 

+  Ic* 

4- 

C*  \ 

3 

aw 

22 

ax^ 

^  r23 

^32) 

ax^ay 

4* 

c*  ) 

a\ 

•4 

c*  ^ 

^  1  t 

^33) 

3x3y2 

T 

^313 

y3  J 

a^w 

+  /c* 

i 

c*  \ 

a’^w 

32 

3x3 

^  \pl2 

^33) 

ax^ay 

t  dt 


+  (C*  +C*  )  -iiw  +c*-^ 

13  31  3xay2  11  ay3 


t  dt 


The  transverse  shearing  stresses  automatically  vanish  at  z— h 
due  to  the  even-symmetric  material  property  distributions. 
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3 .  The  General  Case 

If  the  laminate  is  not  balanced,  none  of  the  operators  in 
equations  (34)  vanish  and  the  in-plane  displacements  and  the 
transverse  displacement  are  coupled.  The  edge  forces  and  moments 
will  be  needed  in  order  to  formulate  boundary  conditions. 


9V  o  dU  1  o  h\J  3V\ 

^21  Ty  ^22  3x  2  23  3x J 

f**  i!  ,  dz  +  T  -  ci,  -  cL 


,1  3^W 


22  33^2  23  3x3y 


o  3V  °  9U  +  1  c°  (11  + 

^11  3y  ^12  3x  2  ^13  \3y  3xy 

-r*  J  .  T,°  m  9^  ^1  3^W 

,  ^1  ■'zZ  T  -  _2  "  ^12  9x2 


a  3^w 
'13  3x3y 


o  8V  o  ^  1 

31  ^  ^32  3x  2  33  I  3y  ^ 


I* 
-h  3 


,  dz  -  F°  T  -  2%  -  cl  ^ 

zz  12  31  3y2  32 


,1  3^W 

'33  3x3y 


pl  IX  +  XU  +  i 

^21  3y  ^22  ax  ^  23 


T  dz  +  T 

zz  22 


.2  3^  p2  3^W 

'22  9x2  “  23  3x3y 


=  ^  11  +  X  ^11  +  XX'\ 

Hi  3y  ^12  3x  2  ^13  \^3y  3Xy 

fh  .  1  2  9 

z  I*  T,,  dz  +  fIt  T  -  c:f,  — ^  - 

,  1  zz  11  11  9y2 


'13  3x3y 


pl  IX  +  pi  11  +  1  pi 

^31  3y  ^32  3x  ?  33 


11  IX 

3y  3x, 


^  z  I*  T  dz  -  F^.  T  -  C^,  1-1  -  C32  1^1 
-h  3  zz  12  31  9y2  32  g^2 


,2  3^W 

'33  3x3y 
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I  (^23  “  ^^23)  ^  ?  H3  ■  ^^33)  3y2 

■’■  J  “  ^^23  ■'■  ^^32  “  ^^^32^  dxdy  ■''  (^22  "  ^^22) 


dt  dz 


,1  \  8  W 
'22/  ■:>^3 


''^  p  r  I*  _L  +  I*  _il  T  dt  dz 

J-h  J-h  I  3  3y  ^  "2  8xJ  ^zz 

-  (c^^  -  hcl^  -  (C33  -  hC33  t  cl^  -  hC2i) 

-  (cL  -  -^2  ^  C^3  -  -L)  iS?  -  H2  -  -^22)  ^ 

=  ^€32  -  11032)  ^  +  I  (C33  -  hC33  +  2cj;2  “  2110^2)  aSTSy 

+  1  {^13  "  ^^13)  I  (^33  "  ^^33)  I  (^13  "  ^^13 

+  2033^  -  2hC3 J  — 

-  I?  -  (=33  -  >>=^3  -  -  ^<^12)  A 

-  (cL  -  -L  ^  cL  -  -L)  -  {<1  -  -ii)  0 


+  2C 


The  stresses  are 


T  =  c*  +  C*  ^  + 
XX  ^21  9y  ^22  8x 


123  /^  . 

2  \9y  9x/ 


V21  9y2 


9^W  . 

— “r  C  o  • 


+  c*. 


23  9x9y/ 


2  -^zz  ^22  ^ 
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C*  ^  +  C*  ^ 

^11  9y  ^  ^12  9x  ^  ^  \^9y  ^  9xj 


fn*  3^W  ^  9^  ,  a^W  \ 

(  ^11  ^  ^12  ^  ^13  93^y  ^  "  ^1 


1  ^22 


+  F*  T 


ic 

n*  ,  ^*  9U  ^33  /  9U  9V\ 

^31  9y  +  ^32  ^  +  “5“ 

-  fc*  +  C*  ^  +  c*  z 

\^31  3y2  ^  ^32  9x2  +  *-33  9x9yy  ^ 


-  I->  T 

3  zz 


-  F?o  T 


^22  3x3  (^^32  ^23) 


+  (c*  +  c*  )  — ^—-  +  i  elj'o  u 

V  12  33/  ^  9y^ y 

+  (j  c*  •  +  (c*  +  c*  A  — I —  +  (2c*  +  ~  c*J  — 

\2  23  3x3  \  21  33;  3x23y  \  31  2  13/3x3 

^  “1  . 


r  * 

t(2-t)  C*  ^ 
Lh  22  3x‘» 


■*■  (2^32  ^23)  ■*■  (^21  ^12  2^33) 

+  (20^,  +  ct^  — ^  +  C^,  W  dt 

\  31  13/  3x9y3  11  3y'+ 

^  P*  9  ^U  ,  f„*  ,  1  „*  \  9  ^U  ,  1  9^U 

j-h  ^  ^  ^  ^ 

,  2..  .  -  -  2..  -  2.. 


.  1  „*  9'^V  /p*  1  ^*  \  9  V  ^  9^V 

7  ^23  3x2  r21  2  ^33/  9x9y  ‘^31  iy2 

^  +  (^23  ^32)  ilaTJ  +  (^21  +  ^33) 

^  ^  — I 


,*  9  W  \  .  ,4.  ^ 

'  ,  — ^  t  dt  + 


V*  9  .  T*  9  1  JO. 

^3  9y  ^2  9x  "^zz 
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1  _*  3^U 

7  ^13  ^ 


+ 


1 


7  ^13) 


9^V  . 

5  x'8  y  ^11  9y2 


8\ 

8x3 


+ 


+ 


r*  \ 

12/  ax^Sy 


+ 


+  c 


31 


8\ 

8x8y2 


1 


dt  + 


z 

^  -h 


AT* 

8x  ^  7y 


zz 


dt 


Let  us  consider  first  the  case  in  which  the  laminate  is 
severely  unbalanced,  that  is  ct  .  and  h  are  of  the  same 

order  of  magnitude.  Then,  W  =  0(h“^)  and,  from  the  first  two 
of  equations  (34),  U,  V=0(h"^).  Thus,  and  P2  are  of  order 

h  and  are  therefore  negligibly  small  compared  to  1^3^  and  LgW 
which  are  of  order  h  In  the  edge  force  and  moment  expres¬ 

sions,  the  Poisson  constraint  terms  which  contain  the  quantities 
I*  are  negligibly  small.  None  of  the  other  terms  can  be  dropped. 

Now,  let  us  suppose  that  the  laminate  is  only  slightly  un¬ 
balanced,  so  that  C^j  =  O(h^) .  Then,  W=0 (h  as  before  and 

U,  V  are  of  order  h° .  None  of  the  terms  in  the  first  two  of 
equations  (34)  can  be  dropped  but  the  terms  L.yU  and  LgV  in  the 
last  of  equations  (34)  are  negligibly  small.  No  terms  can  be 
dropped  from  the  edge  force  expressions  but  the  U,  V  and 
terms  can  be  neglected  in  the  expressions  for  the  moments  and  the 
transverse  shear  forces. 

If  the  equations  for  the  general  case  are  to  be  equally  appli¬ 
cable  regardless  of  the  degree  of  unbalance  of  the  laminate >  then 
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no  terms  can  be  dropped  from  the  differential  equations  (34) . 

The  only  terms  which  can  be  dropped  from  the  edge  force  and  mo¬ 
ments  are  the  terms  in  the  bending  moment  and  transverse 

it 

shear  forces.  The  terms  containing  in  the  stresses  can  be 

dropped  since  they  are  negligible  in  both  of  the  above  cases. 

The  equations  derived  in  this  section,  which  are  the  equa¬ 
tions  of  ordinary,  laminated  plate  theory,  have  been  derived  from 
the  equations  of  elasticity  for  laminated,  anisotropic  plates  by 
making  two  assumptions: 

(1)  The  displacements  are  linear  functions  of  the  z  coordinate. 

2  .  . 

(2)  Terms  of  order  h  are  negligible  compared  to  unity. 

The  assumptions  customarily  made  in  deriving  the  ordinary  lamin¬ 
ated  plate  equations  are  therefore  dependent  only  on  these  two 
assumptions.  The  expressions  for  the  transverse  normal  and 
shearing  stresses  which  have  been  derived  are  those  which  would 
be  obtained  from  the  ordinary  theory  by  using  the  equilibrium 


If  the  "first  order  theory"  of  laminated,  anisotropic  plates  (the 
present  simplified  theory)  proves  inadequate,  then  the  reason 
must  be  due  to  the  above  two  assumptions  and  their  modification 
must  then  lead  to  improvements  in  the  accuracy  of  the  theory. 

An  obvious  modification  is  to  assume  that  the  displacements  are 
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quadratic  functions  of  the  z  coordinate.  It  should  be  remem¬ 
bered,  however,  that  any  representation  of  the  displacements  in 
terms  of  continuously  differentiable  functions  is  inherently 
limited  due  to  the  fact  that  the  displacements  are,  in  reality, 
elements  of  class  C**  (z)  . 

The  expansion  (32)  is  a  formal  one  and  the  relative  magni¬ 
tude  of  each  of  the  terms  is  unknown  at  the  start.  In  facv,  it 
has  been  shown  that,  in  the  general  case,  the  first  and  second 
terms  in  the  expansions  for  U2^  and  U2  are  of  the  same  order  of 
magnitude.  For  a  balanced  laminate,  the  second  terms  are  actu¬ 
ally  larger  than  the  first  terms.  One  can  only  hope  that,  in 
creating  a  higher  order  theory,  the  correction  terms  eventually 
turn  out  to  be  smaller  than  the  first  few  terms  in  the  expansion. 
This,  of  course,  must  be  verified  after  the  fact. 
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SECTION  VII 


A  PLATE  THEORY  CORRECTED  FOR 

transverse  normal  and  shearing  deformation 

As  was  pointed  out  in  the  previous  section,  the  ordinary 
laminated  plate  equations  are  based  on  the  assumption  of  a  dis¬ 
placement  field  which  is  linear  in  the  z  coordinate  and  on  the 
neglect  of  all  terms  of  order  h  .  Even  with  this  limitation, 
however,  the  complete  stress  field,  including  the  transverse 
normal  and  shearing  stresses,  can  still  be  obtained,  at  least  to 
an  accuracy  consistent  with  the  stated  assumptions.  For  moder¬ 
ately  thick  or  highly  anisotropic  plates,  the  previous  assumptions 
are  open  to  question.  In  order  to  develop  a  theory  which  is  cor¬ 
rected  for  transverse  normal  and  shearing  deformation,  both  of 
the  previous  assumptions  will  have  to  be  discarded,  since  as  will 
be  shown,  they  are  actually  not  independent  of  one  another. 

The  linear  displacement  field  given  by  equation  (32)  implies 
that  the  normal  derivatives  of  the  displacement  components  are 
continuous  across  ply  interfaces,  an  implication  which  is  not  sup¬ 
ported  by  equations  (1)  and  (2) .  An  examination  of  these  equa¬ 
tions  discloses,  however,  that  the  discontinuities  are  propor¬ 
tional  to  the  transverse  shearing  strains  and  therefore  must  be 
small  if  the  shearing  strains  are  small.  It  seems  reasonable, 
therefore,  to  develop  a  higher  order  plate  theory  by  assximing 
arbitrary  but  small  additions  to  the  displacement  expressions 
given  by  equation  (32)  which  are  of  class  C°(z).  The  "smallness" 
of  the  corrections  will  be  of  value  in  simplifying  the  resulting 
analysis.  Since  the  corrections  are  proportional  to  the  trans- 
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verse  shearing  strains  and  since  the  transverse  shearing  stresses 
inust  vanish  on  the  lower  lateral  surface ^  the  correction  terms 
must  also  vanish  at  z  =  -h. 

In  accordance  with  the  previous  remarks,  the  displacements 
will  be  taken  in  the  form 

“1  =  u  -  ^ 

u,  =  V  -  Z  +  f''  V*  at  C35) 

2  8y 

j  -h 

* 

u-5  =  W  +  W  dt 

J  “h 

These  modified  expressions  are  thus  consistent  with  equations 
(1)  and  (2)  and  with  the  boundary  conditions.  It  should  be  noted, 
however,  that  the  interpretation  of  the  displacement  components 

U,  V  and  W  as  the  displacements  referred  to  the  middle  sur¬ 

face  of  the  plate  is  no  longer  possible. 


^1 

(x,y,o) 

=  u 

+  1 

J 

rO 

u* 

dt 

^2 

(x,y,o) 

=  V 

+■  ' 

J 

rO 

-h 

V* 

dt 

^^3 

(x,y,o) 

=  w 

1 

J 

'O 

-h 

w* 

dt 

If  the  corrected  theory  is  to  be  easily  reducible  to  the  previous 
theory,  this  inconvenience  will  have  to  be  tolerated. 

This  approach  may  be  compared  to  an  alternative  method  which 
consists  of  a  Volterra  expansion  of  equations  (4)  and  (17)  for 

the  stress  functions  ^22*  ^12  ^33*  Equations  (29) 

and  (32)  are,  in  fact,  the  first  few  terms  of  such  an  expansion. 

The  displacements  u^(x,  y,  -h) ,  U2{x,  y,  -h)  and  U2(x,  y,  -h) 
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are  carried  along  as  unknowns  and  later  determined  by  applying 
the  boundary  conditions  at  the  upper  lateral  face, 

T3,2(x,y,h)  =  Ty^(x,y,h)  =  0 

=  “P 


These  result  in  three  differential  equations  for  the  determina¬ 
tion  of  the  displacements  Uj,(x,y,-h),  U2(x,y,-h)/  U2{x,y,-h)» 

This  procedure  is  thus  directly  comparable  to  that  recommend 
at  the  beginning  of  this  section. 

The  expressions  for  the  strains  which  follow  from  equations 
(35)  are 


^yy  ' 


3U 

z  + 

fZ 

au* 

3X 

^  ^  1 

-h 

ax 

av 

-  ' 

•z 

* 

av 

ay  “ 

3y2  J 

-h 

ay 

2  Y  =  2  A 

^xy 


_  8U  3V  o  ,  . 

12  "  ciy  ^  ^ 


^  r  su* .  9v* 

_h  L  3y-  ^  ^ 


=  W 


ju  3W* 

2  =  U*  +  I -  dt 

XZ  3x 


*  dt 

2  =  V  +  3y 

y2  J-h 

Using  the  fourth  of  equations  (20)  gives  an  expression  for  W  . 

r  Ci  T  T  r\  2- »  ^  Z  c. ,  T  ^  *1 


w*  =  I*  —  -  z 
4  3y 


.2^-  rZ 

P-  at 
5y2  i-h  ’y 
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+ 1*  -  z  ^  —  dt 

^5  8X  9x2  ^  J  .  9X 


^  T*  1  /3U  .  9V 

+  •«-  1  Trrr  "rTTr 


6  ?  v^y 


1  /su*  .  9V* 

7  J.h 


9x9y 


Terms  of  order  have  been  neglected. 

The  last  two  of  the  above  equations  can  now  be  used  to  obtain 
expressions  for  U  and  V  . 


U*  =  2  Y, 


*  9W* 

^  -  ].h  ^ 


Assuming  that  U,  V  and  W  are  given,  the  last  three  equations, 
when  written  out  in  full  (using  equations  20)) ,  together  with 
equation  (17) ,  are  a  set  of  four  integro-dif ferential  equations 
for  the  four  unknown  quantities  U*,  V  ,  W  and  A^^  iu  terms 


of  U,  V  and  W. 


u*  =  2  b*5 


-  (c*  ^ 

\^21  9x 


+  r*  U  -  t 

^  ^31  9yj  Uy  3y2 


“  (^22  ^  ^32  ^ 


(c*  L-  +  c*-  — 

r23  9x  ^  ^33  9y 


1  n\3  ,  9V 

I  I  ^  ^ 


>  t  i 

9x9y  2 


-h  V  3y 


+  I 


2  9x 
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^3  /  ^33 


+  2 


.*  3  .  p*  3  \  /S 
'31  ^  ^  ^11  ‘Sy)  \} 


3V* 

J-h  ^ 


\  /^*  3  .  ^*  3  \  /  3U 

7  “  (^32  lx  +  ^12 


3U*  ,  \ 

J  -h  ^  7 


*  3  ^  3  \  //I  .  av' 

C33  ^  ^13  Ty)  (l7\.^  ^ 


-  t  ^  +  i  r  ^  ds)  +  fl* 

^  axay  +  ?  j  .  Uy  ax  i  ^7 


*3  .  _*  a  ' 

3  ^  ^ 


)  ""A 


1.  r  I*  -  t  ^  ^  as 

ax  4  \^ay  ^  3y2  J^h  ^ 


.  T-*  /  au  .  a^w  ,  au* 


.  1  -»  /su  .  3V  ,.  3^W  .  /3U*  ^ 

I  ^6  (  ay  ai?  "  ax^  i-hV^^ 


)l(sf  •  s^)  '•)  •  ■;  *»]  *« 


V*  =  2  b*5 


-  fc*  ~  +  c*  - 


.  3V*  ,  ' 

/ 

.  au* 


'31  ay  "'21  axM  ay 


.*  i_  +  c*  M(  ^ 

'22  ax  32  ayj  Vax 


av  _  .  3^w 

TTT  *“  ^  *  ’’ '  n ' 


\  'u*  ,  \  _  7,*  1-  +  C*  |-\fiA£+  |X' 
J_j^  Jx~  ^  )  \'^23  ax  ^  *^33  ^  H  2\^ay  ^ 


more 


-  t  22L  +  i 

3x3y  2 


■3  ay 


au*  av* 
ay  ax 


ds  +  I 


2  ax 
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+  2  b 


-  fr*  L^^r*  f  ^  -  t 
rai  Bx'  ^11  3yM  8y  3y2 


t  ^ 
-h  3y 


')  "(^32  lx  ^12  ly)  (H 


.  9U*  j  \  /r,*  9  n*  9  \  /^  1  /  9U  .  9V 

^  "  (^^33  W  +  ^13  ^  j  (  7  (  3y  7x 


.  3^.W  ^1  /3U*  .  3V* 

"  1  [W"  “ 


ds)  +(I 


■^1  ^  )  ^33 


-h  ^  4  3y 


’3  3x 


*  /  3V  *  3^W  ,  f*^  3V*  ,^' 

4  {^-^W  i-h“^ 


r*  -  t  +  f*"  C  as^  a  I 


5  V9x 


1  T*  /3U  ,  ^ 
I  6  I 3y  3x 


.  2t  ^  + 

3x3y 


3U  3V 
3y  3x 


*  \ 

-)ds) 


+ 17  A33  dt 


=  I 


1  ,* 


3V 

ry 

3^ 

9y  " 

io 

Sy2 

•y 

3^ 

3x 

3x^ 

/  3U 

IX  . 

1  9y 

■7* 

3X 

^  3V* 


3U* 


2z  ^  + 
3x3y 


■'■  ^7  ^33 
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A 


33 


(z-t) 


-h 


.*  ■  a- 


+  c 


*  ■  3“ 


v2 


'21  3x2  "  m  3y 


+  2C 


31  8xSy^ 


( 


3V  .  9^W  ^ 

■?rrr  t:  - rr  + 


3y 


+  2C 


ay 

3^  ^  /3U 


"  as'l  +  (c*  ’ 


3y 


+  C 


*  ■  ■  3^ 


22  3x^  ay2 


32  8x3yy 


- 1  ^  + 
9X  3x2 


3H* 


-h 


oX 


ds 


23 


3x' 


+  C 


13 


^  +  2C* 


2  ^ 

\  / 


3y^ 


3 3  3x3 y J  \  3 y 


3U 


3V 


(39) 


2t 


-  I 


2 

^  JL 

Tx3y 

.*  32 


ds 


-h  v"^ : 


3x 


-  iT?  ^33 


dt 


Equations  (36)  through  (39)  are  exact  equations  of  elasticity. 

The  equilibrium  equations,  the  internal  continuity  conditions,  and 
the  boundary  conditions  at  z=-h  have  all  been  identically 
satisfied. 

Equations  (36)  through  (39)  are  integro-dif f erential  equa¬ 
tions  of  Volterra  type.  Solutions  of  these  equations  obtained  by 
Volterra  iteration  will  be  rapidly  convergent  since  the  range  of 
integration  (2h)  is  small.  The  Volterra  iteration  solutions  coin¬ 
cide  with  the  asymptotic  solutions  of  these  equations  in  powers 
of  the  thickness  parameter,  h.  The  powerful  existence  and  unique¬ 
ness  theorems  of  Volterra  equation  theory  guarantee  that  these 

solutions  are  regular  for  all  values  of  z.  Further,  since  U*  and 

“T  **“2 

V*  are  of  order  p  h  and  W*  is  of  order  p  h  ,  the  correc- 

2 

tions  to  equations  (35)  are  of  order  h  compared  to  the  terms 
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U,  V  and  W.  The  one-teonn  Volterra  iteration  solutions  of  equa¬ 
tions  (36)  through  (39)  will  therefore  produce  a  plate  theory  which 

2 

is  correct  to  within  terms  of  the  order  of  h  . 


U*  = 


-  2b 

-  2b 


55 
* 

56 


-h 
z 

-h 

/,* 


L*  U  +  L2  V  -  L3  W 


lJ  U  +  L*  V  -  L*  W 


dt 

dt 


1  8' 


i  aw)  “ 


(40a) 


1  8 


+  I 


*  3 ' 


^  6  9x2  4  3x3y 


V 


-  t 


3x' 


+  I*  9  +  I* 

6  3x23y  4 


3x3y' 


W 


dt 


V*  =  -2b* 

V  <^Dg5 


l1:u  +  V  -  l,  w 


?  ^  -h  L 

z 


-  2b 
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-h 


L*U  +  L*  V  -  L*  W 


dt 


dt 


z 

-h 


(40b) 


T*  1  T*  n  +  -fl*  -  +  1*  —  )  V 

^5  3x3y  ?  6  3y2/  2  \  6  3x3y  4  3y2 


-  t 


+  I 


*  3' 


5  3x23y  ^  3x3y2 


+  I 


3y^ 


W 


dt 


W*  = 


T*  S  +  i  I*  1- 

•^5  3x  2  6  3y 


U 


1  ^  .  j* 

2  6  3X  4  3y 


^2  „2 
«  I  T*  ^  .  T*  ^ 

“  '  ^5  ^6  3x3y 


V 


+  I 


(40c) 


2  1 


4  9y2 


W 
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j  _  (z  -  t)  L*  U  +  L*  V  -  L*  W  dt 


(40d) 


The  operators  are  defined  by 


Li  = 


L*  dz 


They  are  thus  obtained  from  the  operators  following  equations  (34) 
by  replacing  the  material  constants  ct .  with  the  material  prop- 

^  J  ^ 

erty  functions  ,  C^.  with  z  C^.  ,  and  C?.  with  z^  C..  . 

It  is  now  possible  to  write  corrections  (SA^j)  for  transverse 
normal  and  shearing  deformation  to  the  stress  functions  A^j , 


5  ^22 


Z  3U_ 
9x 

-h 


”oU*  3V*  ]  , 


* 

(z-t)  L-  U  ds  +  Lp  V  ds 

-h  [  '  -'-h  ^  ''-h 

(^2  ^  ^1  ^^3  9^^)  ^33 


These  corrections  will#  in  turn#  lead  to  correction  terms  in  the 
differential  equations  (34)  .  Let  the  correction  terras  which  are 
to  be  added  to  the  left-hand  sides  of  equations  (34)  be  designated 


by  '^1  »  ^2  ^3  ^respectively . 
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■  II  [ 


^  *  *  1  * 

U  dt  +  L,  V  dt 

-h  ^-h 


= 


_ 

rZ 

fZ  * 

L* 

U*  dt  +  Le 

V  dt 

*-h 

4  J 

-h  ^  ^ 

-h  J 

dz 


dz 


h 

it 

Z 

f  it  ^  *  f 

(h-  2) 

^  -h 

J 

U  dt  +  Lq 
-h  ^ 

(41) 


V*  dt 


-h 


"  ^^2  ^  +  2I3 


dz 


-h 


(z-t) 


u  +  L* . 

^  -  l;  w 

dt 

i— 

The  full  expressions  for  these  corrections  are,  of  course,  obtained 
by  inserting  the  expressions  for  U*  and  V  from  equations  (40) . 
An  examination  of  the  results  of  the  expansion  shows  that  the 
coefficients  of  the  correction  operators  are  of  the  form 


ijkJljnn 


i, j=l,2,3 
TfPg 

^ijk  - 
i , j=l ,2,3 


h  cz 

P  r*  .  I 


z^  C  .  . 


-h 


ij  J-h  i- 


C*  ds  dt  dz 


mn 


k, £  =  5,6 
rt 


m,n=l,2,3  p,q  =  0,1 


-h 


zP  Ct. 

13 


-h 
k=4,5,6 


ds  dt  dz 
p,g=0,l 


(42) 


13k 

p=0,l 


-h 


zP  I* 
1 


-h 


(z~t)  t'^  dt  dz 


cj— 0  ^1/2  1—1 1^  2  ^  3  3  k— 1  /  2  /  3 


There  are  1296  possible  coefficients  of  the  type 


108  coefficients  of  the  type  K 


pq 

i jk  * 
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The  expansion  of  equations  (42)  can  be  facilitated  by  using 


the  notation 


ijkJimn 


''k,  C 


for  coefficients  of  the  first  type,  and 


KPg  =  kP.  kJ 

ijk  13  k 

for  coefficients  of  the  second  type,  with  the  understanding  that 
the  indicated  operations  are  not  commutative.  The  corrections  to 
the  differential  equations  (34)  may  then  be  written,  in  operator 
notation,  in  terms  of  a  set  of  coefficients  ,  m.  ,  n.  which 

JL  JL 

must  be  calculated  by  a  computer  program. 

p,4  j,4  g4  g4  2  4 

S  L,  =  Jl,  -V  +  5 - 1 - T  +  - T  +  - r 

1  1  9x^  ^  9x^3y  ^  3x^9y^  ^  3x3y3  ^  3y^ 

.4  .4  4  „4  .4 

S  =  if.  -V  +  ^-7  — z —  +  — Z — T  +  - T  +  “IT 

^  ^  3x‘*  ^  3x^3y  ®  3x^3y^  ^  3x3y^  3y‘* 

5  5  5  5 

fi  r  =  a  -J_  +  £  3  +  0  3  +  0  3 

^  3x^  3x^3y  3x^3y^  3x^3y^ 


3x3y^ 


'16 


3  4  9  ^  9  ^  3 

6  L.  =  m,  — -  +  m«  5 —  +  m-  — +  m.  — 

4  1  9x^  2  9x^3y  ^  3x^3y2  4  9^ 


,4  (43) 


- r  +  m^  — r 

9x9y3  5  9y4 


g4  34  9^  9^ 

6  L-  =  ITl^  - P  +  lU-  - -h  iUq  '  ■  ■  ■■  - +  IIIq  ^  ‘T  4*  IUt 

5  6  3x4  /  3x39y  8  3x29y2  9  3x3y3  10  3y^ 

95  g5  g5  gS 

^  ^6  =  "^11  +  ""12  "‘13  ’^14  3x29y3 

I  9^  9^ 

15  9x3y^  16  3y5 
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„  ,5 

7  -  "l  ^  +  "2 


1  :  '  9  +  n  '9  ' 

^  ^  8x^9y9 


■  d'^  '3^ 

ic  - r  +  n/r  — r 

^  9x3y^  °  9y^ 


3^  3^ 

‘7  9j^b  ^8  3x'^9y 


^9  ^10  9x23y3 


3x^9y' 


3x3y 


—  +  n,  «  — r 

^  12  9y5 


SLg  -  ni3  —  +  ni4 


***"'■  4“  X1 1  ^  *1-  X1 1  ^ 

3x^9y  3x‘*9y2  9x99y9 


3^  9^  9^ 

^1"^  9x^9y‘*  ^1®  3x9y^  1^  9y^ 

The  computer  program  for  calculating  the  constants  , 

K.?g  and  the  coefficients  of  the  above  operators  is  described 
in  Appendix  II  together  with  the  results  of  a  sample  calculation 
for  a  five  ply  plate. 

The  coefficients  of  the  uncorrected  differential  equations 
(34)  for  an  "isotropic",  five  ply,  0-90  laminate  are  given  below. 
These  coefficients  were  calculated  by  the  computer  program  PLY. 


DIFF.  EQ.  COEFF.  ARE  (multiply  by  E6) 

(listed  in  decreasing  order  of  X-derivatives) 

FIRST  EQUATION 

1.57975  1.33807  E-6  0.06  1.39452  E-6 

3.56595  E-6  0  -1.48389  E-15  -2.54365  E-10 


0.442785 
•1.84575  E-15 


SECOND  EQUATION 

-5.64515  E-8  0.442785  5.3573  E-6  0.06 

1.57975  2.87856  E-17  -2.39822  E-IO  -6.69603  E-15 


THIRD  EQUATION 
9.47848  E-2  7.6897  E-8 

3.01671  E-2  7.49352  E-7 

-1.20668  E-3  -1.24926  E-8 


3.01671  E-8  3.21438  E-7 

9.47848  E-2  -2.62481  E-3 

-1.16658  E-3 


8.92325  E-6 
■7.02698  E-10 

8.36712  E-8 
■1.8733  E-9 


The  coefficients  of  the  correction  operators  given  by  equa¬ 
tions  (43)  are  shown  in  the  following  table*  These  were  computed 
by  the  computer  program  DIFFCO. 


COEFF.  OF  CORRECTION  OPERATORS;  DELTA  L(l); 
(MULTIPLY  BY  E6) 

(LISTED  IN  DECREASING  ORDER  OF  X-DERIVATIVES) 


FIRST  EQUATION 
-0.190196  1.12672 

E-2 

-2.42698 

E-2 

-4.51758 

e-4 

-1.83759  E-4 

4.37318  E-3  -3.07968 

E-2 

2.65778 

E-3 

-2.38774 

E-2 

-1.83943  E-7 

-7.22392  E-3  6.70288 

e-4 

-2.2455  E 

:-3 

6.30183 

E-5 

-6.03174  E-4 

-4.37178  E-9 

SECOND  EQUATION 
-1.1275  E-9  -5.60855 

E-2 

5.52693 

E-3 

-4.68619 

E-2 

-3.46724  E-3 

-3.67503  E-4  4.14684 

E-4 

-1.08149 

E-2 

-1.10292 

E-2 

-0.078595 

-1.20691  E-10  -2.19585 

E-3 

3.46082 

E-4 

-2.77002 

E-3 

-3.05428  E-4 

-1.90015  E-3 

THIRD  EQUATION 
-4.22148  E-3  2.4875  1 

E-4 

-0.002003 

-2.06286 

E-5 

-1.60055  E-3 

-1.35644  E-4  1.51649 

E-6 

-1.32916 

E-4 

6.70463 

E-5 

-8.88607  E-4 

-4.26205  E-4  -2.87753 

E-3 

-1.89911 

E-4 

1.52315 

E-5 

-3.4868  E-4 

1.19751  E-5  -1.30294 

E-4 

-1.95585 

E-5 

-7.73107 

E-5 

Some  of  these  coefficients  should  be  identically  zero  for  a  bal¬ 
anced  laminate  but  are  not  computed  as  such  due  to  roundoff  errors. 
In  order  to  make  a  valid  comparison  between  the  coefficients  of 
equations  (34)  and  the  correction  operators,  the  coefficients  of 
the  correction  operators  should  be  divided  by  the  square  of  the 
plate  length  since  the  order  of  the  correction  operators  is  two 
higher  than  that  of  the  coefficients  of  equations  (34) .  It  can 
be  seen  that,  even  for  an  "isotropic",  balanced  laminate,  some  of 
the  corrections  are  not  negligible. 

The  corrected  plate  theory  leads  to  a  set  of  governing  par¬ 
tial  differential  equations  of  eighteenth  order  in  the  three  dis¬ 
placement  components  U,  V  and  W.  Once  these  three  functions  have 
been  determined,  the  stresses  can  be  found  from  equations  (18)  and 


(20) .  The  transverse  normal  and  shearing  stresses  are  included 
and  all  of  the  stresses  in  the  plate  are  determined  to  high 
accuracy . 

From  an  asymptotic  point  of  view,  the  solutions  for  U,  V  and 
W  fall  into  two  classes;  slowly-varying  " interior" solutions  and 
rapidly-varying  "edge-effect"  solutions.  The  interior  solutions 
correspond  to  particular  solutions  of  the  differential  equations 
in  terms  of  the  in-plane  variables.  The  edge-effect  solutions 
correspond  to  the  complementary  solutions  which  are  rapidly- 
varying,  a  fact  which  ma^;  be  determined  by  inspection.  Equations 
(36) -(38)  disclose  that  the  rapidly-varying  solutions  for  U,  V, 

W  result  in  U  ,  V  ,  W  terms  which  are  of  order  h  .  In  the  edge 
zone,  therefore,  all  of  the  stresses  are  of  the  same  order  of 
magnitude . 

In  order  to  apply  the  high  order  plate  theory  to  practical 
plate  problems,  it  will  be  necessary  to  develop  the  proper 
boundary  conditions  from  a  variational  principle.  Expressions 
for  the  strain  energy  and  work  done  by  the  external  forces  are 
thus  required. 
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SECTION  VIII 

WORK  AND  ENERGY  EXPRESSIONS 
FOR  THICK  PLATE  THEORY 

In  accordance  with  the  principle  of  virtual  work,  the  varia¬ 
tion  of  the  quantity  u  -  v;  must  vanish  where 


S  is  the  area  of  the  upper  lateral  surface  of  the  plate  and  B 
is  the  area  over  which  the  edge  forces  are  applied.  The  dependent 
variables  of  the  problem  are  the  three  "displacement”  components, 
U,  V,  W.  The  first  task  is  to  develop  the  expressions  for  the 
strains  which  follow  from  equations  (35) .  This  can  be  expedited 
by  writing  the  displacements  U*,  V*,  W*  in  the  form 


u* 

-  x*^u 

+  X*2V  +  X*3W 

V* 

-f  X*2V  +  X*3W 

(45) 

It 

* 

w 

-  X3^U 

^32^  ^33^ 

where  the  are  the  appropriate  differential  operators  in  the 

variables  x  and  y  and  integral  operators  in  z  with  piecewise 
constant  coefficients.  The  term  1^  A^^  in  W  is  of  higher 
order  and  has  been  dropped.  The  strains  can  now  be  written  in 
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the  form 


au  , 

ax  ”  ^  9x2 


2  rZ  r 

9V  9^W  .  9  -  * 

^  "  sv2  J,h 


W  *21“  +  *22''  + 123"  It 


1  r  i£  4.  iX 

2  V  9y  9x 


-  z  -L^  +  ^  [  Jl_lx*  u  +  X*  V  x*,w 

9x9y  2  J_j^\9yj'll  12  +  13 


■*■  4  ^21^  +  ^22^  ^23^ 


I  1*1“  t  I12''  t  *13” 


+  t  ^  +  X*  V  +  X*  W  dt 

2  j 9x  31  32  33 


I  >*!«  +  >■22''  +  Im" 


-If 


X?J^U  +  X*  V  +  X*  w 


131°  t  t  x‘3W 


The  strain  energy  per  unit  plate  thickness  is  thus 


_i  +  _i  +  (  ;i*  -1. 

'xx  9x  ^xy  9y  ^xx  J  '^11  9x 


+  Vy  J_h  '21  ^  ^  ^xy 


.*  ^  \  /14. 

hi  ^  hi  ^ 
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*  f  ^  *  3  * 

+  T  X,.+T  7;^ —  dt  +  T 

X2  11  xz  J_j^  31  3x  yz  21 


+  X*  dt  +  T,,  X*,  6U 

yz  J_j^  31  9y  zz  31 


9  ^  9  7*  9 

^yy  9y  ^xy  9x  ^xx  ^12  9x 


yy  J  ,  "22  9y 


TT^  dt  +  T. 


'  *  8  *  3 ' 

.^12  9y  ^22  ^ 


^  ^xz  ^12  +  "xz 


^32  9x  "yz  ^22 


(47a) 


yz  32  9y 


^32  9v  ■'■  "zz  ^32 


i  2  2  2 

+  i  _z  4  i_  +  9^ 

i  ”  9x2  "xy  8x9y  "yy  9y2 


.  *  9 

"xx  J_i^^l3  9x 


^yy  l_h  ^23  ^  "xy  |_J^  (^13  9y  ^23  w) 


"xz  ^13  ■'■  "xz  ^33  9x  ■'■  "yz  ^23 


+  "yz  1"^  ^  "zz  f 


From  equations  (35)  ,  (44)  and  (45) ,  the  work  done  by  the  surface 


pressure  is 


-  P  6W  + 


X*^  6U  +  X*2  6V  +  X*2  6Wj  dz  dS 


(47b) 
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The  components  of  the  stress  vector  on  the  edge  of  the  plate  are 


■^1  ^XX  ^xy  3J!, 

R  =  "  IZ  3x 

T2  "^xy  8J!.  ”  "^yy  3  5. 

T®  =  7  ^^  -  7 

■^3  ^xz  dz  yz  TT 

where  Z  is  the  arc  length  measured  along  the  edge  of  the  plate 
in  the  middle  plane  and  considered  positive  when  the  plate  peri¬ 
meter  is  traversed  with  the  area  on  the  left.  A  bar  over  a  stress 
component  denotes  an  applied  stress.  The  work  done  by  the  edge 
force  is  now 
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S  in  equations  (47a)  and  (47b)  must  be  transformed  by  means  of 
Gauss's  theorem.  After  this  has  been  done,  the  full  expressions 
for  the  stresses  in  terms  of  U,  V  and  W  can  be  inserted  into 
the  variational  principle. 
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SECTION  IX 


FORMULAS  FOR  INTEGRATION  BY  PARTS 


Gauss's  theorem  is: 


(^1  ^  ^2 


Various  formulas  for  integration  by  parts  can  be  obtained  from 
this  theorem  by  allowing  F^  and  F2  to  be  products  of  two  func¬ 
tions.  For  example f  if  Fj^  =  fu  and  F^  =  0,  we  have 


J  ^  (fu)  ds  =  y  fu  dA.  However,  beginning  with  the  mixed, 
S  L 

second  order  derivative/  it  is  found  that  two  or  more  forms  can 
be  obtained  for  the  surface  integral.  Taking  =  f  ~  and 

=  l^u  for  example,  leads  to  the  expression 

f  f  ff  3  i  /f  II  +  M  u  dJi 

^  j  3x9y  9x9y  y  ■  j  9y  9<i  9x  9  jy 


Taking  F^  =  u  and  =  -  f  leads  to  the  expression 


u")  ds 

f  9x9y  9x9y  J 


^  9x  9S,  ^  3)1 


In  order  to  show  that  both  forms  are  correct,  it  is  necessary  to 
use  the  fact  that  the  integral  of  a  perfect  differential  taken 
around  a  closed  curve  vanishes. 


d(fu)  = 


+  lii  +  u  3  ^  di  =  0 

9A  ^  9y  9«.  )  ^  l9x  9J!.  9y  d  Si  J 


hence 
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f 

L 


Si!,  3y 


di 


The  number  of  possible  forms  for  each  of  the  deritatives  appearing 
in  equations  (60)  is  shown  in  the  following  table. 


rr  a 

Number  of 

(m,n)  in  f  — - — -  ds 

Jj 

possible  forms 

i 

(1,0) 

1  i 

(0,1) 

1 

(2,0) 

1 

(1,1) 

2 

(0,2) 

1 

(3,0) 

1 

(2,1) 

3 

(1,2)  ' 

3 

(0,3) 

1 

(4,0) 

1 

(3,1) 

4 

(2,2) 

4 

(1,3) 

4 

(0,4) 

1 

From  the  enpirical  evidence  presented  in  the  table,  the  number  of 
possible  forms  for  the  mixed  derivative  of  order  m+n  is  m+n. 
The  following  table  lists  all  of  the  combinations  of  functions 
which  can  be  used  to  evaluate  the  mixed  derivatives  through  the 
fourth  order. 


Integrand 

^1 

^2 

^  S^u 

8U 

^  SxSy 

*  w 

ax  ^ 

S^f 

Sf 

£  au 

SxSy  ^ 

”  ay 

"  ax 

. 
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Integrand 

f 


axay" 


u 


i  f  ^ 


ax^  ay 


a**^ 

ax^ay 


u 


a‘^u 

axay^ 

a‘»f 

axay^ 


;a^'u 


a2f 


axay  '  axay 


u 


af 

ax 

-  f 


au 

ax 


a^u 


33f  ^ 

01.  ou  o  a. 

ay  ax  3x3y 

3x2 

ax23y 

a^n  _  af  3^ 

.  u 

^  axay  ax  ay 

3X2 

3  3u 

:  _  af  au 

1 “  ay  ay 

f  a2u 
a'xay 

a2f 

ay2 


u 


11 

ay  ay 


-  f 


-  f 


a^u 
ax  ay 

a^u  _ 
axay 


111-  u 
axay 

af  au 
ay  ax 

u 

axSy 


f 

f 


,  3 

3u 

8u 

a  3f 

ax23y 

^  axay 

ax 

3x2  yy 

3f 

32u 

+ 

3u 

a^i 

ay 

ax2 

axay 

3x 

3x2  3y 

—  g3- 

U 

af  "^a^u 

a3f  . 

3x2 

ay  ' 

ax  3X3y 

3x2  3y 

"33 

u 

3f  a 2u 

—  4. 

32f  u 

3x2 

ay 

3x  axay 

3x2  y 

u 


u 


u 


u 


3x2  ay 2 

3x23y2 


U 


af  a£u 
ay  3y2 

3  2f  ^ 

9y 

- 

. ay^ . 

f  iiu 

ay  3 

f  ^  + 

— axay2.- 

^  a^u 

f  ■—  “ 

ax3y2 

af  32u 


-f 


af  8  2u 
ax  al?^' 

-  f 

ill 

3x2  -Sx 

a  3f 

3x3 
3  3u 
ax3y" 

3  3u 


'u 


ax- 


u 


32f  au^.  3  3f 


axay 


af  a2u 


y  axay'^ 
a3f 


u 


32f  au 

axay  ay 

,  _ 

ax  3y 


82u 


ay  ax3y 
'a^ 


ay' 


au 

ax 


a  3f 

axay2 

33f 

3x3y2 


u 


ax  ay  2 

ay  axay 

axay 

2 

-f 

a  3u  ^ 

af  a 2u 

32f 

au 

ax  ay  2 

8y  axay 

ay  2 

ax 

3f 

32u  _ 

a2f  3u  ^ 

3 

•  U 

3X 

..  '3y^ . 

axay  ay 

ax3y2 

11 

+  a3f 

u 

ax  axay 

3x2  ay 

32f  ^ 

+ 

u 

3x2  ay 

8x23y 

! 

^  a  3u 

32f 

3u 

3x2  3y 

3x3y 

3X 

u 


u 


-f 


3  3u 
3x2  ay 


'll 

ay  ax2 


80 


since  the  natural  boundary  conditions  are  to  be  obtained  by  setting 
equal  to  zero  the  coefficients  of  the  varioiis  derivatives  of  (SU, 

6V  and  fiW,  the  form  taken  by  these  coefficients  is  extremely 
important.  It  seems  clear  that,  in  order  to  obtain  correct  results, 
the  various  possible  forms  of  the  results  of  the  application  of 
Gauss's  theorem  should  be  averaged.  This  procedure  leads  to  the 
following  set  of  formulas  for  integration  by  parts. 

I  *  '  -  II H  “  f  * "  If 

S  S  L 
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f  ^  as  . 


3x3y 


,2 


^  as  ^  1 

3x9 y2  3 


/  3.^U 

\9y^ 


IX  _  2 

dZ 


3  ^u  3x 


3x3y  dZ 


+  M.  iH.  + 

9x  3y  3i; 


3f  /3u  3x 
3y  \^9x  dZ 


9u  3y 
9y  9Jl 


f  ds 


f  ^  ds 


9^f 


3x 


O  !X  _  2  - 

9y2  9£  3x3y  BA 


U 


d5, 


3^f 

ax'* 


9^f 


3u 


u  ds  +  (j 
L 
3 


3x2 

3^f 
ay** 

2 


i_£ 

3x^ 


u 


u  ds  -  (p 


111 

3y2 


9u 

Sy 


aff 

3y  ^ 


u 


3^u  _  ^  a^u 

3x3  9x  9^2 


9^u 

3y  ^ 


3f 


3  u 


3x 

9  £ 


3y  8y2 


d£ 


ax^ay 


ds  = 


a\ 

3xay' 


3^f  1  1 

— - —  u  ds  +  -r 

ax^3y  ^ 


3^u 

3x23y 


a^u  3x 
3x3  3  £ 


+  _?  3^u  ^ 

Txl9x2  3£  3x3y  a£ 


3f  3^u  3y  a^f 

3y  93j2  3^  !S'«r2 


+  2 


3^f 


Su 


3y. 


ds  = 


3x3y  ax  3£ 
--  a^f 


ax 

3 


3u  3y  3u  3x 
3y  3£  “ax  3£ 


ix' 

3£ 


3^f 


ay  ^  3^f 
3x2 3 y  3£  9x3 


3x 

a£ 


d£ 


3x3y- 

a^u 

3x3y^ 


U  ds  + 


L  L. 


'3^u  3y 
^ay3  3£ 


ax 

3£ 


,  ^  3  3x 
3x  9y2  a£ 


af  /, 

3y 


3  3x 
3xay  a£ 


3^u  3y\ 
ay2  ^z) 

„  a^f  au 

3x3y  3y 

3x  3^f  / 
3£  3y2  \ 

^3u 

,3y 

ix 

»£ 

3u  3x'\ 

,  o 

-1-  'i  -  -  ll 

3x  a^f 

r 

ax  aty 

“  J  _  u 

3x3y2 

3)i  9y3 

3£ 

L. 

d£ 
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f  —  ds 


^  ^  u  ds  +  ^ 


2f 


9  Xi  9y 


9x9y 


2  9  Jl 


9^u  ^ 

9x29y 


9f 

9x 


9^u 

9y2 


_2 

9Jl 


9^  9x 


9x9y  dZ 


9f 

9y 


2  2 
9  u  9y  __  9  u  9x 

9x9y  gjj2  9  5, 


9^f  9u  9X 
3x2  8y  3^ 


+  O.  ix  + 

3y2  3x  95, 


9^f 


iE  iz  - 

9x9y  \9y  95,  9X  95. 


) 


-  2 


\9x9y^ 


9^f 

9x29y 


9x 

95 


U 


d5 
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SECTION  X 


BOUNDARY  CONDITIONS 


By  applying  the  preceding  formulas  for  integration  by  parts, 
the  boundary  integral  part  of  the  variational  principle  gives  the 
boundary  conditions  on  the  edges  of  the  plate.  Results  have  been 
derived  for  a  general  form  of  curved  boundary  for  which  the  direc¬ 
tion  cosines  of  the  outward  unit  normal  are  where 

ds  is  the  element  of  arc  length  of  the  boundary  curve  measured  by 
traversing  the  boundary  with  the  interior  of  the  region  on  the 
left.  The  coefficients  appearing  in  these  expressions  are  defined 
by  the  following; 


n 

ijkJl 

rz 

- 

-h 

(t  * 

s^  Cj^^  ds  dt 

-h 

*n 

i  jkji 

*  fz 

.n  *  d 

^  ^kJt  ai’ 

P^ . 

1  jk£ 

P 


n 

i 


’  ^  n  * 

s  I .  ds  dt 

-h  -h  ^ 


P 


*n 

i 


dt 


d 


P 


n 

i 


ic 

Stresses  written  with  a  bar  above  are  applied  stresses  on  the  edges 
of  the  plate. 


6U: 
rh 


-h 


+  T  I^ 
XX  Z2  5 


J  ds 


1  * 

T  +  i  T  I, 
xy  2  zz  6 


*  Stresses  are  designated  by  the  symbol  Tj^  j . 


9x 

3s 
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*7^157.3  -^1^633)  iil  +  -ar  KV532 


p  ^  +  p 

^5522  ^  ^5632/  9s 


1  ^*o 


*o  *o 

^5532  -^5612 


5612 


+  i  P**^  ~  P*°  ) 

^  ^5523  2  5633J  9s 


9y  /,  *o  .  *o 
Ts  ~  r5533  ^5613 


9t 

yz 

3X 


2  (P6?22  ^6632)  H' 


1  „*0  1  „*o  \  9x  , 

1  ^6523  ^  1  ^6633)  9s 


.  1  p*o  .1  p*o  \  ^ 

7  ^6523  7  ^6633)  ^ 


p  4  p 

^6532  ^6612 


p  ^  +  p 

^6532  ^6612 


w  o  ^  o 

p  ^  4.  p 

^6533  ^6613 


3CX 

9x2 


2  P? 


5522 


+  2  P 


5632 


oMy 


4  p®  4  p®  4  ^  p® 

^5523  ^5633  7  ^6 


5  /  9s 

■2 

'‘t 

XX 

9x9y 


2  P^  +•  2 
^  ■^5532  ■  ^  ■^5612 


2  p®  4  2  P*^ 

^  ^5532  ^  ^  ^5612 


4  p°  4  p®  4  ^  p®  )  .iX  -  ^ 

^5523  ^  5633  7  ^6/  9s  1 


i-  iP*^  +  p*^ 

7  V  5533  ^5613 


pO  4  pO  \  ^ 

■^5533  ^  ^5613/  9s 


i  h° 


5522 


5632 


6532 


_ M. 

9x2 


1  T^O 

+  -s- 


2  f P°  +  P°  )  ^ 
^  1^6522  ^  ^6632/  9s* 


1 

+  4  p^- 


6612  2  6523  2  6633 


9x8y 


^  f p®  4  p®  4  p®  4  p®  4  i  p® 

7  1^5522  ^5632  6532  6612  7  6523 


7  ^6633  ^5)  9^  “  7  (^5532  ■*■  ^5612  7  ^5633  7  ^5523 

7  ^6,533  7  ^6613  7  ^5)  Ts  T  7  (^5532  ^5612 


■*■  (^5533  ^5613 


S13  -  I  ^0  -  ^  [  J  ('■5532  ^ 

^5523  ■*■  1  ^5633  ■*■  2  ^6533  ■*■  1  ^6613  ■*■  2  ^e)  Ss' 

2 

\  iiil  4  ^  '^yy  r^-  i/p°  +  p°  ^  H.  ' 

y  8s  3x2  P  \  6522  ^  ^6632y  8s 


8x8y 


ifp°  + 

3  I  6522 


P° 

^6632  J 


_ ^ 

8x2 

\  lx  _ 


^2  P^  +  2  P°  +  P° 

^  ^6532  ^  ^6612  ^6523 
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+  pO  +  pO \ 

^  ^6633  ^  ^5  /  8s 


+  P 


5 
h 

-h 


lx.  -I 

3s  ' 


2 

3  T 


3y^ 

,o 


1  r  ^6532 ^'''  ^  ^6612  ■*■  ^6523  ^6633 


p'^  +  p''  +  i  p*^ 

6533  ^6613  26 


T  iX  _  T  — 
XX  3  s  xy  3  s 


dz 


3x 

3s 


dz 


6V; 


+  i  T 

2  zz  6 


lx  - 

3s 


T  +  T  I. 
yy  zz  4 


1^ 

3s 


*o  +  p*o  ^  ix 

5523  ^5633  j  3s 


*o 


'o 


P""  +P''  +—  P 

^5521  ^5631  ^  2  5533 


1  *o  \  3x 

^  I  ^5613j  ^ 


3T 


xz 


-  2 


3y 

\  3x 


P*o  +  p*o  +  i  P*o  +  i  p*°  \  lx 

5521  5631  2  5533  1  5613  /  3s 


P  +  P  1  _ 

*^5531  ^561lj  3s 


] 


1  r.*0 


3T 

+  — XZ 
3x 

\  3x 


[^6523  ^6633j  ts 

3T 


+  p*v.  lx  - 


+  P  4.5.P''  +i.p''  \ 

^6631  ^  ?  ^6533  I  ^6613  3s 


6521 


XZ 


3y 


*o 


p  +  p 

6521  ^6631 


+  1  p*o 
^  I  ^ 

,*0 


+  i  p**^  _  9  / P*^  +  P*^  ^  lx 

6533  2  6613  )  8s  “  \  6531  661lj  3s 


2 

3*^1 


XX 


3x' 


5523 


+  P  ''  +  i  P® 

^5633  2  ^6 

2 

3^T 


+  P 
4 


3x 

9s 


XX 


9x3y 

.O 


iz 

9s 


1 

3 


^  ^5521  ^  ^5631  ^5533  ^5613 


^  ^5521  ^5533  ■*■  ^  ^5631  ^5613  ^4 


■J  (^5531  ^561l] 

rf^i523  -  ^ 

.o 


\  ^ 


2 

3  T 


lx. 

3s 


XX 


9y^ 


1.  ^ p^  4-  p 

3  1^^5531  5611 


lx 

3s 


/pO  +  pO  ^  lx  -  1 

1^^6523  ^6633^  3s  3 


^  ^6521  ^  ^6631  ^6533 


+  P''  +  P^  +  P^  +  P^ 

5523  ^5633  ^  ^6613  6 


3X 

3s 


2 

3  T 


xy 


3x8y 


2 

1 


(2 


4.  2  P 

6521  ^  ^  -^6631 


p®  +  p®  4.  p®  +  p®  4.  p® 

6533  ^6613  ^5523  ^5633  ^6 


+  p^  +  p^  +  i  p^  +  i  p^ 

^  ^6531  ^  ^6611  ^  1  ^5533  ^  1  ^5613 


lx  -  2 

3s  T 


(^s 


+  p 

5521  ^  ^5631 


2 

3^T 


:i£x 


+  p^  +  p^  +  p^  +  i  p^  +  —  p^  +  p^ 

^5631  ^  ^6531  ^6611  ^  1  ^5533  2  ^5613  ^4 


3y^ 

lx 

3S 


'2  /pO 

T  (^5521 
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+  p° 
5531  ^  ^56 


“  2  (^5531  ^561l)  ^  ^  J  (^6523  ^6633)  H 

"  T  (^6523  -■*■  ^663^  F  “  I  (^  ^6521  ^  ^6631  ‘ 

3  ^ 

-  ^?613  Y  -  -T;F  [I  (^'  ^^521  -  2  .  P° 

■*■  ^6613  +  T  ^e)  F  “  (2  ^6531  +  2  ^6611  ^  ^4)  ^ 

=  [  *  r  T  iz  -  T  ”1  dz 

j_j^  i  xy  as  yy  3s^ 


P  +  p^ 
6631  ^6533 


f  l^il  dz 

yy  3s^ 


3T  ,  aT  *K  - 

xy  ,  1  zz  -_*\  9x 

8x  2  ax  6 )  as 


aT  T  aT  *\  ^  /8T 

xy  1  zz  ay  _  ^ r  yy 

“ay"^  ^  1  -^6  as  ^  av 


.  °"zz  T.*\  ax  . 

'W  ^i)  7^ 


a  ^T 

ax  .  ^  xz 


*1  +  p*^  ^  -  —  ( p*‘^  +  p*^  +  p*^  +  p*^  '\ 

5522  5632 j  as  3  V  5532  ^5523  5612  5633)  as 


axay 


*1  ^  ^  _  i  /'p*! 


5532  "^5523  ^5612  ^5633  as  3  V  5521  '^5533 


+  p*l  '  +  p*l  \  i2i 
^5631  ^5613  )  as 


2 

^  ^xz  \~2  /„*!  ^  ^*1  j.  13*1 

“ay2“  3  I  5521  5533  5631 


*1  \  ay 


5613/  as 


2  +  P*^  l2E. 

^  r 5531  ^5611) as 


yz  r,  /p*l  4.  p*l  \  iX 

^  (^6522  ^6632  )  ^ 


I  (  VI 


+  p*i,A  + 


3  I  6532  ^6523  "6612  6633/  as 


*1  *1  *1 
■'■  ^6523  ■*■  ^6612  ^6633 


SxSy 


^  r4  /  *1 


1  6532 
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^  f  p*^  +  p*^  -  p*^  +  p*^  V  ^ 

1  1^6521  ^6533  ‘  ^6631  ^6613;  9s 


.  Vz  2  f*l  .  _*1  ,  _*1  ,  _*1  V  ^ 

^  ■  a„2  7  V  6521  ^6533  ^6631  ^6613  7  Ts 


9  v*^  4-  D*1  V  3x 

"  ^  1^6531  ^6611  j  Ts 


^  ^5522  ^  ^^32  ■'■  ^57  Is 


^  ^5532  ■*■  ^  ^5523  ^  ^Ll2  ^  ^^33  ^6  7  3f 


9x"^  9y 


T  r  ^5532  ^  ^5523  ■*■  ^  ^5612  ^  ^5633  ■*■  ^6  )  ^ 


7  r  ^5521  ■*■  ^  ^5533  ''■  ^  ^5631  ■‘‘  ^  ^5613  ’'■^4  7  3s 


9x9y 


1(^5531  ^5611  j  3s 


1  _  1  _  1 
„2  2  V  ^5521  ^  ^5533  ^  -^5631 


^  ^5613  ^4  y 


i  fpi 


■S'  r  5531  ^5611  7  3s 


2  P 


U: 


i  [ 

•%  I  r*  r 


6522  "6632/  2  V"  5522  "  5632  6532 


+  p^  +  ™ 

*  ^  t’  /  r\  — 


6523  "6612  "  6633  5y  9s  dx^dy  *2  5522  5632 


3 

3  T 


iJ  (pI, 


+  P^  +  P^  +  P^  +  P^  +  P^ I  -  fp^ 

^6532  ^6523  ^  ^6612  ^  ^6633  ^57  3s  V  5 


j. 

5532  5523 


+  P^  +  P^  +  P^  +  P^  +  P^  +  P^  +  P^'^ 

5612  5633  ^6521  6533  6631  6613  67  3s 


9x9y' 


1  1  1  1  1  1 
P  +  P  +  P  +  P  +  P  +  P 

^5532  ^5523  5612  ^5633  ^6521  ^  ^6533 


4-  4-  +  p^  I  —  f  p^  +  p^  +  p^  -i-  p^ 

^  ^6631  ^  ^6613  ^6  7  9s  2  1^5521  ^  ^5533  ^  ^5631  ^5613 


.  „1  .  „1  .  i^lA  9x 

^6531  ^6611  ^4/  Ti” 


5613  ■  "6531  •  ■"6611  4  7  3s 


xy  i  1  f pi  +  pi  +  pi 

,3  i  2  r  5521  ^5533  ^5631 


1\  9y 


2  P 


5531  "56117  9s 
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1  /pi  +  pi  ^  2E 
1  V6522  ^  ^6632  ds 


[7  ^^6522 


+  p 


1 

6632 


7  r^6532  ^^6523  ^^6612  ■*■  ^  ^6633  ^5)  3s 


3^  r 

I  (^^6532  ^  ^6523  ■*■  ^  ^6612  ^  ^6633  ■*■ 

I  6^6521  ■*■  ^  ^6533  ■*“  ^  ^6631  ^  ^6613  3s' 


+  2  P^  +  2  P^  +  2  +  P-1 

2  ^6523  ^6612  ^6633  5  3s 


i  \  3: 


-  yy  1  / 2p^  +  2  p^'  +  2p^  +  2P^  +  p^  1 

3y3  T  i^^^6521  ^6533  ^^6631  ^^6613  ^6^  3s 

—  f  2P^  +  2P^  +  P^^  dz  =  T  ”  T 

^6531  ^  ^^6611  ^  ^4  ]  3s  I  xz  3s  -^yz  3S 


Other  boundary  conditions  can  be  written  which  involve  the  coeffi¬ 
cients  of  the  following  variations: 


9 

3  3 

3 

3  I 

6U 

Tx  ' 

9y  *’  3x^  '  3y^ 

.2  .2 

9x9y  1 

,>2  “1 

V  \J 

3 

9  9^ 

9 

9 

3x  ' 

9y  '  9p^ 

f  ^ 

3x3y 

U  V 

.2 

.2 

.2 

.3 

3 

9  9 

9 

3 

9 

3x  ' 

8y  '  3x2  '  3y2  ' 

3x3y  ' 

3x^ 

There  are,  in  all,  22  natural  boundary  conditions.  Further,  all 
of  the  natural  boundary  conditions  are  inhomogeneous.  These  are 
all  extremely  complication  expressions.  Because  of  the  length 
and  complexity  of  this  corrected  theory,  its  practical  application 
seems  highly  questionable  and,  for  this  reason,  the  remaining 
natural  boundary  conditions  will  not  be  written  out  in  full. 
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SECTION  XI 


DISCUSSION 


The  Correspondence  Principle  results  are  exact  within  the 
limitations  of  the  small  deformation  theory  of  elasticity.  How¬ 
ever,  the  application  of  these  results  will  be  for  a  finite  number 
of  laminations.  It  was  shown  in  Section  IV  that  good  agreement 
between  the  Correspondence  Principle  and  "exact"  results  can  be 
obtained  with  a  relatively  small  number  of  isotropic  plies.  In 
other  cases,  the  error  depends  on  the  ply  properties  and  the  lami¬ 
nation  scheme.  In  general,  the  error  depends  on  functions  of  the 
type 


fZ 

F*  (2)  f  (z)  dz 

■'-h 


Lim. 

n->“ 


(z)  dz 


where  F*  (z)  is  a  piece-v^ise  constant  function.  If  the  mean 
value  F  defined  over  the  entire  plate  thickness  is  identical  with 
the  mean  value  defined  over  a  comparitively  small  number  of  plies, 
then  good  agreement  should  be  expected.  Considerations  of  this 
nature  limit  the  application  of  the  Correspondence  Principle  to 
plates  which  are  constructed  using  a  lamination  scheme  which  con¬ 
sists  of  a  repeating  subgroup  containing  a  number  of  plies  which 
is  small  compared  to  the  total  number  of  plies  in  the  plate.  This 
rules  out  the  use  of  the  Correspondence  Principle  for  unbalanced 


laminates . 

The  attempt  to  derive  a  plate  theory  corrected  for  transverse 
normal  and  shearing  deformation,  described  in  Section  VII,  although 
technically  successful,  results  in  a  set  of  partial  differential 
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equations  of  high  order  and  an  extreraaly  complex  set  of  boiandary 
conditions  on  the  edges  of  the  plate.  The  difficulty  of  solving 
this  problem  in  any  practical  case  renders  this  attempt  a  practical 
failure.,  Corrections  to  the  classical  plate  theory  are  important 
only  in  the  neighborhood  of  the  edges  of  the  plate  where  the  stress 
distribution  is  three  dimensional.  The  starred  corrections  to  the 
displacement  field  in  equation  (35)  represent  the  rapidly  varying 
part  of  the  displacement  field  near  the  edges  of  the  plate.  The 
high  order  corrections  to  the  differential  equations  (34)  which 
resulted  from  the  inclusion  of  these  additional  terms  contribute 
additional,  rapidly  varying,  edge- effect  solutions  to  the  problem. 
In  the  edge  zone,  the  displacement  corrections  to  equations  (35) 
in  terms  of  U*,  V*,  W*  are  of  the  same  order  as  the  respective 
displacements  U,  V,  W, 

An  alternate  approach  to  tiie  development  of  asymptotic  solU“' 
tions  to  the  stress  field  near  the  edges  of  a  laminated  plate  is 
to  employ  the  procedure  recommended  by  Reiss  and  Locke  (Ref,  3) 
for  the  solution  of  the  generalized  plane  stress  problem  in  homO“ 
geneous  media.  Briefly,  the  recommended  procedure  would  be  to 
perform  a  "stretching  transformation"  on  the  coordinate  normal  to 
the  edge  of  the  plate  in  equations  (4)  and  then  to  develop  regular 
asymptotic  solutions  in  powers  of  the  thickness  parameter  h  .  The 
greatest  advantage  of  this  procedure  is  that  it  would  give  the 
edge  solutions  directly  without  having  to  resort  to  the  simultane-’ 
ous  solution  of  a  high  order  set  of  differential  equations. 


91 


92 


and  foj 


1,4x10 


Fig.  2.  comparison  of  C  and 
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APPENDIX  I 


"BASIC"  COMPUTER  PROGRAM  FOR  THE  CALCULATION 
OF  EQUIVALENT  Mi4TERIAL  CONSTANTS 

This  program  allows  the  rapid  calculation  of  the  equivalent 

material  elastic  compliances  d. .  and  the  equivalent  coefficients 

J 

of  thermal  expansion  ej  defined  in  Section  II,  It  is  written  in 
"BASIC"  language  and  compiled  and  run  on  the  Dartmouth  College 
Time-Sharing  Computer  System. 

The  input  data  is  inserted  sequentially  in  the  DATA  statements 
898  through  910.  The  meaning  of  the  symbols  used  in  the  program 
is  given  in  the  following  chart. 


SYMBOL  MEANING 

El,  E2  Longitudinal  and  transverse  moduli  of  a  single 

ply;  E^,  E^ 

G1 ,  G2  Shear  moduli  of  a  single  ply;  ^ET^  ^TT 

Nl,  N2,  N3  . Poisson's  ratios  of  a  single  ply; 

Tl,  T2,  T3  .  Thermal  expansion  coefficients  of  a  single 

ply; 

M  ...............  Number  of  plies 

Q(K)  . . Principal  direction  of  each  ply  relative  to 

the  X  axis,  degrees,  (K  =  1  to  M) 

H(K)  . . .  Thickness  of  each  ply,  units  arbitrary 


(v^m  represents  the  contraction  in  the  longitudinal 
LT 

direction  due  to  an  applied  stress  in  the  transverse 
direction. ) 

The  program  operates  in  one  of  three  modes,  depending  on  the  value 
of  the  integer  control  number  G0  in  data  statement  898.  If  G0 
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is  0,  the  program  computes  the  compliances  and  thermal  expansion 
coefficients,  using  the  Correspondence  Principle.  It  is  up  to  the 
user  of  the  program  not  to  insert  data  defining  an  unbalanced 
laminate  for  this  part  of  the  program  as  the  results  will  be  mean¬ 
ingless.  If  G0  is  1,  the  zeroth,  first  and  second  order  moments 
of  the  Cj^j  are  computed.  These  quantities  are  needed  to  define 
the  edge  force  and  moment  relations  in  terms  of  the  stress  func¬ 
tion  components,  and  to  evaluate  the  coefficients  of  the  equations 
to  be  solved.  If  G0  is  2,  the  program  computes  and  prints  the 
individual  it  and  thermal  expansion  coefficients  for  each 

ply. 

Several  assumptions  have  been  made  in  writing  this  program 
which  should  be  pointed  out.  First,  the  program  assumes  that  all 
plies  are  of  the  same  material  and  that  the  ply  orientation  angle 
is  the  only  variable  which  changes  from  one  ply  to  the  next. 
Second,  it  was  assumed  that  the  coefficients  of  equation  (9)  are 
constrained  by  the  relations  which  appear  in  statements  35,  45,  55 
and  64.  This  was  done  for  expediency  since  data  giving  these 
coefficients  was  not  available.  All  of  these  assumptions  can  be 
easily  changed,  if  desired. 

The  results  for  a  five  ply  plate  are  shown  below.  The  data 
shown  in  the  listed  program  in  lines  898  -  910. 
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PLV 


PLIES 

r  thic'h"^esses  are 

S02  00^3  3.0a  3o33 

r  AMGLSS  ARE 
93  3  90  n 

-IPLfAMCES  AHE  (MULTIPLY 


.96997 

.955AI 

.25 

.19339 


~J.955Ai  E-2 
R.0699T  F>’8 
-0.55 

»S  .■6,0729  E--6 
0 


•6  SO  IM/LB) 
•!  .90619  E-a 
“!,00659  E-2 

1  o 

Sc0Bll3  E-7 


^HMAL  GOEFFIC  lEivJTS  ARE  (V^VTS  ARE  PER  DEG) 
*065'“iQ  E-T3  U06502  E-5  0*00003 


-1*29663  E“7 
”6.7SS07  E“6 
“^*12569  E-13 

I  O 


"6*52283 


ZEHOTH  MOMEWS 
2*57975 
0«3327f35 
3*56595  E'"6 
FIRST  M0MEMT5  C 
-7.02698  E»10 
-1  E-10 

-1  .8^1575  E-1  5 
SECOMD  MOMENTS 
iol6658  E-3 
a.709B7  E-4 
2«-^i9616  E-9 
ZEHOTH  MOMENTS 
-1  .97468  £-2 
FIRST  MOMENTS  C 
8.3844  E-12 

MOMENTS 
-2.97722  E-ti 
THIRD  MOMENTS  C 
4«33265  E*“14 
ZEROTH  MOMENTS 
-2.09013  E-5 
FIRST  MOMENTS  C 
9oi6281  E-15 


lCIaJ)  AHE  (MULTIPLY  BY  E6> 
.382785  1 c07146  E-5 

.57975  2.78904  E-6 

.64515  E«8  Co  ^12 

ARE  <^?ULTIPLY  BY  E6> 


5o64515  E«8 
C(I^J)  ARE 
1.70303  E-iO 
)ol4n42  E“10 
2.87856  E-17 
^  C(I^J)  AHE 
5.47596  E-4 
>.62481  E-3 
95161  E-ii 


-4 o 85329 
-I o5!268 
-6o95195 
(MULTIPLY 
7.50023 


E-1  5 
E-!  5 
E-I  1 
BY  E6) 
E-9 
E-9 


95161  E-ii  0.306144 
HI)  AHE  (MULTIPLY  BY  E0 ) 
i. 97468  '£“2  6*  17334  E-8 

HI)  ARE  (MULTIPLY  BY  E3) 
).357B5  E-l-2  -2.-97613  E-17 

^  HI)  AHE  (MULTIPLY  BY  E0) 
i. 76203  E«5  4.32133  E-ll 

HI)  ARE  (MULTIPLY  BY  E0 ) 
H39253  E-14  -5*84453  E-20 

'  FCHJ)  ARE  (MULTIPLY  BY  E6) 
H09013  E-5  2o95504  E-11 

FCHJ)  ARE  (MULTIPLY  BY  E6) 
1*22263  E-i4  -lo33835  E-20 


DIFF.  EQ.  COEFF.  ARE  (MULTIPLY  BY  E6> 

(LISTED  IN  DECREASING  ORDER  OF  M-DER I VATI VES > 

FIRST  EQUATION 

1 .57975  1.33807  E-6  0«06  H39452  E-6  0.4427B5 

3.56595  E-6  3  -1*48389  E-15  -2.54365  E-10  -1,84575  E-15 

SECOND  EQUATION 

-5.64515  E-8  0.442785  5.3573  E-6  0.06  8.92325 

1.57975  2.37856  E-17  -2.39822  E-10  -6.69603  E-15  -7.0269B 

THIRD  EQUATION 

9.47348  E-2  7.6897  E-8  3o0i67S  E-2  3.21438  E-7  8.36712  E-B 

3.01671  E-2  7.49352  E-7  9.47848  E-2  -2.-624‘81  E-3  -1  .8733  E-9 

-1.20668  E-3  -1.24926  E-8  -1.16658  E-3 

LOADING  COEFF.  ARE  (DIMENSIONLESS) 

(LISTED  IN  DECREASING  OliDEH,  OF  X-DER I  VAT  I VES  ) 

FIRST  EQUATION 
6.16034  E-3  -1.34343  E-B 

SECOND  EQUATION 
-1.84343  E-B  7.00422  E-3 

THIRD  EQUATION 

-2.22735  E-4  5.96756  E-9  -HGBiOl  E-3  -1 
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r 


w  w 


PLY 


2  READ  Ci^ 

3  REM  IMSERt  DATA  LT\"ES  898  THRU  910  •  MODULI  IM  LIME  900  SHOULD  BE 

4  REM  IN  MILLIONS  OF  LB/SQ  IN  •  ALL  OTHER  DATA  IN  NORMAL  UNITS  . 

5  READ  El#E2irGl>G2>NWN2>N3#TUT2^T3 


10 

LET 

AM  W  ) 

=  1 

./El 

15 

LET 

A(2#2) 

s 

1./E2 

20 

LET 

AC3^3) 

= 

A(2.2  ) 

25 

LET 

A(4#4) 

=s 

1 ./C2.*G1  ) 

30 

LET 

AC5#5) 

1 ./<2.*G2) 

35 

LET 

A(6^6) 

=: 

A<5>5) 

40 

LET 

AC1#2) 

=: 

-N1 /El 

45 

LET 

A<1*3) 

= 

An*2) 

50 

LET 

A(2^i  ) 

5= 

-N2/E2 

55 

LET 

AC3i.l  ) 

S 

A<2^  1  > 

60 

LET 

AC2#3) 

3 

-N3/E2 

64 

LET 

A(3^2)? 

“AC 

2^3) 

65 

READ  M 

66 

FOR 

K*1  TO 

M 

67 

READ  Q(K) 

68 

LET 

PCK)-Q 

CK) 

69 

NEXT  K 

70 

FOR 

o 

11 

M 

71 

READ  H(K) 

72 

NEXT  K 

73 

LET 

74 

FOR 

K*1  TO 

M 

7  5LET  H«H-^H<H) 

76  NEXT  K 

77  FOR  K«1  TO  M 

78  LET  Q<K>*C6. 2832/360* )*Q(K) 

79  NEXT  K 

80  FOR  Kal  TO  M 

81  LET  Q»aCK> 

85  LET  Q1«SINCQ) 

90  LET  Q2=C0SCQ) 

95  LET  Q3=Qit2 


100 

LET 

Q4**Q2t2 

105 

LET 

Q5«ai t3 

110 

LET 

Q6==Q2t3 

115 

LET 

07=01 t4 

120 

LET 

Q8=Q2»4 

125 

LET 

A0*=ACl#a 

5  + 

AC2j  1  >4.2a*Af4i 

^4) 

130 

LET 

A1=A(!^ ! 

>'t* 

AC2ji2>-2.^A 

iC4j 

p4) 

135 

LET 

A2^fK(l»2 

)“ 

ACl^l  )'»'AC4, 

4  ) 

140 

LET 

A3=AC2j2 

)- 

A(2>1 )-A<4# 

4  ) 

145 

LET 

A4=AC3>2 

)- 

AC3^  1  > 

150 

LET 

A5=AC2^ 1 

)- 

AC  1  >  1  ) 

155 

LET 

A6=AC2,2 

)- 

ACl >2) 

160 

LET 

A7=AC2j3 

)- 

AC  1^3) 

165 

LET 

AR=AC 1 , 1 

1  + 

AC2j2)-AC 1, 

2)-A(2*  1  ) 

170 

LET 

A9=AC6^6 

AC5^5) 

175 

LET 

3Cli.l  )=Atl 

, 1 )*Q8>AC2> 

2  )*C7-*-A0i>Q3-*aA 

180 

LET 

3(1#2>=A 

Cl 

,2)^Q8-*-AC2^ 

1  )« 

••Q7+A1  •*Q3*Q4 

135 

LET 

3CU3)  =  A(1 

#3  )*Q4+AC2^ 

3)  *03 

190 

LET 

BCU4  )  =  2 

•  ★ 

A2*Ql*Q6+2. 

*A3*Q5*Q2 

191 

LET 

3C 1 j5)=n 

« 

192 

LET 

BCN6)  =  0 

a 
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i9S  LET  BOi' Slcg''*Q?-- 'iCa^  S  >=»QS->-AI *Q3*Q4  '  ■ 

8550  LET  B(2*a)«A<Cf,jS  >=."Q'J:AC£o2>*Q-^"^A0*Q3*Q4 
803  LET  BCSi.3)~ACl  ,3?*D3-:.ACa^3>*Q4 
81®  LET  BC2*4>aao‘>A3oC,{*g.S{'8o*.<.2*QS«‘Qa 

ass  LET  BC2,5)=Rc,  .  •  .  . 

8ia  LET 

a!3  LET  BOi.  1  )»AC3..t  S  i>Ae3/-2>*g3 
280  LET  B<3..2)sAC33  S  )*Q3':'AC3;>2J*Q4 
825  LET  QO»3>=A?3;3S  , 

83®  LET  B{3i.4  )sS„*,Vr’>OS-sGS 

83S  LET 

232  LET  B<3*6)'-®o 

835  LET  B  (4  ,  I  ) sAS *0 1  ‘^'Gg  <•  A6  AQ 5  aQ'2 ❖A « 4  « 4  )  *Q  1  *02 <•  (  1  .  -8  .  *03  ) 
243  LET  B  <  4  *  2  J  s  A  6  *  Cl  i  ^-O  6  <  A  S  -^e  5  *  QS  «  A  C  4  s  4  )  *  Q 1  *Q2  ♦  <  1  .  -2  .  *Q3  ) 
845  LET  BC4»3>=A7*G:il*Qa 

850  LET  8(4*4  )s2<.*Af;':‘Q4*Q3  >ft':4/>4  >*<  S  <.“3.*Q3)f3 

ass  LET  3<4i.5)=0  = 

252  LET  B(4i.6)=0<> 

aS5  LET  D(5i.5)eA(5^S)*Q«-:^At:6a65*Q3 

860  LET  B(5*6)*A9’:<QS*Ga 
261  LET  B<5*l)>s3. 

868  LET  B(5*2)=0= 

263  LET  B(5*3)«ffo 

264  LET  B(5*4)a;0c 

265  LET  B(6*5)=BCS.-6> 

87®  LET  B(6>6)=AC5i>5>'^-e3-:-*<6i-6i*Q4 
27!  LET  B(6* 1 )«0c 
278  LET  B<6.a>=3o 

273  LET  B<6*3)«0<. 

274  LET  3(6*4)000 

275  LET  T(1  >=Tl*Q4-^-T£*(33 
a7B 

280  LET  T(a)«Tl*Q.3-:’T2«Gi4 

SSS  LET  T(3>=T3 

287  LET  T(4)=«(T2-Ti  )*Q|::^C!£: 

290  LET  Dl=B(l*8)*(BC2.9»  5*3(4p43>‘“B(2*4>*B(4*1)) 

295  LET  D2oB(I*S  )*(BC2**gi‘::BC4-.4)-B(2*4)*B(4*2)) 

300  LET  D3=D(l»4)'-’--(BC2ii2.>'-B(4*i  >-B(2*l  )*3(4*2)) 

335  LET  D=D1-D2+D3 

310  LET  C(  1  *  1  )o(0(  S  S  )«B(  !*  S  )*B(4*4)  )/D 

315  LET  C(2*  1  )  =  (BC  f.  *2  5“S(  S*4  >*B(4*2)  )/D 

320  LET  C(3*  I  )aCBC  !/>  S  )=?'BC4<!a>  =B(  1*2)*BC4*  1  ))/I) 

325  LET  C  (  1  *2  )=(B(2.- S  )--B<'4ii4  )“B(2*4  )*B(4*  1  )  )/D 
330  LET  C(2*2)o{B(a^<<>*:-*BCiK'>a)‘-B«2*2>*B(4*4))/D 
335  LET  C  (3*2  )a(3(2*a  S  ) -B(2»  1  )'«B(4»2>  )/D 

340  LET  C  (  1  *3  )=(B(  S  *  I  )*BCS.-;4  5-BC  f.  *4  )  *B(2*  I  )  )/D 
345  LET  CC2*3)=(B{2*a)*BCS*4)«BiS*2)*B(2»4))/D 
350  LET  C(3»3  )=(BC  S  iig  1  )-B(  2*  2  )*B(2*2)  )/D 

355  LET  1(1  >=3(  I.  .:-2  ) 3  )  *C(  1  *  1  )-*'B(4*  3  >*C(  1  *3  ) 

360  LET  I(2)=D(!*3)*C;C2..2>-.'B(2*3?*C(2*1  )+3(4*3)*C(a*3> 
365  LET  1 (3 )  =  B( 1 *3 )*C(3*2  ><B{2*3  5*C(3* 1 )+B(4*3>*C(3*3  ) 
370  LET  1(4  )=B(3*  2  )‘i>C(2.'  2  3-:EK3^25*C(1*  1  )+B(3*4)*C(3*  1  > 
375  LET  I(5)oB(3*I  )*CC2;.2i-iB(3p2>*C(l*2)+B(3.4)*C(3*2) 
380  LET  I(6)=B<3*  1  >*0(2333  :-B(3.8)*C(l*3)+B(3*4>*C(3*3) 
385  LET  I(7)=,B(3*3)-B(3/-  2  >  *  !  (2  )-B(3*  2  )  *I  (  1  )-B(3*4  )*I  (3  ) 

390  LET  D4  =  -B(1*2)*<BC8*,2  J*B(4i>43“3{2*4)*0(4»l  )> 

391  LET  D5s8(  1  *  I  >*(BC2*2  )*ic!<4*4  )“B(2*4  )*B(4*a  )  ) 


S3 


392 

393 

394 

395 

396 

397 
39H 

399 

400 

401 

402 

403 

404 

405 

406 
410 
412 
414 
416 
418 
420 
422 
424 

426 

427 

428 

429 

430 
432 
434 
436 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 
4  50 

451 

452 
454 
4  56 
4  57 
458 
4  60 
462 
464 
466 
468 
470 
472 


LET  D6=-B( 1 ,4)«(B(2^2)*BC4^ 1 )-BC2^ 1 )*BC4>2)) 
LET  D7=D4+D5>D6 

LET  F1=T(1 )t<DC2# 1 >*BC4#4)-BC2>4)^BC4^ 1 ))/D7 
LET  F2=-T(2)*CB< 1 # 1 >*B(4#4 )-3( l#4)*Bf4^ I ) )/D7 
LET  F3«TC4  >*<3C1  j  1  )*B(2*4>-B<  W4)i'B(2^  I  )  >/D7 
LET  F1-F1+F2+F3 

LET  F2=-TC1 )*C3<2>2)*BC4#4)-B(2#4)^BC4>2))/D7 
LET  F3=TC2)*(B<  U2)*BC4^4)“BC  U4)*B(4>2>>/D7 
tET  F4--TC4  )t(BClj»2)i'B(2^4)-B(  1#4)*B(2^2)  )/D7 
LET  F2=F2+F3+F4 

LET  F3=-TC 1 )*<BC2>2)*BC4^ 1 )-B(2# 1 )*3C4#2) )/D7 
LET  F4  =  T(2>*C3C  U2)*BC4^ I >-BC 1 ^ 1 )*3(4#2 ) )/D7 
LET  F5=-T<4 )*CBC 1>2  >+B(2^ l)-B( 1^ 1 )*BC2#2> )/D7 
LET  F3=F3+F4+F5 

LET  I  C8)==3C3>2)^F1+B(3>  1  >*F2-3C3^4  )*F3+TC3  ) 

LET  JCl  )=3(3>  1  )'♦!  C2) 

LET  J(2)=BC3>2)*I C 1 ) 

LET  J(3)=BOi4)*IC3> 

LET  J(4 )=BC3^ I )*F2 
LET  J(5)=3C3^2)*Fl 
LET  vJC6)=3(3>4  >*F3 
LET  JC7)=F1 
LET  J(B>=F2 
LET  JC9)=F3 
IF  G0=2  THEM  920 
IF  G0-1  THEM  430 
IF  G0=0  THEM  600 
LET  F(3>=-H/2. 

LET  FCK)=F<K-i )+HCK> 

FOH  J=1  TO  3 
FOH  lal  TO  3 

LET  XCI  j  J)  =  XCI>  J)+Ca^  J)*HCK> 

LET  YCI^  J)=yCI^  J)+CC<  I#  J)  +  CFCK)t2-FCK"l  )t2>  >/2 

LET  ZCI#  J)=Z<n  J)  +  CCa^  J>*(F<K>t3-FCK-l  )t3)  )/3 

MEXT  I 

NEXT  J 

FOR  I«1  TO  3 

LET  Ed  )=:E<I  )  +  I<I  )*HCK) 

LET  GCI  )=Gd  )  +  (Id)*(FCK)t2-FrK-l)t2)>/2 
LET  KCI  )  =  KCI  >  +  CICI  )^(FCK)t3-F(K«-1  )  t3)>/3 
LET  NCI  )  =  Md  >-»“<Id  >^(F(K)t4-F(K-l  )t4))/4 
NEXT  I 

FOR  1=7  TO  9 
LET  LCI  )=Ld  )-»»Jd 

LET  MCI  )  =  Md  >  +  C  JCI  >’^CFCK>t2-F(K-l  >t2>  )/2 
NEXT  I 


IF  M-K=0  THEN  460 
GO  TO  680 

PRINT  •’ZEROTH  MOMENTS  OF  CCI#J>  ARE  (MULTIPLY  BY  E6>  *• 
PRINT  XCld  )^X(1j2)^XC1>3) 

PRINT  X(2d  >>XC2>2),XC2,3) 

PRINT  XC3d  )>X(3,2)^X(3>3) 

PRINT  •VFIHST  moments  OF  CCI^J)  ARE  (MULTIPLY  BY  E65  ” 
PRINT  Y(ld)#Y(lj2)>yci^3) 

PRINT  Y(2d)^Y(2>2)#Y(2^3) 
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474  PRINT  YC3i.  1  )*Yf3^2).Y<3>3) 

476  PiUMT  "SECOND  MOMENTS  OF  C(I»J>  ARE  (MULTIPLY  BY  E6)  " 

478  PRINT 

480  PRINT  2(2. 1  )i.Z(8j>2)>Z(a*3) 

482  PRINT  Z(3^ !  )^ZC3/2),Z£3»3) 

484  PRINT  "ZEROTH  MOMENTS  OF  1(1)  ARE  (MULTIPLY  BY  EO)" 

486  PRINT  EC  1 E(2 E<3  ) 

488  PRINT  "FIRST  MOMENTS  OF  Id)  ARE  (MULTIPLY  BY  E0)'* 

490  PRINT  Gd  >^G(2)>G(3) 

498  PRINT  "SECOND  MOMENTS  OF  Id)  ARE  (MULTIPLY  BY  E0)" 

493  PRINT  Kd  )jiK<2)^K<3) 

494  PRINT  "THIRD  MOMENTS  OF  Id)  ARE  (MULTIPLY  BY  E0)" 

495  PRINT  Nd  )»N(a>>N(3) 

496  PRINT  "ZEROTH  MOMENTS  OF  FdjJ)  ARE  (MULTIPLY  BY  E6)  ’* 

498  PRINT  L(7 )>L<8)#L<9) 

500  PRINT  "FIRST  MOMENTS  OF  Fd»J)  ARE  (MULTIPLY  BY  E6)  " 

502  PRINT  MC7 ),M(3)#M(9) 

503  DIM  U(41 ) 

534  LET  H=H/2 

505  LET  Ud  )  =  X(2>2) 

506  LET  U(2)  =  X(3>a)+?t(S^3)/8 

507  LET  U(3 )=K{3»3 )/8 

508  LET  U(4 )=X(3>3 )/2 

509  LET  U(5)  =  X(2d  )+X(3^3)/2 

510  LET  U(6)  =  X(3> 1 ) 

511  LET  J(7)=YC2>2) 

512  LET  U(8)=Y<3>2)+y(2*3) 

513  LET  U<9)=Y(2d  )+Y(3>3) 

514  LET  Ud0)=Y<3^  1  ) 

515  LET  U< 11 )=X<3,2) 

516  LET  Ud2)=X(d2)+X(3>3)/2 

517  LET  Ud3)=X(d3)/2 

518  LET  Ud4)=X(3^3)/2 

519  LET  U(15)-XC3^  l)+X(d3)/2 

520  LET  Ud6)  =  Xd>l  ) 

521  LET  Ud7)=Y<3»2) 

522  LET  Ud8)=YCd2)+Y<3^3) 

523  LET  Ud9  )  =  y(3j.  1  )+Yd  >3) 

524  LET  U<20)=Yd  j  1  ) 

525  LET  U<81 )=H+XC2>8)-Y(2»2) 

526  LET  U(22)=2*H*X(3^2)-2*Y(3>2)+H*X(2>3)/2-Y(2»3)/2 

527  LET  U(23)=HtXd*8>-Yd>2)+H*X(3»3)-Y(3>3) 

523  LET  U(24)  =  (H*XCd3)-y(d3))/2 

529  LET  U(25)=(H*XC2^3)-YC2#3))/2 

530  LET  U(26)=H*X(2» 1 )-Y(2^ 1 )+H*X(3>3)-Y(3>3) 

531  LET  U(a7)=2*H*X(3>l  )-2*Y(3/ 1  )+H*Xd»3)/2-Yd>3)/2 

532  LET  U(28)=H*^X(  d  1  )-Yd  >  !  ) 

533  LET  U(29)=H*Y(2>2)-Z(8>.8) 

534  LET  U(30)=2*H*y(3»2)-8*Z(3/2)+H*Y(2#3)-2(2*3) 

535  LET  U(31  )=H*y(2>  I  )-Z(2#  1  )+H*Yd  *2  )-Z  (  d 8  )+2*Hi'Y(3* 3  )-2*Z(3>3  ) 

536  LEI  U<3a)=2fH*Y(3»l)-2=i'Z(3^1)+H*yd»3)-Zd#3) 

537  LET  UOSlsH'i'YCl^  I  )-Z(  d  1  ) 

538  PRINT 

539  PRINT  "DIFF.  EQ.  COEFF.  ARE  (MULTIPLY  BY  E6)" 

540  PRINT  "(LISTED  IN  DECREASING  ORDER  OF  X-DERI VATI VES )" 

541  PRINT  "FIRST  EQUATION" 
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542  FOR  I*!  TO  10 

543  PRINT  U(I>i 

544  NEXT  I 

545  PRINT  "SECOND  EQUATION" 

546  FOR  1=11  TO  23 

547  PRINT  U(I)# 

54 B  NEXT  I 

549  PRINT  "THIRD  EQUATION" 

550  FOR  1=21  TO  33 

551  PRINT  U<r)^ 

558  NEXT  ! 

553  LET  UC34>  =  -CEC2)/4+Gt2)/f2>t'H)  +  K<2)/{4*H»2)) 

554  LET  UC35  )  =  -<E(3  )/4+GO  )/<2=^'H)+K<3  )/<4*Ht2  )  ) 

555  LET  U(36)  =  -CE(3)//4+G(3)/(2*H)+K(3)/C4*H»8)) 

556  LET  U<37)  =  -(E(1  )/4+G(l  >/C2+H)+K(l  J/CA'i'Hta)) 

557  LET  U<35)=(H+EC2)+G(2)-K(2>/H-NC2)/<H»2)>/4 

558  LET  U(39)=2*<H*E(3)+GC3)-K<3)/H-N(3>/tH»2)) 

559  LET  U<40>  =  H*E{1  >+G(n+K(l  >/H“NCn/(Ht2) 

560  LET  U(41 )=-l 

561  PRINT 

562  PRINT  "LOADING  COEFF.  ARE  <DIMENS lONLESS )" 

563  PRINT  "(LISTED  IN  DECREASING  ORDER  OF  X-DERI VATIVES )" 

564  PRINT"FIHST  EQUATION" 

565  PRINT  U(34),U(35) 

566  PRINT  "SECOND  EQUATION" 

567  PRINT  U(36)>U(37) 

568  PRINT  "THirsD  EQUATION  " 

569  PRINT  U(33>^U(39),U(43>>U(41 > 

598  GO  TO  999 

600  FOR  1=1  TO  3 
605  FOR  J=1  TO  3 

610  LET  R(l^ J>=<C(I^ J)*H(K)>/H  +  R(I/J) 

615  NEXT  J 
620  NEXT  I 
625  FOR  1=1  TO  8 

630  LET  S(I)=(I(1)*H<K>)/H  +  SCI) 

635  NEXT  I 

640  FOR  1=1  TO  9 

645  LET  Wd  )  =  <JCI)*H(K)>/H+W<I) 

650  NEXT  I 

655  LET  V(5x5>=CB{5#5)*H(K>)/H  +  V(5*5) 

660  LET  V(6»6>=(B(6*6)*HCK>)/H  +  V(6j6) 

665  LET  V(5»6)=(B(5i.6)4rH(K))/H  +  U(5j6) 

670  LET  V(3,3)=CB(3>3)*H(K))/H  +  V(3,3) 

680  NEXT  K 

685  IF  G0>0  THEN  999 

690  LET  C1=R<2,1)*(R<1»2)*R{3j3)-R(3*2)*R<1*3)) 

695  LET  C2=R(ld  )*(R(2*2>*R(3»3)“R(3#2)>fH<2>3)) 

700  LET  C3=R<3,1)*<R<2*8)*RC1,3)-RC1*2)*R(2^3))  , 

705  LET  C=C1-C2+C3 

706 

710  LET  D< 1> i >=(R(3d )*R< I<3)-R<1» 1 )*R<3»3) >/C 
715  LET  D(1,2)=(R(2>  1  )*R(3>3)-R(3^  1  )<=r<C2>3))/C 
720  LET  D<1*4)  =  (R<1/ 1  )*R<2^3)-R<2> 1 )*R(1>3))/C 
725  LET  DC2, I )=(R(li2)*R(3.3)-H<3,2)*R(l#3))/C 
730  LET  D(2#2)=(R(2>3)*R(3,2)-a(2#a)*R(3>3) )/C 
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735  LET  D<2^/i)s!<H(2,2)*R{I,3)-;Ul#2)*H(a,3))/C 
740  LET  DC/Ij.  I  )  =  (R(  1  j  1  )*;H3,2)-R(  1»2  )*R(3*  1  >)/C 
745  LET  D(4ji2  )  =  (R(2,2)*;U3»  1  1  )*R<3,2)  )/C 

750  LET  D<4 j4  )  =  CH (2> 1 ) *RC1 <2 )-K< 1  *  1 )*R(2, 2 ) )/C 
755  LET  D(  1  j3  )=S(  1  )  i-DC  I  ,2  ) +S(2)*n<  1  #  1  )+S(3  )*Dn  #4  > 
760  LET  D<2>3  )=SC1  )*D<2>2>+S(2)+D(2*  1  )+S(3)+r)Ca,4) 
765  LET  D<4>3)=  S f 1 ) *0(4^ 2 ) +S (2 ) *D(4j 1 ) +S C3 ) *D<4»4  ) 
770  LET  D<3i  1  )=S(4)i'DC2,  1  )+S<5)*D(l<  1  )+S(6)*DC4, 1  ) 
775  LET  2  )  =  S(4 )*D<2> 2 )+S ( 5)*D< I #2 ) +S< 6)*D(4 >2 > 
780  LET  DC  3*  4  ) -S  (4  )  ♦DC2>4  )+Sf  5)'»'DC1  >4  )+S  (  6  )*D<4j4  ) 
785  LET  DC5,5)=VC5#5) 

790  LET  DC6^  6)  =  V(6i.6) 

795  LET  DC5j.  6)  =  V(5^  6) 


800  LET  D(6#5)=D<5,6) 

805  LET  D(3  j  3  )  =  U(3#  3  ) +D(3i  1  )  *S<2 )  ■•■D<3j  2  )*S(  I  )+D(3*4  )*SC3  )- 

806  LET  E<3)=SC8)-Dt3#l)*WC3)“D<3/2)*WC7)+D<3»4)*W<9) 

807  LET  D1 =DC2> 1  )*DC4>4  )-D(2>4 )*D(4 j 1 ) 

308  LET  D2=DC 1 ^ 1 ) +DC4 #4  )-DC 1 j4 ) *D(4j I ) 

809  LET  D3=DCl,l)*D(2,4)-DCi>4>*DC2#l) 

810  LET  D4=DC2>  2  >'!'DC4  »4  >-D(2#4  )*D(4  j)2  ) 

811  LET  D5=D< 1 > 2 >*DC4 »4 )-DC 1 ,4 >*DC4 j2 ) 

812  LET  D6=-D  C  1 , 2  )=^D<2/4  )-D(  1  >4  )*D<2>2  ) 

813  LET  D7  =  DC2^ 2  )*DC4^ 1 )-DC2> 1 )*D(4 j2 > 

814  LET  D8=D(  1/2)*DC4>  I  >-D<  1>  1  ):*:D(4j8) 

815  LET  D9=DC1^2>*DC2j I )-D( 1 )*D(2j2> 

816  LET  B1=D1 *CD5*D9-D6*D8) 

817  LET  B2=D2’^(D6*D7-D4*D9) 

818  LET  B3=D3*(D4>t:D8-D5*D7 ) 

819  LET  Dl=31+D2+B3 

820  LET  B2=-DC1>2)*D1 


W(  n-W(a)-W(3  ) 


821  LET  U3=D( 1 , 1 )*D4 

822  LET  B4=DC1,4)*D7 

823  LET  B2=B2+B3+a4 

824  FOR  1=7  TO  9 


825 

826 

827 

828 
829 

831 

832 

833 

834 

835 

836 

837 

838 

839 

840 

841 
850 
852 
854 
856 
858 
860 


LET  y( I )=W( I )*B2/B1 

NEXT  I  -  , 

*  ’“'''^^’*^*^^*^'^"’^*5*D3)+yC8)t'CD2%D9-D3*D3>+W(9  )*CD3*D5-D2*D6) 
LET  E(8 )=W(7)+CD4 *D9-D6*D7)+WC8) *CD1 *09-03*07 )+W(9)*CD3*D4-Dl *06) 
LET  K<  4  )  =  Wt  7  )  *(04*08-05*07  )+V;(8)*(Dl  *08-02*07  >+yt9  >*<02*04-01  *05) 

PRINT  m;”PLIES"  V* 

PRINT  "PLY  THICKNESSES  ARE" 

FOR  K=1  TO  M 
PRINT  HCK); 

NEXT  K 
PRINT 

PRINT  "PLY  ANGLES  ARE” 

FOR  K=1  TO  M 
PRINT  P(K); 

NEXT  K 
PRINT 

PRINT  "COMPLIANCES  ARE  (MULTIPLY  BY  E-6  SQ  IN/LB)" 

PRINT  D<l>l),D<l,2),D(l,3),D(l,4) 

PRINT  0(2^1), 0(2^2), DC2/3)> 0(2, 4) 

PRINT  0(3, 1 ), 0(3, 2), 0(3, 3), 0(3, 4) 

PRINT  0(4,1 ), 0(4, 2), 0(4, 3), 0(4, 4) 

PRINT  0(5, 5), 0(5, 6) 
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862  PRINT  D<6#5)*D(6^f ) 

864  PRINT  ••THERMAL  COEFFICIENTS  ARE  CUNIT3  ARE  PER  DEG)^' 
866  PRINT  E<1)/E(8>>E(3)>E(4> 

870  60  TO  999 
898  DATA  1 

900  DATA  25./ 1 .5>  .8 
902  DATA  .05/. 25/. 25 
904  DATA  10.E-6/30.E-6/3O.E-6 
906  DATA  5 

908  DATA  0/9iT/v!/90/0 

910  DATA  .02/ .03/ .02/ .03/ .02 

920  PRINT  •’PLY  C(I/J)  (MULTIPLY  BY  E6)/  PLVJK 
924  PRINT  C(l/l>/C(l/2)/C(l/3) 

926  PRINT  C (2/ 1 )/C(2/2)/C<2/3  ) 

928  PRINT  CC3/1  )/C(3/2)/C<3/3> 

930  PRINT  ••PLY  F(I/J)  (MULTIPLY  BY  E6)/  PLY^'IK 
932  PRINT  J(7)/ J(8)/ J(9) 

934  PRINT  ••PLY  Id)/  PLY  ••JK 
936  FOR  1=1  TO  8 
938  PRINT  1(1); 

940  NEXT  1 
942  GO  TO  680 
999  END 
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APPENDIX  II 


EVALUATION  OF  THE  CORRECTION  COEFFICIENTS  FOR 
TRANSVERSE  NORMAL  AND  SHEARING  STRAIN 


The  coefficients  Kf ^ 


i jklmn 


and  are  of  the  form 


F  (z) 


G  (t) 


S^  H  (s)  ds  dt  dz 


To  obtain 

*  * 

P  =c,. 

*  * 

*  * 

mn 


and  to  obtain  ,  insert 

*  * 

P  =c.. 

ic 

G  =  1 

*  * 

«  ="k 


The  functions  F*,  G*  and  H*  are  assumed  to  be  piece-wise  con¬ 
stant  functions  of  z,  defined  within  each  ply  by  constants  P^, 
G^p  H^.  Let 

d^  =  -h  +  (hj^  ^2  *  *  *  ^i^ 

be  the  distance,  measured  in  the  z  direction,  from  the  middle  sur- 

‘th 

face  of  the  plate  to  the  top  of  the  i  ply*  Lst  h^  be  the 
individual  ply  thicknesses.  Then, 


S^  H  (s)  ds  = 


k  Hi 

'iii 


X  1“1 


q+1  k 
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where  d.  <  t  <  d,  ,, 

1  “  -  X+x 


Now, 


The  third  and  final  integration  can  now  be  performed. 


N  is  the  total  number  of  plies.  The  integrations  indicated  may 
now  be  carried  out  and  the  result  written  in  the  form 
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The  initial  indexes  have  been  changed  to  1,  2,  or  3,  as  appropri¬ 
ate  to  avoid  requiring  the  convention  that  a  summation  is  not  to 
be  carried  out  unless  the  upper  limit  equals  or  exceeds  the  lower 
limit. 
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APPENDIX  III 


COMPUTER  PROGRAMS  FOR  THE  EVALUATION  OF  THE  CORRECTIONS 
FOR  TRANSVERSE  NORMAL  AND  SHEARING  STRAIN 

The  integral  which  was  evaluated  in  Appendix  I  has  been  pro¬ 
grammed  for  computer  evaluation,  using  BASIC  language.  However, 
since  BASIC  does  not  allow  the  use  of  a  subscripted  variable  with 
eight  levels  of  subscripting,  it  was  necessary  to  regard  the 
coefficients  and  as  functions  of  their  indexes. 

The  program  is  called  DIFFCO  and  is  listed  on  the  following  pages. 
The  individual  summations  in  the  integral  I  (Appendix  I)  are 
written  as  subroutines  FNA,  FNB,  FNC,  and  FND  in  lines  910-1000, 
1010-1130,  1150-1270,  and  1300-1500,  respectively.  These  sub¬ 
routines  do  double  duty  for  the  evaluation  of  the  coefficients 
^ijkmn  ^ijk  replacing  b*^  by  1  and  C*^  by  I* 

the  latter  calculation.  The  coefficients  of  the  correction  oper¬ 
ators  defined  by  equations  (56)  are  also  calculated  by  this  program. 
In  order  to  carry  out  the  programming  within  the  limitations  of 
BASIC,  it  was  necessary  to  label  the  eight-level  and  five-level 
subscripted  coefficients  as  singly  subscripted  variables,  using 
the  re-labeling  program  MATCH,  which  is  listed  immediately  follow¬ 
ing  DIFFCO.  The  output  of  MATCH  lists  the  subscripts  i,  j,  k,  £, 
nt,  n,  p,  q  in  the  left-hand  columns  and  the  corresponding  single 
subscript  in  the  right-hand  column.  Following  this  table  is  a 
list  of  the  indexes  i,  j#  k,  p,  q  of  the  five-level  subscripted 
coefficients  and  the  corresponding  single  subscript.  The  coeffi¬ 
cients  of  the  correction  operators  m^  and  are  computed 

by  tlie  subroutines  DIFFCO  1,  DIFFCO  2,  DIFFCO  3,  and  DIFFCO  4, 


108 


which  are  listed  at  the  end  of  the  program  DIFFCO.  In  order  to 


check  the  calculation  of  a  specific  coefficient  m^^  and  n^^ 

in  DIFFCO,  it  is  necessary  to  look  up  the  single  subscript  of  the 


constants 

quantities 


„  and  K?^,  in  MATCH  and  then  locate  these 
ijkJlmn  13k 

in  DIFFCO.  For  example,  K235g33  which  is  involved  in 


the  calculation  of  £g  is  referred  to  in  DIFFCO  as  K(306). 


Using  the  data  listed  in  lines  1600-1620  of  DIFFCO  (an 


"isotropic",  five-ply,  symmetric  laminate),  all  of  the  coefficients 
^ijkdmn  ^ijk  been  computed  and  are  listed  following 

MATCH,  The  left-hand  columns  show  the  pertinent  indexes  and  the 


right-hand  colvimns,  the  value  of  the  coefficient  with  those  indexes. 
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DIFFCO 
i  READ  HS 

6  read  Eiii2,Gl,Qa.Ni,N2;N3 
10  LET  Atl.l)ai/61 

19  LET  A(2|2)«l/E2 

20  LET  A<3i3)“A(2*2) 

29  LET  A<4i4)»l/(2*Gi) 

30  let  A|5i5)«1/{2*G2) 

39  let  AC6i6)=  A(5,5)  _ _ _ 

AO  LET  A(li2)«»-Ni/El 
45  let  A(i,3S9A<i.2) 

90  LET  A(2i1)»-N2/E2 
95  let 

60  let  A(2i3)B-N3/E2 

65  LET  A(3,2)«!  A{2.3)  _ _ _ 

70  READ  Y 
75  for  X»1  TO  Y 
80  READ  Q(X) 

65  let  P(X)aO(X}  _  _  „  _ 

fO  NEXT  X  ■  ■  . . 

?9  fqr  xbi  t8  y _ : _ _ _ 

ioo  read  h(x) 

105  NEXT  X  _ _ _ 

110  LET  Ha  0 

115  FOR  Xai  TO  Y _ _ _  _ 

120  LET 

122. NEXT  X  _ 

124  FOR  Xal  TO  Y 

..126  LET  O(X)a(6e2832/360>»Q(X)  _ _ 

128  next  X 

129  let  M2a0  _  _  _ _ _ 

130  FOR  Isl  TO  3 

132  FOR  .Isl  TO  3 _ 

134  FOR  Ka  1  TO  2 

136  FOR  Lsl  TO  2  _  _ _ 

138  FOR  Hsl  TO  3 

140  FOR  Nsl  TO  3  _ _ _  . 

145  FOR  Xsl  TO  T 

150.  let  OsO{X)  _ _ _ _ 

155  LET  QlaSlN(Q} 

160  LET  Q2aCOS<0}  .  _ 

165  Let  Q35Q1*2 

170  LET  04sQ2f2  _ _ 

175  let  05aQl«3 

180  LEt  06=02»3  . . . 

165  let  07s(}it4 

19q  let  Q8sQ2f4  ,,  ,  _ _ .: 

195  let  Ala  A{l,l)*A{2,2)-2'*At4«4> 

200  let  AOa  A(1.2)*A(2»lh2*A(4K4)  _ _ 

205  let  A2a  A(i,2)“A{lil>*A(4,4) 

210  let  A3»  A(2,2)-A(2il)-AH»4) _ 

215  let  A4a  A<3i2)-A(3ii) 

220  let  A5aA(2<l)-»A{i,i)  . . . . 

225  let  A6S  A<2,2)-A(l,2i 
230  let  A7a  A52»3)-A{i,3 j 
235  let  AOa  A(l,l)*A(2i2)-Atl,2J-A(2;i> 

240  let  A9a  A{6,6>“A(5j5j 

245  let  B{t,i)s  A{lil)*aB*A<2,2j»07>AO'»034Q4 
250  let  3(1,2)*  A(i#2)»a3*At2»l)*QT*Ai*Q3«04 


255  |,ET  B{1,4>«2«A2»Q1«A6*2hA3*Q5*32 

HQ  UET  8(2. 1)s  A{1i2)‘*Q7*A(2i1)»05*A1*Q3»Q4  _ 

265  LGT  B{2.2)a  A(l» 1 )«U7*A»2. 2) »Q34AO»03«04 

2^0  UET  B{2.4)«  2*A3»Qi*Q6*2«A2«Q5*Q2;  .. 

275  let  B{3,1>«  A{3,i)»Q4*A(3,2)*Q3 
280  let  B{3,2)s  A<3il)ttU3^Al3.2i*Q4 
B(3.4>s  2*A4#0i*Q2 

0(4,1 )«  A5*Qi*02*A6*Q5*Q2*A(4,4 j*Ql*Q2*(l-2*03j 


285  LS' 


LiT 

250  LET  -  --  «-  -  -• 

295  lET  B(4,2>e  A6*0i«06  +  A5*05*02"A{4»4}*Q1*02«><1-2«U3J 

300  let  0(4. 4)«  2»A6«Q4*Q3*A(4M>«(1»2»Q3Jt2 
805  LET  B(5.5)a  A(5»5)«(J4*A(6,6)*03 

310  let  B(5.6)  =  A9*Q1»Q2  _ _ 

815  let  B(6,5)8  B<5.6) 

820  let  B(6.6)a  A(5»5)*Q3*A(6,6J*Q4 . .  . . . 

125  let  DlaB(l,2)'»(B(2,l)*BM.4>-B{2.4)*B(4»l)> 


325 

330 


LET  D2a  B(lil)*<B(2»2)«Bt^*4>'3{2*4)«B(4.2)| 
335  let  D3aB{1.4)«(8(2,2)*Bi4,l)-.B(2;i)*B(4.2)) 

340  LET  Dfl0l-D2+D3  ^  . 

345  let  C(1,1)b  tB(l,4‘)«‘S(4yi)-0(l,l)*3(4|4}  )/D 
350  LET  C(2.i)a  (B(i,2)*B{4;4)-a{j..4)*3(4.2))/° 
355  let  0(3. i}a  (B(1.1)«3(4,2)-’9(1,2>,3(4.1)  )/D 
360  let  C(1,2)s  (B(2,l)*B<4f4)-B{2i4)*3(4.1) )/D 

865  let  C(2.2)a  (8{2,4J*£'(4i2)*a{2.2)»3f  4.4)  )/0 
370  let  C|3,2)a  (B(2,2)»B{4;'l)-0{2a»»3«'*»2))/D 
375  LET  0(1, 3)a  (B(i,iJ*3(2»4)»'8{ii4)*3(2»l>)/D 
380  let  C(2.3)s  (8(2.2)»3(1.4)iB(1.2>»3(2»4))/D 
585  I  ET  C(3.3)c  (8(i,2>*8(24l)-8(lil>»3(2i2))/D 
386 


405  let  B(i.l)a  B(5»5) 

410  LET  B(l,2)n  B(5.6) 

415  let  B(2.1)3  B(5i6) 

420  U^T  B(2.2>=  B(6i6) _ 

500  let  D(0)«  -H/2 

.510  LET  D(X)«  D(X~1>*H(X)  e  v. 

530  LET  ZaZ+rNC(liJ.K,L.M.N,R,S,X)4FND{|.^»K,L.H.N,R|S,XI _ _ 

'535 

540  next  X _ ^ _ _ 

845  let  TaT*i 

550  IF.  M3>0  THEN  564  . . . . : _ _ _ ^ _ _ 

Jlo  RRiNT^"VALUeS^OF  K{I  .'jiK.L*«.N»>.Qj  ARg.  IHUCTlPL-^..®t_6^J.  ."j. 

561  PRINT 

562  RRINT  I»JaK^4iL*4»l1lN|Rj8|TA3135^jZ: _ _ _ 

‘563  IF  M3a0  THEN  580 

.564  let  K{T1*Z 
570 


580  NEXT 
590  NEXT 

.600  next 

610  NEXT 
620  NEXT 


630  next  1 
640  let  Ml8Mi*l 
650  JF  Mial  THEN  680 
660  ir  Nla2  THEN  700 
670  IP  Mla3  then  720 
675  IF  Hl>3  THEN  1650 
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fl 


A8Q  b,ET  S«i 

690  BO  TO  130  _ 

700  LET  R«1 

705  LEI  S,0  .  -  . . •  .  • 

710  00  TO  138 
720  LET  S«1 
725  60  TO  130 

910  DEr  FNA(  1  »J»K.L«M»N*«.S,X)  - - -  - - - - - 

9i2  if  H2«b  then  ^20 

914  LET  J»  J)»I(K)/(S'»1> 

916  QO  TO  930  .  - 

IIS  U?  ^ 

9«0  let  B3a{0U'-l»*<S+2)  }»<D{XH<R*l)“0<X*l>f  (B*1)  )/US*2)«|R*1  )  J 

960  lET  05a{D|X-l>*(S  +  2h<*{DIX)T(R*l)-O(X«J:>rCR*t)>:/tR*l> 

970  LET  B6a91»(E!2-B3“B4'»-B55 
972  IF  X>i  THEN  976 
974  let  CO.O 

976  let  C0»B6*C0  _ _ 

’980  let  PNA,C0 

1000  fn  end  _ _ _  .. 

1010  DEF  FNBe I  * JiK»E*M»N,R«S»XI 

1012  IF  X>1  THEN  i050  . . 

1014  LET  F3:»e 

"1050  LET  a?=ca».i)*<D^)f<Rar^^^ 

’~10i7  LET^^BSsiIkiM  DI  j()HS*-2J-DIx"’lj*<S+2)I/I  ^S*2)*(S*1)) 

1070  LET^08ab<K#L>«0(M,N)«l0{X)t(S*2}-0<X»l>MSi2))/((S*2)«(S*l)) 
1080  LET  B9s(D<X)»tS*i)  )»HCK)/^S*1) - - — - - 

1082  IF  m2«0  then  1085 

1083  LET  B9aB9«l{K)  . .  .  - 1 - 

1084  60  TO  1090 

_l0eS  let  B9,09»B(K,Lj«C«M.Ni _ _ _ 

in9o  LET  B68B0-89 
_1100  let  .F3507*F4^r3 _ : - 

1105 

-1110  LeT  F4«08*F4  . . 

1120  LET  rN8=F3 

1130  FNENO  .  ,  . . 

“itSO  OEF  FNC<  I»Jif^»E»M,N,R,S»XJ 

_1155  IF  X>1  THEN  1182 _ _ _ _ _ 

1160  LET  C3^0 

._1165  LET  C43  0  - - - - - - 

1184  lE7^CiSc{i^j7MD(X>»(R*2)*D{X»1)*IR-»2)?/IR42) 

1190  LET^Clic?l,J>*8(KiLj«(D«X)HR*2)-0<X-lM|Ri2n/(R*2) 

1200  LET  Q2sD{X^l)«{D(X)*(R+iJ*0(X-l)»{R+l) J/(R*1) 

1201  IF  M2»0  then  1204  . 

1202  LET  C2=C2*Ctl»J» 

1203  Go  To  1205 

1204  LET  C2aC2»C{I»Jj*B(KiL) 

1205  LET  QlaCl-C2  .  . . — 

1216  IF  M2«o  Then  1220  . 

1217  LET  02jHK)«<D(X)t(S*l><'DU'’i>*<S*l)|/<S*lI 
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1218  60  TO  I2w0 

1220  LET  C2,C(M»N>*(0(X)i(S,«lJ*0(X-l)*(S*l»)/{SIU 
""  1230  l,ET  (;3sCi*C4*C3 

1240  LET  C4!;C2*Crt  .  . . 

1250  lkt  rNc=(;3 
1270  FNEND 

1300  DEF  Ff4D(|*J.K,UfM,N,R|S.XI 

1310  IF  Kvl  THEN  137{)  _ _  _  . 

1320  LET  FOsO 
1330  LET  FlaO 

1340  LET  F2a0  .  .  -  . 

1370  Let  c53c<J  + 

1302  IF  M2»0  THEN  1400 

1394  LET  C68l(K)<»(D(X)t(S*l)-D(X-l)t(Sfl)f/^S*l)_  _ 
1396  GO  TO  1432 

1400  LET  C6^C(M,NJ«<D(X>t(S*lJ*0<X-l)»{S^l>7/(Sll) 
1432  IF  M2«0  THEN  i440 

1434  LET  g75H(X)  _  _  _ _ _ . _ 

1436  qO  TO  1450 

_ 1440  LkT  C7^8(K,L)»H|X)  _ 

1450  let  r0»C5»Fl*F0 

1460  LET  FlaC7*F2^Fl  _ _ _ 

■  ■  1470  LET  F2^C6^F2 

1480  LET  FND*F0  _ _ : _ _ _ _ 

"■  1500  FNEND 

__  1595  data  1  __  _ _ _ _ 

1600  DATA  25. f  1.1  15.  ,2 

1605  DATA  .05,  *25,  .25  _ _ 

■  1610  DATA  5 

1615  DATA  0»90.0i90»0  _ _ _ 

.  1620  data  i02>  .03i  .02*  .03002 

_ 1650  LET  M2,l _ : _ •• 

1652  LET  TaO 

1653  FOR  R»0  TO  1  _ _ _ _ _ ^ _ _ 

1654  FOR  S«0  TO  1 

1655  FOR  lal  TO  3  _ J _ 

1660  FOR  .Jal  TO  3 

_ .1665  FOR  K»1  TO  3  _ _ 

1670  FOR  X«l  TO  Y 

1675  LET  0«Q(X)  _  _ _ _ _ _ _ _ 

■1680  LET  QlsSINIQ} 

1685  LET  Q2sCQS(Q)  _  _ _ _ 

■'■■■:i690  LET  Q3sQlt2 

_ JL695  LET  04s02t2 _ ' 

1700  LET  Q3:,Qlt3  ■' 

l7o5  let  Q6«Q2»3  _ _ _ _ _ _ 1 _ 

■  '1710  LET  Q7«01t-| 

1715  LET  Q0«O2t4  . ^  _ _ 

’  1720  LET  A0sA(i,2)4A(2,1)*2i»A(4«4) 

_ 1723  LEy  AlaA|i,l)*A(2,2)-2»»A{4»4|  _ _ _ 

1730  LET  A2>A(l,2)-A(l,l‘)*At4,4) 

1735  LET  A38A(2.2)-A{2.1)-A«4.4)  _ _ 

«740  let  A48A(3i2>-A(3;i) 

1745  LET  a5sa(2,1)-A(1,1»  . . . . 

1750  LET  A6BA(2,2»-A(t,2j 

1755  LET  A7aA(2,3>-A(l,3)  ^  _  _  _ _ 

1760  LET  A8,A(l,l)*A(2,2j-Ati*2|-A|2.1J 

j;763  LET  A9sA{6,6>-A(5,5>  .  . 


L  l77o 
I  1779 
1780 
X789 
1790 
1795 
■  1800 
_  1805 
‘1810 
1819 
■  1820 
1825 
1830 
__1835 
1840 
-  1845 
■  “1850 
_  1855 
‘I860 
1865 
"“1870 
_  1875 
1880 
1885 
1890 

_ 1895 

1900 

_ 1905 

1910 

__l9l5 

1920 

_ 1925 

1930 

__1935 

1940 

_ 1945 

1950 

_ 1955 

1957 

_ i960 

1962 

_ 1965 

1966 

_ 1970 

1971 

_ 1972 

1975 

_ 1900 

1985 
^1990 
1995 
_  1998 
2000 

2005 

2006 
_ 2007 

2010 
.  2015 


let  8{l»i)=A(i* t )»Q84A(2,2}*Q7*Ao*Q3»Q4 

let  B{1.2)«A«l*2)«08«>A<2»l)«Q7*Al»a3«04  _ _ _ _ 

let  9tl,4)=2»A2*Ql*A6’«-2‘‘A3*Q5'‘Q2. 

let  B(2»l)fBt 1»2)*Q7^A12i1}*Q3*A1»Q34Q4  ,  _ _ 

LET  B(2,2)sA(lil)«U7>A<2,2J«Q3*A0»Q3«Q4 
LET  B(2»4)e2*A3*Qi*06*2«A24a5*Q2: 
let  C(3il>rA(3rl)*Q4*At3i2}'»Q3 

LET  B{3.2)-*A<3il)*a3*At3,2)*Q4  _  _  _ 

let  B< 3» 4) ?2*A4«qi*q2 

LET  B(4il)»A5oQl«'Q2*A6Ji05«Q2*A{4»4)«Ql*02«<l«2«Q3l 

let  8<4  *2)?A6*lJi*Q6*A5«Q5*02~A(  4»  4)*Q1*Q2*<1'«2*Q3) 

LET  B(4»4)=2*A8*Q4»03*Ai4i4)«(i*2»Q3)*2 
LET  B|5,5)RA(5,5)«04+At6,6)«Q3 

let  8(5.6 )aA9»Ql*Q2  _  _ 1 - 

LEy  B(6.5)aB(5i6) 

let  B(6.6)sB(5ii5)*Q3*A,)|6,6)*Q4  . .  .. 

LET  0l=B(li2)»(B(2ii)«3(4j45-3(2.4)«B(4#l)J 

LET  D2aB ( i>« (B ( 2» 2 )*8 '4* 4 )*3 ( 2» 4)*B(4» 2)) . . . 

LET  D35BU»4}*<B(2e2)«B(4Al>‘-3(2»l)‘»B(4,2)7 

LET  DBD1-D26D3  .  ,  , . „ - 

LET  C(l.l)s{B(l,4')«B(4»t)-8(l»l>»3(4»4T>/0 

LET  (i(2.1)s(y  (1.2)*3(4»4)^B(1i4)»3(4i2)  >/D  _ 

LET  C(3.1)f (0(1.1 )*3<A;2i*B (1.2) *3<4fl)J/0 

LET  C(i.2)?(0(2.l)*3(4;;4)*B(2»4)*3(4ii)»/0  - - 

let  C(2.2)s(8(2,4  j*3(4',;2)w8(2i2)*3(4i4)  )/D 

LET  C(3.2)s(B(2,2)*B(4a)*B(2,l)*3(4.2))/D  _ 

let  C(1.3)a(B(la)«a(2;4)-0(li'4)»3(2ilT)/D 

LET  C(2.3)«iB(2,2)«3(i;4)*0(1.2)»3(2.4n/D  _ _ 

LET.CI3.3>-(8!l,2j«B(2*l)-B(l.l)*3<2i2))/D 

let  |(i)«B(3»i>*C(2,i5>8(3r2)«C(i.  j.)*8(3.4y«C{3il) 
let  1(2)?  B<3«l)*C{2»25+B(3»2)eCCl»2)»B(3»4)«C{3j2) 

lET  1(3)4  B(3»lJ»C(2»3}*^i3»2S»C(li3J*B(3*4)«C(3i3) — 
LET  0(0)?-H/2 

LET  D(X)r.D(«-l)'*H(X| . 

LET  Z*FNA( i I J. Ki L»  MjNj  S>I)*FN3( 1 . J. K» L. M?N» R. S» X ) 

LETZ?Z*FNc(i»  J*X»L«H.NJ,R.  S«X)tFNO(  I .  J.K.L.H.N.R.S.X) 

next’x 

LEt  T«T*1  - - 

IF  M3>0  then  1972 

IF  T>1  THEN  1970 . . . .  -  —  - — 

Pr|nT  sVALUES  OF' kU.J«KjP>,!3)  are  (MOLTSPLT  BY  E6)  » 

PR  I  NT 

print  I»JjK43)R|SjTAB(30ifZ:„_ _ ; - 

IF  M3»0  THEN  1975 

LET  J(TJ?Z  - - - - - - 

NEXT  K 

NEXT  1  . 

NEXT  S._ - - - - - - - 

NEXT  R 

DIM  X(i30O}»d(iiO)«L<i6l»N<l6)»N{l9>  . _  . 

IF  M3*Q  then  4000  ^  , 

print  »c06FF,  of  correction  OPERATORS)  DELTA  LiUl" 

PRINT  n{LlSTED''ON^OECReASIN3  ORDER  0^  j^-DERIVATIVES)* 


.  L »  M.N.R.S.X) 


PRINT 

PRINT 

PRINT 

PRINT 


bFIRST  EQUATION" 


R8020  lET  HbH/2  .  , 

202S  GAtU  «OIFFGOl««K()»JU»l.(»t{H!__ 
■■'2030  FOR  T»i  TO  i6 

2035  PRINT  UJT),  ,  . .  . 

8040  NEXT  T 

8045  PRINT  . 

8050  PRINT  hSECOND  EOUaTION^ 

2055  CAl,l»  nDIFFC02wU((  )»  jnZH(UH!  . 
”  20^0  *^0^  ’“1 

’■  2065  print  M«T)i . .  . . 

8070  NEXT  T 

207S  PRINT  . 

2060  PRINT  »TH1RD  EQUATION" 

■  _  2065  GALL  '*DIFFC03'MK(  )  *  JC )  TNO«Hi 

2090  CAU  >DIFFCQ4"{K(  )»JOjN|lrH: 

"  _  8095  FOR  T«1  TO  l9  ,  ..  _ _ 

3OOO  PRINT  N<T)» 

3005  NEXT  T  . . . .  . . . 


4000  end 
5000  SUS 

6001  let 

5002  LET 
6003  LET 
5004  LET 
6005  LET 
5006  LET 
'6007  LET 
5008  LET 
6C09  LET 

6010  LET 

6011  LET 
5012  LET 

"5013  LET 
50i4  LET 
"5015  LET 

5016  LET 

5017  LET 
_5018  LET 

5019  let 

5020  LET 

5021  let 

5022  LET 
"5023  LET 
.  5024  LET 

5025  UET 
_6026  LET 
"5027  LET 

5028  LET 

5029  LET 
.  5030  LET 

5031  LET 

5032  LET 

5033  LET 

5034  LET 

5035  let 
.5036  LET 

5Q37  LET 
6038  LET 


-  . . . . . 

L(l)?-2«<l?ll4y)-*K(l6l)tJU4)/2>  ^ 

Ut2J«-2‘‘K{152)-K(150)*2*K(l55)-K|l62)-Utl5)/2-2«»K(2571  . . 

L{2)t.L(2)-2«K(269)^J<23)-<(l85)-K(l97)-d(l7>/2 

L:2)??l.{2>“'<(206j-'Kl20-1)/2-K(209)-K(2i6j[/2-2*K(i3ll'‘2»K(i43) 

L{2)?L(2)-K|3llj-Kt323}-J(l7)/2 

L(3)»-KC153)-K(i56)-2<»K{260)-<(258)-2*K(263)-K(270)  _ 

L<3)?*U(3>-K{ia8>-'K»l86)/2-K(l9l)-K(l98)/2-K(293)-'K(305) 

.L(3)-l13)-j|24)/2-J(18J/4-J(27)/2-K<207)/2-K(210)/2.-24KU34} 

L<3)iL(3l-Kli32)-2«K(i37)-.X(i44)’iX(3i4>"Kl3l2l/2"»^<3i7) 

L!3)=-L(3)-Kiy4)/2r2*K|239)*2.K(25uij<l8)/4-Jua’J<26)/2 
L(4)?-K(261)-k<261>.-K(i99)/2-<(l92)/2-K(296>-K(294)/2 
L(  4)8Lt4)’'Ri  299  j-Kv306)  /2— J{27)/4»K(l35>»K(13e)-K{315l  /2 
L«4>*Ll4)-K<3l8)/2r2*KC242)#K|240)-2*K<245)-Kt252l-J<l2l/2 

LU)*=L|4  >-J(27)/4-'Ut20l  _ _ . . 

L<5)n»'K(297>/2'*K(300)/2*K(243)’»K|246j»J(21)/2 
L|6)«-K(l50)-K<i62)-J(l5)/2»K<204)/2-K<2l6)/2  ,  „ 

L(7)B-2*R(l<8)'’'<|153>-2*K(l60)-K{l56J-K{256)«K{270)-K{l66J/2 

L(7)aL(7)-K<i98)/2rJ<13)-JI24)/2-J<l8>/4  . .  . . . 

L<7)*L<7)-K<202i-KI207}-K(2l4»-K<2lO)/2’*K«l32)-Kll44J 

l.(7)*L(7)-K(3l2»/2-K(324>/2*J(lBJ/4  .  . . 

L(  8)s-2*t'<15l  )‘'2*K  1154  j'.2»K( 256) "K( 2611  "K{268)*2 

L(8)«L(8)-K{264)-K*i84)«K(L99)/2*Kil96J"K(l92)/2-»K|294)/2 . 

L(8)8L(B)-K(3o6)/2-’J(22)-JI16>/2-J(Z7-)/4 

L{8>*L{0>-K(205)-Kl2{)9)-2«Kli30>‘T2*Kli35>-2*K<l42l  ^ _ _ _ 

L(8)aL<8>-X(l38)-Kit3lO)-K(3l5>-Kl322)-H(3lB)/2‘»K<240) 

L(8)*L<8)-K{252)-vJ{16)/2**J|l2)/2-d(27)/4 

l,,9)..2«K(259>-2«K^262)-K{l87)'Ku90)“K(292)'-K(297)/2 

L«9)3L(9)'K(304)-Kl300>/2"J(25>/2-2«K(l33)-2«X(l36)-K{S13) 

L(9)?l<9>-K(3i6)-2<»K(233)*2»K(243>-2«K<250I 

Ll9)»L(9)-K{246>-JnO)rJ(25)/2-Jl2l)/2  _ 

L(10)»-'<{295)rK(29G)-2'?<t241)-2«X{244)-J(l9) 
L(ll)#2*K{473>*2*Kt405)*J<41)*2BX(527)+2«t<{53«) 
L<12)«2eK<476)*2*Kl474i*2#K{479)*2*K| 486)*2bK(581)*2»K(593) 
L<l2)=Lli2»*l<<509):K}52l)*J<42)*d<50)*J<44)/2  .  ,  ■ 

L(l2)sL{12I*K{530>»K|528)*K{533)*Ki540F*K<455)+K(467)t*X|635) 

LU2>sL(12)*K(647)*j(44)/2  «  *  - 

L{13)»2«'<(  472>*2*K(477)^2»X<434)*2*K|480)*2*K<584)*2*K<982) 
L<l3)BL(l3)*2eKj587l*2fKI594)4K(5l2)*«510>*K(5l5)*K(322) 
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LET  LCl3)sLtl3>*K(6t7>*K{629)*j(4o)*J{5l)<.J(45)/2*J(53)/2 

let  La3)  =  L(13l*K(526)+K<531)*K(538)*K(534T*2«K(458>«2»K  456  . 

LPT  L  l3l=L{13)*2•^^(46lJ♦2»<:(469)♦K(638)♦K<636)♦K(641)♦K(648) 

LPT  r  13  .L  13U2«K  563i»2«KJ575)*JM5)/2>J(38)*J<53)/2 

LET  L(14)»2***<(475)*2»Ki5470)*2»K(53O)  *2oKt585)*2*K(592)*2»K(588) 

LET  L(l4)-L(14>*K(508>*K(5l3)*K{520)*K(5i6f*K<620>iK{6l8) 

LPf  Lti4)sL(l4)*K(623UK<630)*J{49)+J(43)/2*Jl94)/2>K(529)*K(532> 

LPT  L(  •'  4 1“L  ( 14  >*2*K  <  4?4 )  ■••2*<  (  459>i2*K{466)i2*K(462)*K(634]l*i<(639) 

LET  L<14Ul(14>*K{646>«K(6  42)*2'‘K<566)42«<<564)*2«K(569I*2*<(576) 

IEt  L(l4)aL{l4UJ(43)/'2*J{39)'t-J(54)/2*J{47) 

LET  L(l5J‘=2*‘K<^83)*2eK^586}iK{511> *K(51<>*K<616J «>K(621)^Kt628} 

LET  L(l5)aL{15)*K(624)4J{52)/2  .  » 

LET  L(  1 5 )®L ( l5  )‘*2**^  (45?  >*24<:t  46 0 )  *K  ( 637  )*K<64 0 )*2*K  <562? 

LET  L(15)«L  {  (567  )*2«K  ( 374  )  <-2®K  (570  )  ( 37>i'J(52)/2'»-y  (  48)  . .  . 

LET  Ul6i=K(6l9)iK(622)«2»K(565)  +  2»Kt56e)*j(46) 

SU3END  ,  ■  . -  . 

SUB  nDIFFC026:K(),J()»M()»H 

LET  M#3>«-K(261)~K(264)-2*K(44)‘*K{42)-2«l<{47)-K(54)"J{6)/2-K|296) 

LET  mI3)iH{3)-’K(294!/2-!<{299)-K(306>/2-J(27W4-K  77J«K(89)«J(8>/2 

I  PT  Mf  -iisM/  (  3i5)/2-K(  3i9  )-'  J(  27  )/2''2’‘K{  2.42)“*^  ( 24o)"2*'<  {245) 

let  H(3)=M{3>-K{96)-K(96)/2-K(101 j-K(108)/2«J(85/2-K(252)«J(20> 

ti;  . 

LET  H(4)sM{4)-K(90>/2"U<9>/4-4{243)-K(246)-J(21>/2-K{99l/2 
let  H(4)4M(4)'’K(l02)/2“J{9)/4-2eK(26)-K(24)-'2«K{29)-t<l36)-J(2) 

LET  M(5)?’-K{8i)/2nK(3»|/24K(27)»K(30)”J(3)/2 
LET  ul  6) s-K  { 258 )  "K(  396 ) ”  J<  24  ) {  312 ) /2*K(  324 ) 

LET  M#7)s-2«l((256)"K(261)"’2®^(250?-*<l264)-J(22)''K(42)-K454)-J(6)/2 
LET  S!7LH{7)-K?294?/2"Ki3065;2«JC27)/4-K(310)-K(3l5)/2-Kt322) 

LET  H(7)fM(7)-K(3ia)/2«J(27)/4-X(240)*K(252>''*<<®6)/2”K(i0a)/2 
LET  H  8)*-2»K<259r2<»K5262)-.2«.K(40>-K(45)-2*K(52)-K(48,-J{4) 

LET  H(8)3M(8)-KC292)-Kc297)/2-K(304)-K(300)/2-J(25)/2-J<78)/2 

LET  M{8)5M(8)’'K(243j<»2eX(290)''X(246)’'J(2l)/2»K(94)*'K(99)/2“K(10  ) 

LET  H(8)?M<8)*K|102)/2"J(9)/4"'X(24)-'K(36)  ^.00. 

let  M(9js-2<U({43)-2eK  (46)*K(295)-<(298)-K(76>-K(81)/2*'I<<8w) 

LET  M|9)aH{9)-K(  84)/2-U{7)/2r2»K( 241)-2«K( 244 )- J( 19 ) •K< 97 )-<(  lOO ) 
LET  M  llaM  9!-^  7>/2-2.K(22>-<{27)^2«K(34)-K(30)«J<3)/2: 
let  H|lC)  =  -K(79l-K(82)t72*KC25)-2*<(28>-J{l7 

let  M(12!=iAK!98^)*|lK'562!*2lK?537)42!K«594)*J(5l>^2»K{365) 

LET  M«i2)=M{12)*2»K(377)«U{32)*K(6l7|*K(629)4j(53)/24X(638)*K(636> 
let  M!l2)=M(12>*K{641)«t<{648)*J<53)*2*Kt563>-»2«K(575)*l<<419|iK<43l) 
LET  |1{13 Jc2*X(5eo)*2*X;585)«2»K(592)*2*K(588)*vJ(49)*2*K<36B) 

I  PT  Mti-?  =H(i3)*2«K(366)*2T<(S7i>^2«X<378)4J{33»*K(620)*K{6i8) 

LET  H(i3)aM(i3)*t<(623)*K(630)*J(54)/2*K(401)tK{4l3)*J(35)/2*K(634) 
LET  M(l3)=M(i3)*X(639)6K(646)*K(642>*J(54)/2’*2*K(566)'%2*K{564) 

LEt  M(13)=M(13)*2»K(569 J-»2*<{576)iJ{47)*K(422)*X(420)*K<425)*K(432) 
LET  M(l3)=M(l3)*J(35)/2*2»K(347)*2<»Xt359j  OaKfiTo) 

LET  Ml  14  » -2 «K< 563  596 ) *2sK  (364  )<^2*K(  369  )*2*K( 376  )*2»K(  372) 

let  M(l4>sM(l4)42*K(p62)*2«<( 567)*2»K(574)i2*K(970)*J(48)4K(4l8) 

LET  H( l4)=M(i4)*K(423)4K(430)*K<426)4j{36)/2*2»X(350)*2*Kf343) 


5097  LET  Ha4)it|i(14)*2»Kj353)  +  2^«‘{3a0)*d{29! 

5098  LET  H(l9>=2«K(367H2aKU70UK<6l9)4K(622)*K(4OO)*K(4O5}*KUl2) 

5099  UET  Hii5}.-:K«i‘)UK?4O8)<j.J{34)/2*2*<(565>*2*K(560)*d(46)*K<42l) 

5i0C  LET  MCi5>=M{15}*K(424»<>J?34)/2*2»<(346>t2»K{35l>*2*K(35a) 

5101  let  M{l5>=M(i5)*2«K(354)*JJ30) 

5102  LET  H?l6?«iK(403)*KU06|*2OK(349)42»K?352)*J<2B} 

5103  SgBEND 

5104  SUB«0lFF-CC3n!K{i.J«i»N<J,H 

5105  LEV  N{U‘--2«!’(*K<i49}-2»H»K(i6iJi*H»d(l4)*2«K{797)^-2«K{0o9)|j{68) 

5106  LET  NUJ*t?J{i)^HftK<2Q3)-H*K(215)^K?S5iHK{863) 

8107  LEV  N(25s"2«H‘'K<15iJl”H«K(150)-2«H*iUliJ5>-H'»1t(162')»n(^JU5)/2 

5108  LEf  NJ2}sN(2)+2aK(8oO}*K(79B)»2*K(Oo3)-».K(BiO)*J<69)/2-4«H9K(257) 

5109  LET  N{2)?Nt2)~4»H«K(269)-2*4oJ«  23)*4*'Kf905>*4‘»K(9l7)+2*J(77> 

913.9  LET  N(2)»N(2)»H(»K«lQ5)-K»Kfl97)«H».U4J')/2^.K(a33L)*JtCXi-45)*J{71)/2 
gill  LET  N<2)aK42)'»HoK{2d6>-H<»K{204)/2-H**K(299)^H*K(2l6>/2-H»Jfl7)/2 
§112  LET  N(2)sN(2)+'<(854)*KtG52)/2*K(857-)*K-fe-6-4-f/^*Jf7l-}-/^2«H»K(131J 
§113  LET  NS2)=N(2)“2»H«K(i43>f2*<(779)+2*K<79l)^2*H«K(311)"2«H?K{323> 

HU  LET  N{2>t»N|2>-^2«K(959)^2»Ki971) 

5ll9  LET  NJ3J?-H»K<l53)»H»Kfi56)*K(801}*K{O04)-4'»H»K{2&0)''Z‘‘H*Kl25e) 

5116  let  N{3>8N(3)-.4»Kt263)n2«H*<(270)<'H»J{24)+'4*l<(908)+2*K<906) 

5117  LeT  N{3)8WS3>  +  4»K{9ti)*2*K(9i3)*J('76}-H*K<lQ8>-H«K(l86)/2-H«K(19l j 

5118  let  Kl«3>»iv!<3>"H»K(198)/2-H«J(18)/4<rK{036)*K<334)/a^K{839);iiK(646>/2 

5119  LET  N«3)?NS3)*J{72)/4»2*H«KM1)«2*H*K(53)^2«K(689t<>2*K<701)*J(59) 

Sl2Q  LET  N<3)aN(3)-2«H*K{293)-2*H*<(3o5)-H«»J{287A2»K<94l)*2»K<953)*J(8o> 

5121  LET  N(3)=N('3)'*Hi»K<207)/2-H«/<(2lO)/2-H»J(lB>/4*K(a55)/2*K(e5B»/2 

5122  LET  N{3}sN(3)*J(72)/4-2*H»K(134)-H«K(13'2)-2«H«K{137)-H*K(144)-H»J(11) 

5123  LET  N«3>«N(3)f2«K(762)<,K(780)*2*K(785)>K(792>*J(65)-2«H«K7314) 

6124  LET  N{3)«N/3)-HoK#3l2>-2«Ho‘<(3l7)^t4«K(324)-H«d(26)<'2»«/962)*K(960) 

5125  LET  Ne3)sNl3)*2*K<965)^;l<{972)TJ(80>-4»K{23934H-.4#M«Kt251)44«K(887i 

5126  LET  N(3>*N(3>+4«K(899)-H*Kf95)-H*<(i07>*‘'(743)*K<755) 

5127  LET  N(4>a-2«HttK{26lj-2»H9KC264»*2*K<909)*2,,K(9l2)«H*K(l89j/2 

5128  LET  fa(4)«N<4)-H«K(192)/2+K(337)/2*K(840)y2-2«H*K(44)-H«K<«2J 

5129  LET  N(4)»N<4)-2»H*K/47)-H*K(54)«H*d{6)/2j^24K(692>*K(690J*2««{695| 

5130  LET  M{4),Ni4)*K{7o25^g460>/2r2*H»<{296)-H»K{294)-2«H*K<299)-H«K(306) 

5131  LET  M(4>»N(4 J-H«J(27)/2*2*K<944)+<<942)*2*»K<947)*K{954)6jt8i)/2 
6132  let  N{4>4N<4)_H*K{77),.H*K<89).H*J(6)/24K(725)«.K(737)*J(62i/2 

5133  LET  N{4>sN(4>-HftK(135)-H*Kil3B)»H»jU2)/2^K{T83>*K(786) 

5134  LET  NMJ*N(4)+J<66)/2“H*^K(3l5)-W*'<{318)-H*J(27)/2*K(963f*K(966)  . 

'5135  LET  N<4>t:N{4)*J(8i>/2''4*H'»K{242>-2*H«K(240l"4*H*K(245>-2«H»<(252l 

5136  LET  N14>4N(4)-2»H*J(20)*4*K(890)>2»K(B88)+4#K«893)*2«K(900> 

5137  let  Nj4>«N(4)+2flJ(74)-H*ia93)-Ht»K(96)t/2-'H«K(lOl)-'H*K{108)/2 
9138  LET  N<4)»N{4)-HftK{e)/2*K(746)*K(744)/2*K{749)*K(756)/2*J(82|/2 

5139  LET  N(4)aN(4)-2«H»K(23)-2»H*K(35)*2*K(6.7j,)12*K(683) 

5140  let  N<5}a-H*K(45)-H«K<48)*K<693)+<<696J-HttK(297)">H«K<300)4K(945) 

5141  LET  N{5»»N{5)*K(948)-HoK(80)-4«<(78)/2-H«<<83)’*H*K(90>/2-««J(9)/4 

5142  let  N<5»aN{5»*rt(728)*K|726)/2*K(731)*K(73B)/2*J{63)/4-2*H4X(243j 

5143  let  tU5)aH{5)-2*H»K('246>-Haj<2l)4^2«Kt89l)^2«K<894>  +  J(75HHeK(99»/2 
6144  LET  N(5»'N(5)-HttK(ld2>/2-H4J<9)/4*K(747>/24K(750)/2^J(63)/4 

5145  LET  N(5)sN(5)-2*H*K{26)-HaK(24)*2ftK*K(29)-HaK<36>»H*d<2)*2»K<674| 

5146  let  N{5)aN{5)-»K(672)»?«K{677)*K(634)*J{£6) 

6147  LET  N{6)a,H«K(8i)/2-HcK(84J/2*X(729)/2+K(732)/2-H»K{27>-Hi?K{30) 

5146  let  N(6jsN(6)-Hoj{3)/2*K<675)*K(678)Tj(57)/2 

5149  LET  N(7)a-H«K(150}-H*K(.i62)-H»J(l5>/2*K(798)*K{ai0)*Jl69)Jf2 

5150  LET  N{7)cN{7)'-H«K(204)/2"H«<'(2l6)/2+K(852)-»K{e64) 

6151  LET  Nt8Ja-2»H*K(l48)-H«K(153>^2*H»K(160>-H4K(i56)-»H*JU3)42oK(796) 

5152  let  Ne8)»N(8)+K(80l)*2«K<603)*K{P04)*J(67)-2*H«K<253)-2»H<'K<270) 

5153  LET  N(8)aN<8)-H«J<24)*2eK(906)*2»4{9ia)«J(78»-H*K(186)/2-«*<<l93)/2 

5154  LET  N(8),N(8)-HftJ(lB}/2«Kte34J/2*<{846)/2*J(72)/4-HaKI?.02j 


5i55 

6156 

6157 

8158 

8159 

8160 
6161 
6162 

6163 

6164 

6165 

6166 
5167 
6168 

6169 

6170 


LET  N<B>l‘N(8)»HeK<287)/2»Hf<Ua4NH®K(210>»HaJ(le}/4'^K(e58)*K(855)/2 
LET  N(8 ) sN{ 8  j  862  « 053  S /26j<  72  (132  )-"«««( 144  )+K  (780) 

LET  N(8>sN(B)«K{792S-H9K(3i2)’'H*K(324}*K(96D>*K(9725 

LET  N(9Js~2«H6K(  151  5-26H«K  (154  h2»K  C799542«K(802)"'4»HeK<256> 

LET  N(9  )sN(9  )”2»H«K(  26i  265  )-’2»H«K(  264  >*'2«H»J(22S+4*K(904) 

LET  N(9)«N(9)«2*K(9o9)4  4»KC9l6  5*2»KC^12J*2«J(76)'-H»K{184>’*HeK(ia9)/2 

let  N(9)«:N{9>-H«K(196)'-H»KCl92}/2«H®vSn.6  J/2*K(632S4K(637)72*K(844) 

LET  N(9)»N(9)+K(840j/24j|70J/2”-H«<(42)-H®K<54>-H«»J(6}/24K{690)*4(702) 
LET  N(9)BN(9)4J{60)/2'H»K(294)-H»'<i306j'"H«j(27>/2*K(942»4j(ai)/2 

LET  N(9)8N(9)’H«K{2o5  )-H»K{208  J(16  )z24H{853  5  4K ( 856  ?  <- J  J'/O  > /2 

let  N{9)'*NJ9)-2»HoK{l3o}~HoK{135)'"2*»H«K(l42)"H«'K<l3S>"H®J#l2J/2 
Let  N(9)aN(9)42«K(776)aK(7a3)*2»K(790)*K{786Js..J(665/2<=»2<»H^K{3l0J 

LET  N(9)bN(9)-H*K(3i5)-2«H»<{  322)”H«K(3l8)'fH->«J^27  J/2<-2el<$959  J 

LET  N(9)«N£9)«.K(963  j-!>2»K(97C)^K(956)iJ(8i)/'2'’2«H«K(240  } 

LET  N(9)eN{9)-.2«H»K(252)'*-294'{0BS}4goKC9eSK;H«i<S96|/2'>H»KUOB§/2 

LET  N(9)sN(9)4K{744j/2*KC756)/2  .  : 

SUBENO 

SUB  >IDIrFC04»tjK(),Jc)5MCnh5  ^ 

LEt  N(lO)a‘’4*W*KC259i'’4»H»K(252J'»4®K(907)44eK{9iOJ“’H»K(S:87) 

LET  N/18)=N{iO  J-'H»K(l9o  J+K(635j  «-KjB38j-2»H«K{40)'^H«K(45l"a«-l4»K(52) 

LET  N(i0)aN(l0)“HoK<48}~H»J(4)*2<^<(6835«'KC6935*2®KC700S4Kf696} 

LET  N(iO)=N(iO)*J(58  5”2*^''®'*^<2^2>’'^®’^^2^7)~2«H®l<C364J'»H«K{300) 

.6T  NllO)aN{10»'>H»J(25J«'2®K{940Hi<(945|*2eKf  952)*K(94e}<>J(79) 

,ET  w(10ioN<10)’’H«K(73i/2*H5>K{90)/2''H®J(9j/4*K(725>/2*K(738}/2 
,ET  N|10)=N(10)*J(635  /4”2®H*K(l33|  ''2aH«K{136)  "H»J{  10 ) t»i2^K (7015 
,ET  N(  10  )  ~N  (10  )4  2»K  (  784  J*  J£  64 )  ™2si4«’K  ( 3l3  )  "2®H»K  ( 3l6 )  "H«  J  C  25  ) 

-ET  N(i0)sN(i0)«2®K(9  6i>+2*<£964)'«.J(79  7*-4«H»K{236l-2»>^**<C243) 

VFT  N(10)=N(iO)"46HeK(250>'r26H^K(246>"HeJ<21)*4nK(686> 

.ET  N(10)=w(10>*2eK(e9i)+4^<<59S5-!>2»K(  894).<fJ|75>"H»K{94T 
.ET  N/loiaN(3,0)-HoKi  99}/2*H6K{3D6)®H®K^102>V2-H*J(9)/4->K(742) 

LET  N(105  =  N(i0)'iK(74?}/2*Ke754  5«.K(75aj/2i..J|62{/4-.H«K(24l 
,ET  N(i0)=N(10)-H»K(36  ^^K{672^•'■KC654y  .  „  « 

.ET  N(ll)=’-2®H6K(43}*fS4HeKC46  5*2«4(69lH2oK(694>-»2*6H*Kc295) 

,ET  Nlil>t:^J(ll)-2»H«K(258H^»<(943S■^2^»K(946)■=-H«K(76>-H0K(8l)/2 
.ET  N(ll)=Nfll)''K®K/6Si'“-H«»K(64S)/2«H»Jc7J /2'6K(724)!S.K(729J/2*t<(736) 

.ET  N(ll>aN(ll)*KC  732i/24-JC6lS/2“4»H«Kee4lT'*49H«K<244»-2«H*J(l9) 

,ET  N|.«)aN{li>*4«K{889S->'4«4'{S92J*2»J(V3>"HfK(97}-H*K(iO0) 

,ET  tJ(ll)=N(ll)"H«tJ(7S/2->K^745^'3'K(7485.e.J(61>Z2-2eH»K(22  5’-MeK(27} 

,ET  N(ll)=f^{ll>-2®H45KS34J*'H6K(30)”H«4S3J/2+2»K(67o)‘^K(675} 

.ET  N(llMN(ll)*2*K|  682S<-K(67SJ«’J{57J/2 

,ET  N(12)K^H*K{79J'^H®t<{82H'4^(727Hf^(730l'-24H®K(25H2®H»K(28  J-H»J{1 ) 
,ET  N(i2)=N{l2)*2**^<6^3J*2'?!<{676UJ£55>  , 

^  E  T  N  f  1 3  J  *5  2  *  f  ^  ^  5  2  »  H  45*  K  f  4  B  5 1  H  «i>  J  ■C4iKT;2^Kjll2i}w2«^K?li33^cjJ(^5) 

ET  N(l3)=N(13)<-K«Kf52IM'H«<(539,)-:-KCll755'"K<ll®^^ 

IeT  N(i4i=2*‘U®*^(47^i'2‘'*MeKC4745<'2®H®K(479)?2®Hfi-K(4C6)*H®Jf42J 
LET  N ( 14 )  =N' ( 14 ) -2«K { 1124  ii»36K(  1122  S-ZsK?  1127  )-'2®K«  1134} 

LET  N'{i4)=.N{l4J'’J(96)*4»H*®KC53i)«'4«H»K(593)®2*H«J(50)'«4*K(i22^) 

LET'  N(i'.')sN(14)s'4e!£(l241}?;2i»-J(104}*H«K(5tl9}*H«K  (521)'t'H»U(44)/2 
LET  N(l4)=fj(14)-K(li57}-K(ll69}®J(98JZ2o2»H»K(455}*2«H«K(467) 

LET  N(j,4is.N{14)-2®K(HQ3>»’c«’f<(1113?*2®H«K(6355 ♦2oH<»K(647}-’2oK(12S3} 
LET  f^{i4jgN(l4)»2oK(l295  5«-HoK(‘i3DS*H*^K(528T*H»K(533}«H*K(b4  0) 

LET  N(i4)sN(i4)*H»U{44}/2^KCij78)^K(li76)"K(  tl0l)-K<ll80)-'J(98)/2 
LET  N(l5)=2®HsK(  472}<'26H«Kf47  7  }»2eHBK(  404  )42»H«K(48O)®J<4O} 

let  N(15  )  =N  ( 15  >- 2»K{  112 0  )’’2®K{  1125  5"2eK{  1132  )"’2»K(  1120  )"jf  94) 

LET.  ^J{15)eN(15)*4'»^^»K(5S4}^4®4ftf<(582>*4®HoK(5B7>«4»H»K(594) 

LET  N(l5)aN{15)*2oH*J(51J'»46K(12325“4«?K(1230)"4fi|<(1235>-4«K(l242) 

LET  N(i5)sN(l5)'»2«U(i05S<»He<(5i2)*HoK(5i0)4K6K(5i5)*H«K<522> 

LET  N(l5j=N{15)*H»J(45)/2’;K(1160}"K(ll58)-K(1163)-K(1170)-J<99)/2 
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I22S 

ii224 

§2?JJ 

§226 

J‘227 

S22B 

§229 

i23G 

S23i 

^232 

|233 

5234 

5231; 
5236 
gg3"/ 
.S23B 
523$ 
§24  0 


L  e  ?  m  i  s  i «  N  n.  5 }  6  2  »  H  ff.  K IJ  S  6  s  H  2  ®  M  K  U  ’/  >’  H  o  J  ( 3  2  5  -  g  ®  K  S 1 0  i  S'  'r“  2  «  K  U  0  2  5 } 

Ce?  Ki?  1.5 1  e(-ff  155  asvM  ?  ’'>>2®H'-’K{  617?)  c-goHas?  ^  6g9|  aH^tM  531  *"2®!'?  i265S 

Lev  i'\l  f  1  ‘5 )  c 15  5  “'2 cK  f  i277  j  =’ J  C  10 ?  | ? 526 1 ? 531,  I^H^K  C  ^38  j  534  5 

Lev  £ Ls j  =L'*i Us) 45 ^ /2«k « n'n  1  i n'/s  1  “K c ii86 f 3.162 1 

LE'''  N  ?  1 ‘i'f  ef'ij  1 5  S.-.J?  95  H'58  I  ■f^g-sHeK  ?  456  J  4  2«H»K  H6l  I  <-gs'4<jK  ?  4SSS 

LEV  IrjSsf<ill5?e  H«'.n30f ‘•2'''K«iia5)«'2«iKU.i04  5-2sK<ii09^-2<^Kllil6J 
Let  H(’  1 5  { “4?  15  i  '"4  C  J  *2*HsK  1 638  J  ajsHoi-CC  636  a2®H<^K  ( 6415  a2®4©K  J  643  5 

IE?  4j'i5UN{i55i-HoJjg3i-»2®4f i2865«2«‘'K'^12e4i’«’2fiK{ig89|~2®l»:^S.29S» 
iEV  •'! S  j  5  Jsfv'f  15  f "  J C 10?  563 5 C575  l”4»K 1 1211  li«4*>K  f  I22SS 

Lif  i.i«'L5Js.H!l5?eH»t?.C4i9S'H'!«4l63}JoK?lC:'6“Ji»KU0?^^  .  . 

('  E ?  N  S  i’ 6 1  e ■? a H a, K  i*  4 7§y^2.4 !■•! ..•■  K  |  4  7'  6 1  ■■' g k 1 1 2 3  5  “ 2 e i<  f  i  1  6 1  ^ 4 s K  ‘i  5 e si  S  o H 

C  E  ?  N 1  i  6  J  K>  f  1 6  S  4^  4  H  c  K  {  5  8  S  j  1 4  »  4  »  K  (  5  9  2 1*  a  4  ®  H  »  K « 5  0  B }  &  2  «  H  ®  J  H  9  [, «  4  « !  t  f  1  g  2  8 1 

■( .  E  T  e  1 6 1 « ri  f  1 6  I  '■■  4  » '-■  f  1. 2  5  3  ^  «  4  »  K  1 1 2  4  0 1 "  4  a  K  e  1 2  3  6  ^  -  2  «  d  e  1 0  3 1  ■:■  H  «  K  ^  5  0  a  S 

L  E  'f  V'H  1  6  1  •■  (■'  1 1 6  )  *  l-i  o  K  i  5 1 3 1  <■  'i «  4  :<!  5  2  fi  5  *  .H ;;  K  f;  5  i  6  >  4'  H  a  JC  4  5  I  /  ?: «  K  ( 1 1 5  6  .W 14  J  i  1 6  i  ? 

IE?  N  i «  6 1  c  i\! !  1  6  ?  -  K  “  1 4  6  8  J  -  K  f  J  -s  6  4  I  u  -.If  7  ?  2  *  2  ®  H  !<  S  3  6  a  J  2  »  H  a  K  i  3  6  6 1 

t  E  f  M  C 1  a  « 1 6  H  2  «■  H  ti  K  ?  3  7  :U  4  2  »  4  «.  ri  C  3  7  8  ?  *  H  e  d  { 3  S  K.  2  »  K  1 1 3 1 6  J  -  2  e  K  n  0  i  4  J 

I.EY  N  S 16  S  ( :\6 }  -'ac'li  J  i(il9  J  «’.'2«'KU026 1-4?  S?  S  ‘>2®K«^‘''  ?  02G }  fiS  I 

Ik?  N«i6J=NC16)a2«Ht>KU>23f-v2®^»'<?630?<HeJ(54  5--2®Ka268)-"2«l<U266l  , 

LG'f  N|*l6)aKs^i6 Jp2eK;l,27i J'=2siC{l279J'-J|103HH®'K5  40!l)'«'H«K^4^S1 

LET  N  C 16 1  sN  i  16 } (35 1  <'2“  K  U.D 49 1  “K  ( 106 1 J  >» J  C  B9  J  /  2aH®K  ?  529 1  c*!4otC  ?  S32^ 
LET  Li  C  4  6  >  sfH  ?  6 }  ■S’H'SvH  43 1  /  g’”!*  C  3.<  77  J  “*K  { i3.0p  1  '*  J ?  97  5  /  gAg-oH^sK  C  4§4 1 
^  T  ^  ^  6 )  —  M  ( i  6 1  ^  2  i*  H  {.^  K  C  4  5  9  p  1'  2  o  H  <v  It  C  4  6  6 }  a  2  45  H  ft  iC  ^  4  6  g )  a  11 J  ^  3  9 1  ™  2  ^  K  ^  1 1 0  ^  I 

Le?  NSl6)s?'.iU6>''2*KCli07^-  2»!<n.il4  J~2OK(lii0  ^ 

lEf  NC16  isW^C-si  c-2c>Hc;K(  6391  aSO'tasf?  546  J  a2®H»KC  642)  <>N-»dC,54 
Let  N.l6l"iM;16?-2®K;i237  5-2<5K^1290J-'dcl08H;4«H;jK;566|-¥4‘‘H®iq564J 
(jS?  M?{6isNC16)-i>4t-H6K?'5691'V4®H®i<C576f*2»H»PC47}'^46KCl(?l4J . 

LET  fl?  tfiJsNf  1224}  “’2o  J?  16£5*H®K(  gSgf^H^K  C  4i!?0S  *H»K  C  425> 

{.ET  N«  i6  j-sN(l6Si.H'»K  I  432)4HyJ4  35}/2-"KC1070)"K(106H}>*K  tiD73|«tj|j89};'g 
lET  N$i6}sNCl6}«2»HeK(3‘v7}'42®Ha!v?3595“2»K{995}'’2»K?1007}  . 

t,E?  Nf i7.)n4«H«K{583|-e-4cH<>l'C 5855e4sK«12315-4®K(l234?<'H®Ki{511}  . 

UE?  W(i7}aN?s7},,,H«K(5i45“KCll595''<Ul62}*2«H»K{364>,;.2«-H»!U369s 
LET  LjS17}KMf l75a2eHais:(376^'1^2»‘^)»ir.C372HHo^K31>“3»Ka0i2}-2sitU017} 

LfT  r'^?l7}s(M(l7}-26K{102‘1}«g<5!C  (1-020} -J? 05}  *9»H«K« 63.6)  621) 

IsE?  N|  i7)bN  j  l7)e2®HaKS  6g4)*H®j!:52)'’2^’K§1264}“2®KU269) 

let  NS^7)BM{l7)-2sKci276}*2R'KCi272)-J?ie6)^H®K{4C4)<s.H®K|4eaS  , 

UET  NC 17 }  eN«  1?  5  4-H»K  f  407  )  414 )  *HavH36 } /2«K«  1052  )~KU05fl)“>K  1105^1  . 

lET  N(  17  )“W?  17  I'-K  ?  18(^2  s»,J(?0 } /2'!-2»H»K  (457  C  460  Jig?  ) 

Li?  Ma7j3N(l7)'-2oit(ll05)'^2eK<li0S)*'JC9l}-!-2®H«KC637>*2aH»tt«e40) 

Lev  N(i7ssN{i7l*"*-!^2)'=*2®KU2S5H'2®l'?123B)“>.)(106)o4®H&K556p 

Let  N  f  1  7 }  sN  j  i7  }  4 4C’H®K  ( 567 )  4  ( §741  a4®H®K  (570 ) J  ( 4 S  } 

{,ET  WC  17)  sNIiT}^ 4«Kf  i2|;fj}<»4»K(l21.5l-4«K{  i222}'-'4«Kll2l8}”2»JU02} 
kEV  N(57)-'M(37S'!'K'eKC4i8)<>K^iL{4?S}*!^®^?4S0)-SH«'*^?4ge)}»Ha,J<3S};^2 
I ET  N  ( 1 7 ) aW  f  17  J  -K  ( i D 6S  }  ■■[(  ( 10 7.1 )  «(<  ( 1078 )  »K  ( i 074  ) «J  ( 90  ) »K  « 350 ) 
IE?  Mei7)KM(i7)-t2»HDK(  348}->2«4«!t(333?'^2«HoKC360).&HoJ(29  5-»gsKf 996) 

Let  Li7)sN/l7)*'>2®*<|996)'«'2?'t'(100i)"’29!<(10085«d{B3| 

(oE?  M  ( 18  5  ygsHoK  j  369  }-!.2<sU»it  §370  }  t)2»K  { 1015  5  “2®K  1 1010  f  319 1 

tET  L!Ci8)3K«in}-:*2®HeKC622)‘;2«-^«  l267>‘'2eK{l270J'^^l**<?4BG}-:'H«in4O5) 

LET  Ma6  5£;NCl6.-*H«K(4i^)■^.H^^4.?4(|B)*H^^J?34  3/2“K£lO48)•“tt^lO53)-K(lO60) 

LET  Lj( ie)a;\|(18}«KClo56}'=J(88  5/2*4®M«H(565)44®H«Ki(568)'J’2®H®Jf 46)  . 

let  NM'8}aM(i0J«4aK  (1213 }®4^K(  1216  5-20  J  {100  )'-'>W^‘K(i523.)‘i>H<)K§4§4) 

LET  l\iiia}aN{l8)4K&U?34}/2.'K{1069)'»K{l072)-J(3B)/S 

J  E?  N  { 5  0 )  sN  { ■!  8  )  -ig'sHaK  {  346  )  ■52®'°)®'^  ?  S^l  ^  3501*2®^’®®'?  SS'H’S'HsJtgfi ) 

U;;T  L!{iO)«M{i0)“2c.K(994)»'2t?i<'C99  9)-f2«K{lO06f’"2{iK{10C2)-J<84} 

LET'  M{l95sH®K(403  (40  6  )  •’K  { iOSi  )“K  { 1954  §  349  )42eM»U  {352 ) 

LET  Nci9SaN{i9?t5l-!oJf28)-2®Kc997)“2fflK{18SS}‘-J{62) _ _ _ _ 

SU^IEMD  . 


MATCH 

RiO  FOR  P«8  TO  1 

20  FOR  Q«0  TO  1.  . 

30  FOR  U1  TO  3 

40  FOR  JeX  TO  3 

50  FOR  Kai  TO  2 

i^O  FOR  laX  TO  2 

70  FOR  Mai  TO  3 

60  FOR  Nal  TO  3 

. '  90  LET  TaT*l 

^  100  PRINT  I  jJjK«‘4lL’f'4?M<NlPlQjT4B(SS)lT 

. 110  NEXT  N 

120  NEXT  M 
130  NEXT  U 

_  140  NEXT  K  _  . . . . . 

150  next  J 
160  NeXT  I 
170  NeXT  Q 
180  NEXT  P 
190  print 

__20O  LET  T?0  . . . . . . . . . .  ... 

210  FOR  PsO  TO  1 

220  FOR  Q?0  TO  i  . 

230  FOR  lal  TO  3 
_  _240  FOR  Jgl  TO  3 
250  FOR  Kal  TO  3 
___260  let  TsT*!  .  ' 

.270  PRINT  I; J;K+3JPJ0iTAB(35) S T 
'  _  2®0  ^ 

:290  NEXT  4 

300  NEXT  I  . . . . 

310  NEXT  Q 

_  320  NEXT'  P- _ _ _ _ _ _ 

330  END 
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MATCH 


115  5  110  0 

1155120  0 

115  5  13  0  0 

11552100 
11552200 
1155230C 
115  5  3  10  0 

115,  §3200 
11,5  53300 

115  6  110  0 

11561200 
11561300 
Il562l00 
11562200 
11562300 
11563100 
11563200 
11563SOO 

116  5  110  0 

1  1  6  9  1  2  0  0 

11651300 
11652100 
11652200 

1165  2  300 

11693100 

I  1  6  5  3  2  0  0 

11653  3  00 

11661100 

II  '6 '6 1  200 

11661300 
1  1  6  6  2  1  0  0 

1  1  6  6  2  2  0  0 

11662  3  00 

1166  3100 

1  1  6  6  3  2  0  b 

11663300 
1  2  5  5  1  1  0  0 

12551200 
1  2  5  5  1  3  0  0 

12592100 
I  2  5  5  2  2  0  0 

1255  2300 

12553100 
12593200 
1  2  5  5  3  3  0  0 

I  .2  5  6,  1  1  0  0 

1  2  561200 

1  2  5  6  1  3  0  0 

12562100 
1  2  5  6  2  2  0  0 

12562300 
1  2  5  6  3  1  0  0 

1256320  0 

125  6  3300 
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13652111 
i  3  6  S  2  2  1  1 
13652311 
156531  1  1 
13  6  5  3  2  11 
13653311 
1  3  6  6  1  1  1,1 
13  6  6  12  11 

1  3  6  6  1  3  1  1 
1  3  6  6  2  1  1  1 
1  3  6  6  2  2  1  1 
13  6  6  2  3  11 
1^36,6.3  1,11 
13  6  6  9  2  11 
13663311 

2  1  5  5  1  ill 
2  15  9  12  11 
2  1  5  5  1  S  1  1 
2  %  5  5  21  1.1 
2  15  5  2  2  1  1 
2  15  5  2  3  1  1 
2  15  5  3  111 
2155  3  211 
21553311 
2  1  9  6  1.  1  1  ,  1 
2  15  6  12  11 
21961311 
2  1  5  6  2.  1  1  1 
2  1  5  6  2  2  1  1 
215623  11 
21563111 
2  19  6  3  2  1  1 
2  1563311 


10  41 

1042 

1043 

1044 

1045 

1046 

1047 
104B 

1049 

1050 

1051 

1052 

1053 

1054 

1055 

1056 

1057 
1059 
1059 
1050 
1061 
1062 

1063 

1064 

1065 

1066 

1067 

1068 

1069 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 
1076 
1079 

1090 

1091 

1092 

1093 

1094 
1085 
1096 
1087 

J088 

1099 

1090 

1091 

1092 

1093 

1094 

1095 

1096 

1097 

1098 
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2  1  6  5  1  i  %  1 

216512  11 

2  1  6  5  13  1  1 

21662111 
21652211 
21652311 
21653111 
216532  1  1 

2  1653311 

2166  1111 
2166  1211 
216613  11 

2166  2111 
2  1  6  6  2  2  1  1 

2  1  6  6  2  3  1  1 

21663111 
2  1  6  6  S  2  1  i 

2166331  1 

2  8  9  5  1  1  1  1 

22  5  51211 

2  2  5  5  1  3  1  1 

2  2  5  5  2  111 

2  2  5  5  2  2  1  1 

2  2  5  5  2.  3  1  1 

2  2  5  5  3  1  1  1 

2  2  5  5  3  2  1  1 

2  2  5  5  3  3  1  1 

2  2  5  6  1  1  1  1 

2  2  5  6  12  15. 

2,  2  5  6  1  3  1  1 

2  2  5  6  2  1  1  1 

2256  2  211 

2  2  5  6  2  3  1  1 

2256311  1 

2  2  5  6  3  2  1  1 

2  2  5  6  3  3  1  1 

2  2  6  5  1  1  1  1 

2  2  6  5  1  2  1  1 

2  2  6  5  13  11 

2  2  6  5  2  1  1  1 

2  2  6  5  2  2  1  1 

2  2  6  5  2  3  1  1 

22653  111 

2265  3  211 

2  2  6  5  3  3  1  1 

2266  1111 

2  2  6  6  1  2  1  1 

2  2  6  6  1  3  1  1 

2  2  6  6  2  1  1  1 

2  2  6  6  2  2  1  1 

2266231  1 

2  2  6  6  3  1  1  1 

2  2  6  6  3  2  1  1 

2  2663311 

2  3  5  5  1  1  1  1 

2  3  5  6  1  2  1  1 

2  3  5  5  1  3  1  1 

2  3  5  5  2  1  1  1 


1Q99 
1100 
1 1  ti  J 
1102, 
1103 
110  4 
1105 
110  6 

1107 

1108 

1109 

1110 
111! 
1112 

1113 

1114 

1115 

1116 
1117 
lllB 

1119 

1120 
112.1:' 
1122. 

1123 

1124 

1125 
1125^ 

1127 

1128 

1129 

1130 

1131 

1132 

1133 

1134 

1135 

1136 

1137 

1138 

1139 

1140 

1141 

1142 

1143 

1144 

1145 

1146 
114  7 
1143 

1149 

1150 

1151 

1152 

1153 

1154 

1155 

1156 


1  4n 


n  ?  3  5  5  2  2  1  i  1157 

r  I  2  3  5  5  2  3  1  1  1158 

23553111  1159 

2  3  5  5  3  2  1  1  1160 

2  3  5  5  3  3  1  1  1161 

2356  1111  1162: 

2  3  5  6  1  211  1163 

23561311  1164 

23562111  1165 

2  3  5  6  2  2  1  1  1166 

2  3  5  6  2  3  1  1  1167 

23563111  1168 

2  3  5  6  3  2  1  1  1169 

23-563311  1170 

23651111  1171 

2  3  6  6  1  2  1  1  1172. 

2  3  6  5  1  3  1  1  1173 

23652. Ill  1174 

2  3  652211  1175 

2  3  6  5  2  3  1  1  1176 

2  3  6  5  3  1  1  1  117  7 

23653211  1178 

23653311  1179 

2  3  6  6  1111  1180 

2  3  6  6  1  2  1  1  llBl 

23661311  11B2 

2  3  6  6  2  1  1  1  1183 

'  2  3  6  6  2  2  1  1  119^ 

23  662311  11S5 

2  3  6  6  3  1  1  1  1186 

2  3  6  6  3  2  1  1  1199 

23663311  1188 

31551111  1189 

3  15  6  12  11  1190 

3  15  5  13  11  1191 

3  1  5  5  2  1  1  1  1192 

31562211  1193 

3  1  5  5  2  3  1  1  1194 

3  1  5  5  3  1  1  i  1195 

31553211  1196 

31953311  1197 

3  15  6  1111  1198 

3  1  5  6  12  1  1  1199 

3  15  6  13  11  1200 

3  1  5  6  2  1  1  1  1201 

31S62211  1202 

31562311  1203 

31563111  1204 

3  1  9  6  3  2  1  1  1205 

315633  11  1206 

31651  111  1207 

3  1  6  5  1  2  1  1  1208 

316513  11  1209 

31652  111  1210 

3  1  6  5  2  2  1  1  1211 

31652311  _  1212 

31653  1  11  1213 

3  .1  6  5  3  2  1  1  1214 
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1 

6 

9 

3 

3 

i 
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1 

6 
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1 

1 

1 
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1216 
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X 

6 

6 
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o 
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X 
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X 

X  . 

1264 

3 
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X 

i 

1259 

3 

3 

9 
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X 

1 

1266 

3 
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1258 

3 

3 

9 

9 

3 

-y 

O 

X 

1 

1269 

3 

.  3 

9 

6 

1 

1 

X 

1 

1270 

3 

3 
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o 

X 

1 

1271 

3 

3 

9 

6 

1 

3 

X 

1 

1272 

n,  3  3  5  6  2  1  1  1  1273 

'I  3  3  5  6  2  2  11  1274 

3  3  5  6  2  3  11  1275 

3  3  5  6  3  1  1  1  1276 

33563  211  1277 

3  3  5  6  3  3  1  1  1278 

3  3  6  5  1  111  1279 

3  3  6  5  1  2  1  1  1230 

336  5  1311  1231 

'  3  3652111  1202 

3  3  6  5  2  2  1  1  1293 

3  3  6  5  2  3  1  1  1284 

3  3  6  5  3  111  1285 

33. 6  53211  1296 

3  3  6  5  3  3  1  1  12S7 

33661  1  11  1288 

3  3  6  6  1  a  1  1  1239 

3  3  6  6  13  11  1290 

3  3  6  6  2  1  1  1  1291 

3  3  6  6  2  2  1  1  1292 

336  6  2  311  1293 

33663111  1294 

3  3  6  6  3  2  1  1  1295 

3  3  6  6  3  3  11  1296 

1  1  4  0  0  . 1 

1  1  5  0  0  2 

1  1  6  0  0  S 

1  2  4  0  0  4 

1  2  5  0  0  5’ 

1  2  6  0  0  6 

1  3  4  0  0  .  .  7 

1  3  5  0  0  8 

1  3  6  0  0  9 

2  1  4  0  0  10 

2  1  5  0  0  11 

2  1  6  0  0  12 

2  2  4  0  0  13 

2  2  5  0  0  14 

22600  .15 

2  3  4  0  0  16 

2  3  5  0  0  17 

2  3  6  0  0  18 

3  1  4  0  0  _  19 

3  1  5  0  0  20 

3  1  6  0  0  21 

3  2  4  0  0  22 

3  2  5  0  0  23 

3  2  6  0  0  24 

3  3  4  0  0  25 

3  3  5  0  0  26 

3  3  6  0  0  27 

114  0  1  26 

115  0  1  29 

116  0  1  30 

1  2  4  0  1  31 

12  5  0  1  32 

1  2  6  0  1  33 


143 


■  1  ?  6  0  X 

2  14  0  3, 

2  1  5  0  1 

2  16  0  1 
2  2  4  0  1 

2  2  5  0  1 

2  2  6  0  1 

2  3  4  0  1 

2  3  5  0  1 

2  3  6  0  1 

3  14  0  1 

31-50  1 

3  16  0  1 

3  2  4  0  1 

3  a  5  0  1 

3  2  6  0  1 

3  3  4  0  1 

3  3  5  0  1 

3  3  6  0  1 

11  4  1  0 

1  1  5  J  0 

116  10 

12  4  10 

.  1  2  5  1  0 

12  6  10 

1  a  4  1  0 

13  5  10 

13  6  1  0 

2  14  10 

2  1  5  1  0 

2  16  10 

2  2  4  1  0 

2  2  5  1  0 

2  2  6  1  0 

2  3  4  10 

. 2  3  5 . 1  0 

2  3  6  1  0 

3  1  4  1  0 

3  15  10 

_  .3  16  10 

3  2  4  1  0 

. 3.  2  .5.1  Q 

3  2  6  1  0 

3  3  4  1  0 

3  3  5  1  0 

3  3  6  1  0 

114  11 

-.1  1 . 5  1  1 

1  16  11 
12  4  11 

12  5  11 

12  6  11 

1  3  4  11 

. . .1  .3  .5  .1  I 

1  3  6  11 

2  14  11 


34 

35 

36 
3’/ 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
40 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
33 

84 

85 

86 

87 

88 

89 

90 

91 
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^  V! 


2  16  11 

2  2  4  11 

1  1 

1  1 
1  1 

1  1 

2  3  6  1  1 

3  14  11 

3  15  11 

3  16  11 

3  2  4  11 

3  2  5  1  1 

3  2  '  6  11 

s3  3  4  1  1 

3  3  5  1  1 

3  3  6  11 


92 

93 

94 

95 

96 

97 

98 

99 

100 
101 
102 

103 

104 

105 

106 
107 
109 


R 

DIFFCO 

VALUES  OF  K(  1  ,J,K#L»M,N*P.Q)  ArE  <'iULTlPLY  BY  E6> 


115  9 

115  5 

11  5  5 

115  9 

115  5 

115  5 

115  5 

I  1  5  9 

115  5 

II  5  6 

1  1  5  6 

115  6 

1  1  5  6 

I  15  6 

II  5  6 

115  6 

115  6 

115  6 

11  6  5 

116  5 

116  5 

I  16  5 

116  5 

1  1  6  9 

116  5 

116  5 

I  1  6  5 

1  16  6 
116  6 
1  16  6 

II  6  6 
.116  6 

116  6 
116  6 
116  6 
116  6 

12  5  5 

1  2  5,5 

1  2  5  5 

1  2  5  5 

12  5  5 

12  5  5 

12  5  5 

12  5  5 

12  5  9 

1  2  5  6 

12  5  6 

1  2  5  6 

12  5  6 

I  2  5  6 

12  5  6 

12  5  6 


110  0 
12  0  0 

13  0  0 

2  10  0 

2  2  0  0 

2  3  0  0 

3  10  0 

3  2  0  0 

3  3  0  0 

1  1  0  0 

12  0  0 

13  0  0 

2  10  0 

2  2  0  0 

2  3  0  0 

3  10  0 

3  2  0  0 

3  3  0  0 

110  0 
12  0  0 

13  0  0 

2  10  0 

2  2  0  0 

2  3  0  0 

3  10  0 

3  2  0  0 

3  3  0  0 

1  1  0  0 

12  0  0 

13  0  0 

2  10  0 

2  2  0  0 

2  5  0  0 

3  10  0 

3  2  0  0 

3  3  0  0 

1  1  0  0 

12  0  0 

13  0  0 

2  10  0 

2  2  0  0 

2  3  0  0 

3  10  0 

3  2  0  0 

3  3  0  0 

110  0 
12  0  0 

13  0  0 

2  10  0 

2  2  0  0 

2  3  0  0 

3  10  0 


4.72616  E-2 
1,56731  E-2 
"1,70966  E-4 
1,12845  E"? 
6.91198  E-2 
"3.42554  E”3 
1,04722  E"? 
«1, 55702  E-'9 
4,38228  E’»3 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

4.72616  E"2 
1,66731  E"2 
■  1,70966  £’*4 
1,12845  E-2 
6,31190  E-2 
-3.42554  E-3 
1,04722  E-7 
-1,55782  E"9 
4,38220  E-3 
9,46502  e-3 
6,0569  E-3 
-3.2572  E-S  , 
2,17294  E-3 
2,44995  E"2 
-6,52624  E~4 
1,99513  E-B 
-3.15843  E-10 
0.00129 
0 
0 
0 
0 
0 
0 
0 
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12  5 

12  5 
12  6 
1  2  6 
12  6 
12  6 
12  6 
12  6 
12  6 
12  6 
12  6 
12  6 
12  6 
1  2  ‘6 
12  6 
12  6 
12  6 
12  6 
12  6 
12  6 

13  5 

15  5 

1  3  5 

13  5 

13  5 

13  5 

13  5 

13  5 

1  3  5 

13  5 

13  5 

1  3  5 

1  3  5 

1  3  5 

13  5 

13  5 

13  5 

1  3  5 

13  6 

13  6 

13  6 

1  5  6 

13  6 

1  3  6 
13  6 

13  6 

13  6 

13  6 

13  6 

1  3  6 

13  6 

13  6 

13  6 

15  6 

13  6 

13  6 

2  1  5 

2  1  5 


6  3  2  0  0 

6  3  3  0  0 

5  110  0 

5  12  0  0 

5  13  0  0 

5  2  1  0  0 

5  2  2  0  0 

5  2  3  0  0 

5  3  10  0 

5  3  2  0  0 

5  3  5  0  0 

6  110  0 

6  12  0  0 
6  13  0  0 

6  2  10  0 
6  2  2  0  0 

6  2  3  0  0 

6  3  10  0 

6  3  2  0  0 

6  3  3  0  0 

5  1  1  0  0 

5  12  0  0 

5  13  0  0 

5  2  10  0 

5  2  2  0  0 

5  2  3  0  0 

5  3  10  0 

5  3  2  0  0 

5  3  3  0  0 

6  110  0 

6  12  0  0 
6  13  0  0 

6  2  10  0 
6  2  2  0  0 
6  2  3  0  0 

6  3  10  0 

6  3  2  0  0 

6  3  3  0  0 

5  110  0 

5  1  2  0  0 

5  13  0  0 

5  2  1  0  0 

5  2  2  0  0 

5  2  3  0  0 

5  3  10  0 

5  3  2  0  0 

5  3  3  0  0 

6  110  0 

6  12  0  0 
6  13  0  0 

6  2  10  0 
6  2  2  0  0 
6  2  3  0  0 

6  3  10  0 

6  3  2  0  0 

6  3  3  0  0 

5  110  0 

5  12  0  0 


0 
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2  3  5  6  2  2 

2  3  5  6  2  3 

2  3  5  6  3  1 

"2  3  5  6  3  2 

2  3  5  6  3  3 

2  3  6  5  1  1 

2  3  6  5  1  2 

2  3  6  5  1  3 

2  3  6  5  2  1 

2  3  6  5  2  2 

2  3  6  5  2  3 

2  3  6  5  3  1 

2  3  6  5  3  2 

2  3  6  5  3  3 

2  3  6  _^6  1  1 

"2 '3  6612 
2  3  6  6  1  3 

2  3  6  6  2  1 

2  3  6  6  2  2 

2  3  6  6  2  3 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2.20371  E-S 
9,24799  E-4 
►  7,3814  £-^6 
5,21456  £•'4 
0,003765 
►1.57914  E-4 
4.32758  E-9 
►7,6424  E-11 
2,26696  E-4 
’7.67054  E-5 
•2,96745  E»5 
2,75874  E-7 
•1,32311  E-S 
•1,21001  E-4 
5.5275  E-6 
-1,68981  E-10 
2,67508  E-l2 
-7,50906  E-6 

0  . 

0 

0 

0 

0  . 

0 

0 

0 

G 

0 

0 

0 

0  . 

0 

0 

0 

0 

0 

-7,67054  E-5 
-2,96745  E-5 
2.75874  E-7 
-1,82311  E-5 
-1,21001  e-4 
5,5275  E-6 
-1,68981  E-10 
2.67508  E-12 
-7,50906  E-5 
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3  13  5  1 
3  15  5  1 
3  1  5  5  1 
3  15  5  2 
3  1  5  5  2 
3  15  5  2 
3  15  5  3 
3  15  5  3 
3  15  5  3 
3  13  6  1 
3  15  6  1 
3  15  6  1 
315  6  2 
31-562 
3  15  6  2 
3  15  6  3 
3  15  6  3 
3  15  6  3 
3  16  5  1 
3  16  5  1 
3  1  6  5  1 
3  1  6  5  2 
3  16  5  2 
3  16  5  2 
3  16  5  3 
3  1  6  5  3 
3  16  5  3 
3  16  6  1 
3  16  6  1 
3  16  6  1 
3  1  6  6  2 
3  16  6  2 
3  16  6  2 
3  1  6  6  3 
3  1  6  6  3 
3  16  6  3 
3  2  5  5  1 
3  2  5.  5  1 
3  2  5  5  1 
3  2  5  5  2 
3  2  5  5  2 
3  2  5  5  2 
3  2  5  5  3 
3  2  5  5  3 
3  2  5  5  3 
3  2  5  6  1 
3  2  5  6  1 
3  2  5  6  1 
3  2  5  6  2 
3  2  5  6  2 
3  2  5  6  2 
3  2  5  6  3 
3  2  5  6  3 
3  2  5  6  3 
3  2  6  5  1 
5.2  6._.5  1 
3  2  6  5  1 
3  2  6  5  2 


1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 
3 
1 
2 
3 
1 
2 
3 
1 
2 
3 
1 
2 
3 
1 
2 
3 

1 

2 
S 
1 
2 
3 
1 
2 
3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 

1 


1  0 

1  0 

1  0 

1  0 

X  0 

1  0 

i  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1  0 

X  0 

1  0 

1  0 

1  0 

X  0 

1  0 

1  0 

1  0 

1  0 

1  0 

1..  0 
1  0 

a,  0 

1  0 

10 
1  0 

.10 
1  0 

1  0 

1  0 

1  0 

1  0 

.  1 . 0 

1  0 

1  0 

1  0 

1  0 

1  0 

1..  0 
1  0 

1  0 


2,34495  E-9 
9,071^4  E-IO 
-8,43371  E-12 
5,57339  E-10 
3,69911  E-9 
-1.68981  E-10 
5.1658B  E-15 
-8,17797  E-17 
2.29558  e-10 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2,34495  E-9 
9,07174  E-10 
-8,43371  E-12 
5.57339  E~10 
3. 69911  E-9 
"1,68981  E-10 
5,16588  E-15 
•8.17797  E-17 
2,29558  E-iO 
-3,71223  E-li 
-1,43612  E-11 
1,33512  E-13 
-8,32308  E-12 
-5,95596  E-li 
2,67508  e-12 
-9,17797  E-17 
1,29463  E-18 
-3.63407  E-12 

...  0 . . 

0 

0 . 

0 

0 

0 

,  ,0  . : . . 

0 

0 . 

0 

. 0 . . 

0 

0 
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3  2  6  S  2 

3  2  6  5  2 

3  2  6  5  3 

3  2  6  5  3 

3  2  6  5  3 

3  2  6  6  3, 

3  2  6  6  1 

3  2  6  6  1 

3  2  6  6  2 

3  2  6  6  2 

3  2  6  6  2 

3  2  6  6  3 

3  2  6  6  3 

32-663 
3  3  5  5  1 

3  3  5  5  1 

3  3  5  5  1 

3  3  5  5  2 

3  3  5  5  2 

3  3  5  5  2 

3  3  5  5  3 

3  5  5  5  3 

3  3  5  5  3 

3  3  9  6  1 

3  3  5  6  1 

3  3  5  6  1 

3  3  5  6  2 

3  3  5  6  2 

3  3  5  6  2 

3  3  5  6  3 

3  3  5  6  3 

.  3  3  5  6  3 

3  3  6  9  1 

3  3  6  5  1 

3  3  6  5  1 

3  3  6  5  2 

3  3  6  5  2 

3  3  6  5  2 

3  3  6  5  3 

3  3  6  5  3 

3  3  6  5  3 

3  3  6  6  1 

3  3  6  6  1 

3  3  6  6  1 

~3  3  6  6  2 

3  3  6  6  2 

3  3  6  6  2 

3  3  6  6  3 

3  3  6  6  3 

_  3 . 6 . 6  _  3 

115  5  1 

115  5  1 

115  5  1 

1  15  5  2 

115  5  2 

-1-1- . 5  5  2 

115  5  3 

1  15  5  3 


2  10 
3  10 

110 
2  10 
3  10 

110 
2  10 
3  10 

110 
2  10 
3  10 

1  1  0 
2  10 
3  10 

110 
2  10 
3  1  0 

110 
2  10 
3  10 

110 
2  10 
3  10 

110 
2  10 
3  1  0 

110 
2  10 
3  10 

1  1  0 
2  10 
3  10 

110 
2  10 
3  10 

1  10 
2  10 
3  10 

110 
2  10 
3  10 

1  1  0 

2  10 

3  10 

110 
2  10 
3  10 

110 
2  10 
3  10 

111 
2  11 
3  11 

1  I  1 

2  11 

5  11 

111 
2  1  1 


0 

0 

0 

0 

0 

-3,71223  E-11 
-1,43612  E-11 
1,33512  E-13 
-8,32308  E-12 
-5,95596  E-li 
2,67508  E-12 
-9,17797  E-17 
1,29463  E-18 
-3,63407  E-12 
1,24924  E-4 

5.11829  E-5 
-4,47543  E-7 

2,96001  E-5 
2,08468  E-4 
-9,96712  E-6 
2,74133  E-10 
-4,3397?  E-12 
1,26625  E-5 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

0 . . 

0 

0 

0 

1,24924  E-4 

5.11829  E-5 
-4,47543  E-7 

2,96001  E-5 
2,08468  E-4 
-8,96712  E-6 
2.74153  E-10 
-4.33972  6-12 
1,26625  E-5 
-2,07126  E-5 
-1,32791  E-5 
6.81961  E-6 
-4,59427  E-6 
-7,36859  E-5 
1,3664  E-6 
-4,1772  E-11 
6,61282  E-13 
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Pi 


1195  3  311 

11561111 
115  6  12  11 

11561311 
11562111 
11562211 
11562311 

115  6  3111 

11963211 
1  1  563S11 

1169  1111 

11691211 

116  5  13  11 

11-652111 
11652211 
11652311 
11653111 
116  5  3  211 

1  1  6  5  3  3  1  1 

1  1  6  6  1  111 
1166  1211 
116  6  13  11 

1166  2  111 
11662211 
11662311 
11663111 
116632  11 

1  i  6  6  3  3  1  1 

12551111 
I  2  5  5  1  2  1  1 

12551311 
1  2  5.5  2  1  1  1 

12552211 
12552311 
1  2  9  5  3  1  1  1 

12553211 
12553311 
1  2  5  6  1  1  1  1 

12561211 
12361311 
1256211  1 

1  2  5  6  2  2  I  1 

12562311 
1  2  5  6  311  1 

1256  3  211 

125  65311 

126911  11 

12651211 
126313  11 

1.  2  6  9  2  1  1  1 

12652  211 

1  2  6  5  2  3  1  1 

12653111 
12653211 
12653311 
1  2  6  6  1  1  1  1 

1  2  6  6  1  2  1  1 

1  2  6  6  1  3  11 
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12  6 
1  2  6 

12  6 
12  6 
12  6 
13  5 

1  3  5 

1  3  5 

1  3  5 

13  5 

13  5 

1  3  .5 

1  3  •_  5 

13  5 

.1  3  5 

13  5 

1  3  5 

1  3  5 

13  5 

"^1  3  5 

13  5 

1  3  5 

13  5 

13  6 

.1  3  6 
13  6 

. 1  3  6 

13  6 

1  5  6 
13  6 

_.,1 . 3  ._..,6„ 

13  6 

. 13  6 

13  6 

_  1 . 5  6 

1  3  6 

_  1.  3  6 
15  6 

1  3  6 

13  6 

1  . 3  6 

2  1  5 

. 2 . 1.5. 

2  1  5 

. 2  1  5 

2  1  5 

2  I  5 
2  1  5 

.2 . 1  .  5 

2  1  5 

2  1  9 

2  1  5 

2  1  5 

2  15 

2  -  1  5 

2  1  5 

2  1  5 


6  2  11 
6  2  2  1 
6  2  3  1 

6  3  11 

6  3  2  1 

6  3  3  1 

5  111 

5  12  1 

5  ■  1  S  1 
5  2  11 

5  2  2  1 

5  2  3  1 

5  3  11 

5  3  2  1 

5  3  3  1 

6  111 

6  12  1 
6  13  1 

6  2  1  1 
6  2  2  1 
6  2  3  1 

6  3  11 

6  3  2  1 

6  3  3  1 

5  111 

5  12  1 

5  13  1 

5  2  11 

5  2  2  1 

5  2  3  1 

5  3  11 

5  3  2  1 

5  3  3  1 

6  ill 

6  12  1 
6  1  3  1 

6  2  11 
6  2  2  1 
6  2  3  1 

6  3  11 

6  3  2  1 

6  3  3  1 

5  111 

5  12  1 

5  13  1 

3  2  11 

5  2  2  1 

5  2  3  1 

5  3  11 

5  3  2  1 

5  3  3  1 

6  111 

6  12  1 
6  13  1 

6  2  11 
6  2  2  1 
6  2  3  1 

6  3  11 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

1 

1 

X. 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


»4. 32481  E-6 
-3,12943  B-5 
1,30965  E-6 
-4.00373  E-ll 
6.33019  E-13 
-1.30204  E-6 
6.31961  E-B 
6,5364  E-8 
-2,21358  E-10 
1,49608  Er-B 
2,63304  E-7 
-4,4352  E~9 

1.35588  E-13 
-2.14646  E-15 

1,25906  E-B 
0 
D 
0 
0 
Q 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
e 
0 

6.81961  e-8 
6,5364  E-8 
-2,21358  E-10 
1,49608  E-B 
2,63304  E-7 
-4,4352  E-9 

1.35588  E-13 
-2,14646  E-15 

1,25906  E-8 
-4 , 59427  E-6 
•4.32481  E-6 
1,49608  E-8 
-1,01041  E-6 
-1,74242  E-5 
2,9976  E-7 
-9,16392  E-12 
1,45072  E-13 
-8,36277  E-7 
0 
0 
0 
0 
0 
0 
0 
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2  1  9  6  3  2  1  1 
21563311 
216  51111 

2165  1211 
21691311 
2  1  6  5  2  1  1  1 
2  1  6  5  2  2  1  1 
216523  11 
2  I  6  5  31  1  1 
2  1  6  5  3  2  1  1 
2  1  6  5  3  3  1  1 
21661111 
2  1  ■  ‘6  6  1  2  1  1 
21661311 

2  i  6  6  2  1  1  1 
21662211 
2  1  6  6  2  3  1  1 
2166  3  111 
2  1  6  6  3  2  1  1 
21663311 
22  5  51111 
2  2  5  5  1  2  1  1 
22551311 
22552111 
2  2  5  5  2  2  1  1 

2  2  5  5  2  3 . 1  1 

2  2  5  9  3  1  1  1 
22553211 
2  2  9  5  3  3  t  1 
2256111  1 
22561211 

2  2  5  6 . 1 . 3  1  . 1. 

22562111 
2  2  5  6  2  2  1  1 
2  2  5  6  2  3  1  1 
2  2  5  6  3  1  1  1 
22563211 

2  2  5  6  3  3. . 1 . 1 

2  2  6  5  1  1  1  1 
2  2  6  5  1  2  1  1 
2  2  6  5  1  3  1  1 
2  2  6  5  2  1  1  1 
22692211 

2  2  6  5  2  3 . 1  -1 

22653111 
2265321  1 
22653311 
22661111 
22661211 
2  2  6  6  1  3  1  1 
2  2  6  6  2  1  1  1 
2  2  6  6  2  2  1  1 
2266  2  311 
2  2  6  6  3  1  1  1 
22663211 
2  2  6.6  3  3  1  1 
23551111 
235  5  1211 


0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

Q 

•4,59427  e:-6 
•4,32481  E-6 

1,49608  E-6 . 

•1,01041  E-6 
•1,74242  E-5 
2.9976  E^7 
•9.16392  6-12 
1,45072  E-13 
•  9.36277  E-»7 
•7.36859  E"5 
-3.12943  E-5 
2.63304  E-7 
•1.74242  E-5 
•1,27376  E-4 
5.27565  E-6 
-1» 61281  E-10 
2,5532  E-12 
•7,63647  E-6 


-7,36859  e-5 
-3,12943  E-5 
2,63304  E-7 
•1.74242  E-5 
•1,27376  B-4 
5, 27565  E-6 
-1,61281  E-10 
2,5532  E-12 
-7,53647  E-6_ 
1,3664  E-6 
1,30965  E-6 
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2  3  5  5  1  3  1  1 

2  3952  111 

. ”  2  3  5  5  2  211 

2  5  5  5  2  3  1  1 

23553111 
2  3953211 

2  3  5  5  3  3  1  1 

23941111 
—  2  3  5  6  1  2  1  1 

235  6  1311 

.  2  3  9  6  2  1  I  1 

23562211 
2  39  6  2  3  1  1 

2  3  ■  5  6  3  1  1  1 

2  3  5  6  3  2  1  1 

235633  11 

23651111 
2  3  6  5  1  2  1  1 

2  5  6  5  1  3  1  1 

2  3  6  5  2  1  1  1. 

2  3  6  5  2  2  1  1 

2  5  6  5  2  3  1  1 

2  3653111 

2  3  6  5  3  2  1  1 

236933  11 

2  3  6  6  1  1  1  1 

23661211 
23661311 
2  3  6  6  2  1  1  1 

2  3  6  6  2  2  1  1 

2  3  6  6  2  3  1  1 

__2  _3  .  .6_6__  3  1.  1  -I- 

2  3663211 

_2  3  6  6  3  31  1 

'3  1  9  5  1111 

3  19  5  12  11 

.  3  1  5  5  15  11 

3  1  5  5  2  1  .1  1 

3  1  5  5  2  2  1  1 

3  1  5  5  2  3  1  1 

3  1  5  5  3  1  1  1 

3  1  5  5  3  2  1  1 

"  "  3  1  5  5  3  3  1  1 

3  1  5  _  6, . 1  1  1  1 

3  1  5  6  12  11 

3  15  6  1  3  11 

31562111 
3  1  5  6  2  2  1  1 

319623  1  1 

3  1  5  6  3  1  1  1 

”  3  1  5  6  3  2  1  1 

3  l  5  6  3  3  11 

'  3  16  5  1  1  1  1 

31651211 
S  i  6  5  1  3  1  1 

3  1  6  5.211.1 

3" . 1  6  5  2  2  1  1 

3  X  6  5  2  3  1  1 


"4  I  4352  £**9 
2.?976  E"7 
5i27565  £-6 
-6,98552  E-B 
2,71669  £-12 
-4,30071  £-14 
2,5227  E-7 
0 
0 
0 
0 
0 
0 

0  .  , 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

D 

1,3664  e-6 
1,30965  E-6 
-4,4352  £-9 
2,9976  £-7 
5,27565  £-6 
-5,98652  E-8 
2.71669  £-12 
-4,30071  £-14 
2.5227  £-7 
,4.1772  E-11 
-4.00373  £-11 
1,35508  £-13 
-9.16392  £-12 
-1,61281  £-10 
2.71669  £-12 
-8,30515  £-17 
1,31477  £-18 
-7.71212  E-12 

0  . . . . . 

0 

0 

0 

0 

Q 

B 

0 

0 

0 

0 

0 

0 

0 

0 
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3  16  9 

3  16  3 

3  16  6 

3  16  6 

3  16  6 

3  16  6 

3  16  6 

3  16  6 

3  1  6  6 

3  16  6 

3  16  6 

3  2  5  5 

3  2  ■  5  5 

3  2  5  5 

3  2  5  5 

3  2  5  5 

3  2  5  5 

3  2  5  5 

3  2  5  5 

3  2  5  5 

3  2  5  6 

3  2  5  6 

3  2  5  6 

3  2  5  6 

3  2  5  6 

3  2  5  6 

3  2  5  6 

3  2  5  6 

3  2  5  6 

3  2  6  5 

3  2  6  5 

3  2  6  5 

3  .2  6  5 

3  2  6  5 

3  2  6  5 

3  2  6  5 

3  2  .6.5 

3  2  6  5 
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